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ABSTRACT 

This study treats the longitudinal resonance of waves and 
energetic electrons in the eartn's magnetosphere, and the possible role 
this resonance may play in generating various magneto spheric phenomena. 
The first part of the study Is concerned id.th the derivation of 
time-averaged nonlinear equations of motion for energetic particles 
lon'^ltudlnally resonant with a whistler mode wave propagating with 
non-:;ero wave normal. It Is shown that the wave magnetic forces can be 
neglected at lower particle pitch angles, while they become equal to or 
larger than the wave electric forces for a>30°. The time-averaged 
equations of motion were used In test particle simulations which were 
done for a wide range of wave amplitudes, wave-normals, particle pitch 
angles, particle parallel velocities, and in an Inhomogeneous medium 
such as the magnetosphere. It was found that there are two classes of 
particles, trapped and untrapped, and that the scattering and energy 
exclv^nge for those two groups exhibit significantly different behavior. 
The trapped particles are characterized by a bounded ptiase variation 
Cwlth respect to the wave) which Is less then 2ir, whereas the phase 
variation of untrapped particles la unbounded. It Is also found that 
the trapping of the particles requires that the wave amplitude exceed a 
certain threshold value, and that the trapped electrons become space 
bunched due to the interaction. The full distribution simulations 
Indicate that the expected particle precipitation Is considerably 
smaller (one order of magnitude) compared to gyroresonance-induced 
precipitation for w^ves of comparable amplitude, which shows that the 
scattering efficiency of the longitudinal resonance Is small. The 
amplitude threshold effect, together with the space bunching effect, was 
found to support one of the mechanisms suggested to explain whistler 
precursors. 
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I, INTRODUCTION 


A. ORGANIZATION OF MATERIAL 

This study treats longitudinal resonance Interactions between 
energetic electrons and VLF waves In the earth's magnetosphere. The aim 
was to derive suitable analytical methods for test particle studies, and 
then to use those methods to Investigate various aspects of the 
longitudinal resonance process. 

The first part of the study Is concerned with the derivation of 
equations of motion and their applications to the longitudinal resonance 
for a wide range of magneto spheric parameters. The second part gives 
the results of the numerical simulation of wave-particle Interactions. 
The numerical simulations are done using a test particle approach to 
determine the perturbations of pitch angle for various wave functions. 
Also Investigated are the perturbations of the full particle 
distribution and the energy exchange process. 

In conclusion the longitudinal resonance Interaction Is compared 
to the cyclotron resonance Interaction, and Is related to phenomena 
observed In the magnetosphere. 
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B. WAVE-PARTICLE INTERACTIONS IN THE MAGNETOSPHERE 


The magnetosphere, a magnetized region extending from about 1000 
km altitude up to distance of roughly 100,000 km from the earth, Is 
filled with both 'cold* and 'hot' plasma; the cold plasma consists of 
electrons and protons with energies In the 0.1-1 eV range, while the hot 
plasma consists of energetic particles with higher energies In the range 
from 100 eV to tens of MeV. The cold plasma together with the earth's 
static magnetic field determines the wave propagation properties of the 
magnetosphere. The hot plasma Is a source of energetic particles which 
participate In the wave-particle Interactions that result In radio wave 
emissions. As seen from both ground and satellite observations the 
magnetosphere supports numerous modes of wave propagation. It can be 
shown that the hot plasma, due to Its very low density, does not affect 
Che wave dispersion properties of Che magnetosphere, l.e. the dispersion 
of waves can be explained assuming Chat only cold plasma Is present. 

It Is known that very-low-frequency waves can propagate In the 
magnetosphere with phase velocities much smaller Chan the velocity of 
light, and chat chose waves, called whlstler-mode waves, can undergo 
interactions with energetic particles both through longitudinal 
resonance and cyclotron (gyro) resonance. In longitudinal resonance the 
particle parallel velocity is matched to the wave phase velocity, 
whereas In the cyclotron resonance the doppler-shlfted frequency of Che 
wave (shifted due to the particle parallel velocity) matches Che 
gyrofrequency of the energetic particle. Both types of Interactions may 
Induce perturbations of the energetic particle distribution through 
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pitch angle scattering, and may also result In different types of radio 
wave emissions, wave amplification (growth) and wave attenuation. The 
purpose of this study Is to Investigate the longitudinal resonance 
Interactions of energetic particles with whistler mode signals 
propagating at an oblique angle to the static magnetic field. The 
approach taken Is to use a test particle analysis and to study how the 
resonance process depends on various parameters. The particle 
trajectories are then used to estimate other effects such as wave 
growth/damping and particle trapping and precipitation. The trajectory 
calculations were done using a set of nonlinear equations of motion 
which are averaged over one gyroperlod [inan and Tkalcevlc, 1982]. 


C. PREVIOUS WORK ON LONGITUDINAL RESONANCE 


The longitudinal resonance process has been Invoked by many 
authors to explain various magnetospherlc wave phenomena. One of the 
early works considered the travellng-wave-tube type of process as a 
generation mechanism for VLF emissions [Gallet and Helllwell, 1939], and 
this process was also considered for amplification of whistler mode 
signals fBrlce, 1961 ]. The travellng-wave~tube mechanism was also 
considered by Dowden [l962] as a possible mechanism of hiss generation. 
Bell [l964.1 derived linearized solutions for the trajectories of 
longitudinally resonant particles, but these have not been extended to 
cover the nonlinear regime. The various emission-generation theories 
have been reviewed by Brice [l964], Including both Cerenkov radiation 
and the traveling wave amplification hypothesis. The Cerenkov mechanism 
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Is a process In which charged particles radiate electromagnetic waves as 
they travel through a medium. The necessary condition for the existence 
of this type of radiation, called a coherence condition, is easily 
found, and is the same as the condition required for the longitudinal 
interaction between the wave and the particle. Therefore, it is evident 
that those two processes, the longitudinal resonance interacrions and 
Cerenkov radiation, are based on the same physical principle. 

The Cerenkov radiation mechanism has been suggested by many 
authors [eIIIs, 1959,1960; Dowden, 1960; McKenzie, 1963] in order to 
explain VLF hiss. The problem of stability of whistler mode signals, 

. . .. the possibility of wave growth, accounting both for longitudinal 
and gyroresonance effects, was discussed by Kennel and Fetschek [l966]. 
Kennel and Thorne [l967], and also by Brlnca [l972]. The work on 
radiation from moving charged particles, which Includes Cerenkov 
radiation. Includes the analysis done by Llemohn [l965], Mansfield [l967], 
and Seshadrl [l967]. A good review of work done on Cerenkov radiation, 
along with additional analysis of the hiss power density spectrum, was 
given by Taylor and Shawhan [l974]. Their work gives examples of the 
power spectral density of hiss, both measured [curnett, 1966; Gurnett 
and Frank, 1972], and calculated [jorgensen, 1968; Lim and Laaspere, 

1972]. Swift and Kan [l975] showed that an electron beam can excite a 
whistler mode instability near the resonance cone through the 
longitudinal resonance interaction. Maggs [1976] and Kumagal at al. 

[l980] investigated beam amplification due to Cerenkov radiation from 
longitudinally resonant electrons, and considered this type of beam 
instability as a generating mechanism of VLF hiss. The whistler 


5 


precursor generation mechanism of Park and Helllwell [l977] was based on 
modlf tcaClons of the particle distribution function achieved through 
longitudinal Interaction between whistlers and energetic electrons. 

Most of the above studies were primarily concerned with wave 
growth calculations using the wave dispersion relation. On the other 
hand, the detailed nonlinear motion of longitudinally resonant particles 
was studied only for the case of electrostatic waves [Nunn, 1971; 1973]. 
Palmadesso [l973] derived equations of motion for a case of oblique 
propagation, and used particle trajectories to estimate thc>. nonlinear 
time dependent Landau damping rate of the wave. 


D. CONTRIBUTIONS OF THE PRESENT W)RK 


The motion of electrons longitudinally resonant with a whistler 
mode wsve propagating at an angle to the static magnetic field Is 
represented by a simple set of equations motion which are averaged over 
the cyclotron period. It Is shown that these nonlinear equations are a 
very accurate representation of the electron motion for a wide range of 
magnetospheric parameters. 

Using the time-averaged nonlinear equations of motion In 
numerical simulations involving whistler mode signals propagating In an 
Inhomogeneous medium It wrs found that the effects of wave magnetic 
forces can be neglected for low pitch angles, high wave norm&l angle, 
and/or high normalized wave frequency. At the higher pitch angles the 
wave magnetic forces become very important and it Is necessary to 


6 


Include Che additional force terms as derived. 

The sample calculations indicate that there are two classes of 
electrons, distinguished by the behavior of their phaces with respect to 
the wave. In a case when Che phase variation is bounded, l.e. less then 
2ir, Che electron is said to be crapped, «diereas unbounded phase 
variation characterizes the uncrapped electrons. The scattering and 
corresponding energy exchange for Che trapped and untrapped electrons 
exhibit significantly different characteristics. 

It is also found chat the crapping of electrons is easier under 
conditions of spatial amplitude variation of a na;rowband signal rather 
than for a constant amplitude. Analysis was done for a constant 
amplitude CW signal, a CW signal amplified at Che equate »r through 
gyroresonance, and also for a spatial amplitude variation of the pulse 
formed by a nonducted signal. 

It is also shown that the longitudinal resonance process 
Involves a wave amplitude threshold effect, i.e. the crapping of 
electrons is possible only if the amplitude of the wave parallel 
electric field Eh exceeds a certain value. The trapped electrons also 
become space bunched and temporarily Increase the electron density over 
a particular range of parallel velocities. 

The full distribution results show that the expected 
precipitation is small when compared to gyroresonance-lnduced 
precipitation for waves of comparable amplitude. In general, the 
results indicate that the longitudinal resonance scattering efficiency 
(scattering vs. amplitude) is considerably smaller, i.e. the 
efficiencies of the two processes differ by as much as an order of 
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magnitude. 

The amplitude threshold effect was tested on whistler 
precursors, and It was found that the whistler amplitudes are well 
correlated with the occurrence of precursors, l.e. only whistlers with 
amplitudes above a certain threshold resulted in precursors. This 
provides support for the whistler precursor generation mechanism 
suggested by Park and Helllwell [l977], which Involves longitudinal 
resonance Interactions, and therefore It should exhibit a threshold 
effect as Indicated by the measurements. 



II. BASIC PHYSICS AND TIME AVERAGED EQUATIONS OF MOTION 


A. MOTION OF CHARGED PARTICLES IN EARTH’S MAGNETIC FIELD 

Motion of the charged energetic particles in the magnetosphere 
Is governed by the earth's magnetic field. The earth's field In the 
Inner magnetosphere can be approximated by the dipole model with the 
magnetic field s«*rength given as 

Bq -0.312-10 '*(Rq/R)^( 1 + 3sin^x/^* Wb/m^ (2.1) 

where X Is the geomagnetic latitude » R Is geocentric radius » and R^ Is 
the radius of the earth. The axis of the magnetic dipole Is Inclined 
with respect to the rotation axis by 11^. 

The motion of a particle In the magnetosphere Is uniquely 
described by either the parallel and perpendicular velocities of the 
particle, v„ and v^ respectively, or by the parallel (perpendicular) 
velocity and pitch angle a ■ arctan( Vj, / v„ ). Fig. 2.1 shows a 
typical geometry with the definitions of v„ , Vj, , and a . 

It can be shown that for a spatially changing magnetic field, 
such as the earth's magnetic field given by Eq.2.1, charged particles 
will bounce forth and back along the field line between the mirror 
points [Northrop, 1963; Buneman 1980]. This Is so because the particle 
perpendicular velocity must change In order to satisfy adiabatic 
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Invariants > while the total kinetic energy of the particle must remain 
constant. The first adiabatic Invariant is the invariance of the 
orbital magnetic moment, given as 

Wj,/B * constant (2.2) 

where Is the perpendicular kinetic energy of the particle. 



FIGURE 2.1 DIPOLE GEOMETRY AND SYMBOLS USED FOR PARTICLE IDENTIFICATION. 

Note that the z-axls is aligned with the magnetic field line 
and that both the wave and the particles travel in the + 2 , direction. 
Particle orbits are described In terms of equatorial values of Vu and a. 
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The second adiabatic Invariant requires that the magnetic flux 
through the circle described by the particle gyrating around the field 
line remains constant , or 


r„ X B ■ constant 
n 


(2.3) 


where r_ Is the electron gyroradlus. 
n 

Thus if the magnetic field increases, the perpendicular 
kinetic energy must also Increase according to Eq. 2.2. Furthermore, 
the parallel energy W„ of the particle must decrease so that the total 
energy W., 4- Wx remains constant. Therefore, the particle pitch angle 
a“arctan ^ increases as B Increases up to the point where 

a ■ 90®. At this point the parallel velocity of the particle has been 
reduced to zero, and the particle begins to travel in the opposite 
direction along the same field line. When the particle reaches the 
conjugate point where again a 90®, the process repeats. Hence the 
particle bounces back and forth along the magnetic field line between 
the two mirror points where v„ ■ 0. 

Finally, the motion of a particle trapped along a field line can 
be described by the following equations 


2 

dvii ^ _ "^x ^ dBo 
dt " ” 2Bo dz 


(2.4) 


dvx , ^ vii va , dBo 
dt 2Bo dz 


(2.5) 


which can be derived from the first adiabatic Invariant and the law of 


n 


energy conservation. 


B. LONGITUDINAL RESONANCE 


The bounce motion of the particles can be affected by resonant 
Interactions between waves and the particles. The resonance condition 
Is satisfied whenever the doppler-shlfted frequency of the wave seen by 
the particle Is equal to an Integral multiple of the particle 
gyrofrequency, l.e. 

0 ) - h„v„- moj m “ 0,±1,±2,±3,... (2.6) 

H 

where cj Is the wave frequency » k„ Is the wave number In the direction of 

the static magnetic field, and o) is the particle gyrofrequency. 

H 

The resonance condition given by Eq. 2.6 can be furth divided 
into three subgroups according to different values of the parameter m. 
For m>0 we have the resonance condition for the mr-th order 
gyroresonance; m<0 is the resonance condition for the m-th order 
anomalous gyroresonance; m~0 yields the resonance condition for the 
longitudinal or Landau resonance. The last condition is given as 

0 ) - k„ v„- 0 (2.7) 

or 

Vp«- v„ 


( 2 . 8 ) 


12 


where Is the wave phase velocity measured In the direction of the 
static magnetic field. 

Before discussing the longitudinal resonance we should note that 
this resonance (m*0) is fully separable from the gyroresonances (m/0) . 
since the longitudinal resonance is possible only when the wave and the 
particles travel in the same direction, while the gyroresonance 
condition is satisfied only if the wave and the particles travel in the 
opposite direction. This separability of the different resonances makes 
their analysis much simpler. It is still possible for the same particle 
to Interact simultaneously in both resonances with two different waves 
that satisfy corresponding resonant conditions. In this report we shall 
limit ourselves to discussion of the longitudinal resonance, although a 
comparison with the gyroresonance mechanism is given later in the text. 

The condition given in Eq. 2.8 is the necessary condition for 
the longitudinal resonance. However, in order for the particle and the 
wave to exchange energy through the particle trapping process, the 
parallel component of the wave electric field must have a non-zero 
value. Therefore, even if the particle parallel velocity matches the 
wave phase velocity there will be no energy exchange between the 
particle and the wave if E., - 0. The direction of the energy exchange 
(whether wave or particle gains energy) depends on the initial velocity 
of the particle v„ . In the case when v„ is initially less than the 
phase velocity v^^^ the particle will gain energy; if the initial v,i is 
larger than v^^^ the particle will lose some of its energy. We shall now 
present a simple analytical model for the longitudinal resonance and 
trapping process similar to that given by Seshadri [1973]. 
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Let us assume that the longitudinal component of the wave 
electric field, propagating In the homogeneous medium. Is given by 


E,i (s,t) " E,^sln(k„«s - o)»t) 


(2.9) 


where s Is the space coordinate. Eq. 2.9 Is written In the laboratory 

coordinate system, but It Is useful to do the analysis In the wave frame 

which moves at the phase velocity v . In this case a new space 

P" 

coordinate z Is defined as 


z 


8 - Vp„t 


( 2 . 10 ) 


Now, Eq. 2.9 can be rewritten as 


E„ (s,t) - E„ sin [k„ (s - ^ t)] (2.11) 

O Kii 


and using Eq, 2.10 and v_ - 

“ " K|i 


Eq. 2.11 simplifies to 


E„(z) - E,,^ sln(k„ z) 


( 2 . 12 ) 


The electric field given by the Eq. 2,12 Is static In the wave 
frame and It is possible to derive a corresponding scalar potential 
4(z),by Integrating En(ri) where n is a dummy variable. 

z 

<Kz) - - (E(n) • dn (2.13) 


o 
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z 

<t(z) ■ - |E„^sin(k„n) dn (2.13a) 

o 

E«o 

«(z) ■ T — (coa(k„z) - 1) (2.13b) 

*V|| 

Next we consider an electron (a similar derivation Is possible for other 
types of charged particles) and Its potential energy Wp(z) which. In the 
wave frame Is given by 

Wp(z) ■ - e*<Kz) (2.14) 

W (z) • (1 - cos(k„z)) - (l-cos(k„z)) (2.14a) 

P Pmax 

The constant of the Integration is chosen such that the minimum 
potential energy given by Eq. 2.14a Is zero. Thus, the potential energy 
of the Electron Is a periodic function, as shown In Fig. 2.2. 

It can be shown that the possibility of an electron being 
trapped depends on the Initial kinetic energy of that electron measured 
in the wave frame. In a case when the initial kinetic energy of an 
electron, placed at z at the time t"0. Is larger than the potential 
energy given by Eq. 2.14a, , there is no net Interaction between 

the wave and electron, regardless of the electron Initial velocity. The 
electron simply slides up and down the potential well as it moves either 
forward or backward through the wave, and there is no net energy 
exchange when averaged over one wavelength. 

However, if the kinetic energy of the electron In the wave 
frame, Wj^(t*0), Is less than the potential energy given by Eq. 2.14a, 
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Wn(z) 




7T/2 t=0 


FIGURE 2.2 PARALLEL ELECTRIC FIELD AND THE CORRESPONDING POTENTIAL 
ENERGY. Both the parallel electric field E„ and poten- 
tial energy Wp of the electron are periodic functions in a reference 
frame moving at at the parallel phase velocity Vp„. In (b) zg indi- 
cates the bottom of the potential well. 
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^Pmax ^^8* 2.2 the electron is trapped in the potential 

well. The trapping condition is then given as 


2 

(Vii-Vp„)^ < w 

(2.15) 


^max 


1 

2 “ 


(2.15a) 

iv„ 

- mk„ 

(2.15b) 


Rewriting the Inequality of Eq. 2.15b as 


J 2eEiiQ 
nk„ 


2eE 


"O 


< V.. < v_ + / • i 

Pm mk„ 


(2.16) 


we have a range of velocities for which It Is possible to trap an 
electron. Therefore, all electrons with parallel velocities that 
satisfy Eq. 2.16 are trapped In the wave potential well. The trapping 
velocity bandwidth v^ Is given as 


▼t 



(2.17) 


Furthermore, It can be shown that the total energy, exchanged 

between the wave and electrons during the trapping process Is 

Vp..+Vt 

AW - lf(v„)AEdv„ (2.18) 

Vp„ - Vt 

where f(v„) is the electron distribution function; AE is the amount of 
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energy exchanged through trapping of a single electron, and It Is 
expressed as 


^ ■ 



7 “ 


(2.19) 


AE 


(v - v_ ) 

6 p II •• P II 


(2.19a) 


wnere On Is a time-varying periodic function describing the oscillation 
of an electron at the bottom of the potential well. Expanding f(v„) in 
a Taylor series around v„ ■ Vp^^ we obtain 


£(v„ ) - f (Vp j + (v„ - VpJ 


3f(v„) 


( 2 . 20 ) 


Pi. 


and finally substituting Eq. 2.20 In Eq. 2.18 the total energy exchanged 
In the trapping process, AW, Is given as 


AW - -|-m V v^^ f (■'fc) (2.21) 

3 Pii c p„ 

The result derived In Eq. 2.21 shows that the net energy 
exchanged between the trapped electrons and the wave depends on the 
slope of the distribution function at a point where the electron 
velocity Is equal to the phase velocity of the wave. In the case when 
the number of electrons moving faster Is larger than the number of 
electrons moving slower than the phase velocity, the wave gains energy 
and Its amplitude grows. Similarly, If the number of slow electrons Is 
larger than the number of fast electrons, the amplitude of the wave Is 
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reduced. 

The above analysis, using a longitudinal plasma wave and 
one-dlmenslonal distribution function f(v„), has demonstrated that It Is 
possible to have wave damping In the absence of collisions, also known 
as Landau damping. It was also shown that the wave amplitude grows If 
the slope of the distribution function Is positive. However, the 
expressions for the energy exchange were derived assuming that the 
particles are already trapped. It was also assumed that the medium is 
homogeneous, and that both the wave and the distribution function are 
one-dlmenslonal . 

In the magnetosphere Eq. 2.18 Is still valid, but the trapping 
process Is governed by the particle equations of motion. Thus In order 
to find the energy exchanged between a wave and particle ( AE In Eq. 
2.18) It Is necessary to derive the equations of motion for a single 
particle when It Is in longitudinal resonance with waves In the mag- 
netosphere. 


C. NONLINEAR EQUATIONS OF MOTION FOR LANDAU RESONANCE 
INTERACTIONS t/ITH A WHISTLER MODE WAVE 

Now we consider an elllptlcally polarized wave propagating In 
the cold plasma of the magnetosphere with a static magnetic field 
The wave frequency f Is assumed to be less than the electron 
gyrofrequency fg ; In that case there Is only one propagating wave 
[Ratcllffe, 1959; Budden, 1961], which Is called a whistler wave. 




ORIGINAL PAGE iS 
OF POOR QUALITY 


In the most general case all Cartesian components of the wave 
electric and magnetic field have non-zero values. All of these 
components can be expressed In terms of through the cold-plasmi 
dispersion relation. Without any loss of generality the wave vector k 
Is confined to the x-z plane, at an angle 6 from the sL ^tlc magnetic 
field. The coordinate system used Is shown In Fig. 2.3. 



FIGURE 2.3 COORDINATE SYSTEM FOR THE EQUATIONS OF MOTION. The wave 

vector E* Is at an angle 0 from the static magnetic field 7^. 
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Me also assume propagation as exp l(u) t - k«r). Using a plasma 
dispersion relation [Stlx, 1962] 


Ej. - n^cos^S - lEjj n^sln6 cos0 




Ej. - n‘ 


n*slu9 COS0 0 E,i - n^sln^0 


Sy 


0 ( 2 . 22 ) 


all electric field components can be expressed In terms of 5 as follows 


fig ■ E„ cos( tut - ^*r) 


(2.23) 


« . g' ; E„ cos ( (Ut - k*r) 

n^sln0cos0 


(2.24) 


- E~ n^sln0-Eii „ ^ ^ 

Sy - E„ sin (tut - k*r) 

n*-Ej.n^sln0cos0 


(2.25) 


where e„« 1 - Ex“ 1 - — » £»■ ~ — • The refractive Index 

(U^ 


2 2 
-“h 


2 2 
(U U) - tUg 


n can be derived from Eq. 2.22 as (QL approximation) 

^2 


n^« 1 + 


f(f{jCOS0 - f ) 


(2.25a) 


3 B 


Using Maxwell's equation 7 x E the wave magnetic components are 


a _ kcos0 u 
*x " ^ ®y 


(2.26) 
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»y - 


k cos 6 
w 


«x - 


k sin 6 

OJ 




(2.27) 


3* - 


k sln9 

b) 




(2.28) 


which can be also expressed In terms of using Eqs. 2.23, 2.24, and 
2.25. 

The variation of th% total electron velocity v is governed by 
the Lorentz force equation 


m ~ - q C^+ V X (B^ + Bq )] 


(2.29) 


where m and q are electron mass and charge. For the case when 
|BwI« ! I * electron gyromotlon can be assumed to be unaffected 

by the wave to the first order, so that the Cartesian components of the 
electron velocity vary as 


vz - v„ 


(2.30) 


- V* cos (oigt + Sq) 


(2.31) 


■ Vx sin (Wgt + 3^) 


(2.32) 


where u>{{ is the electron gyrofrequency and Is the initial cyclotron 
phase. Furthermore, as long as the wave field is much sn<aller chan the 
earth’s magnetic field, it is permissible first to derive the force 
applied to an electron by the wave fields and then to superimpose the 
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adiabatic variation of v^ and v,i. Therefore, the perturbation of the 
electron motion induced by the wave fields only is given by 


dv — — 

“ dt " ^ * »w] 


(2.33) 


It is useful to examine each Cartesian component in Eq. 2.33 separately. 
These three components are given as 


“ *lC®x ^y®2~ ^z®y] 
Fy - q[Sy + VgSx- V^Sj] 
Fz - q[Sz + v^Sy- VySx] 


(2.34) 


(2.35) 


(2.36) 


Before investigating those equations we simplify cos( o)t - k*r), which 
can be expressed as 


cos( wt - k COS0* z - k sin9» x) 


(2.37) 


or letting y " cat - k cos9 z in Eq. 2.37 we have 


cos( Y ~ k sin9 x) 


(2.38) 


Eq. 2.38 can be further simplified using the fact that 


Wjj 


X 


sin(coyt + 6 q) 


(2.39) 
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which is derived by integrating Eq. 2.31. Finally, replacing x in 
Eq. 2.38 by (2.39) 

cos( (lit - k«r) ■ cos(y - n sin<j) ) (2,40) 

Vxk sin6 

where d» - (Out + 6 and n ■ 77 

Now, using the result derived in (2.40) we can rewrite three 
Cartesian components of Che Lorentz force as 

9 [ExSin(Y - H sin<|> ) + v^, sin4> By sin(Y - H sin<<> ) 

- v„ cos(y - n sin(|> )] (2.41) 

Fy» q [EySin(Y - H sin(J> ) + v„ B^ sin(Y - 0 sln(j> ) 

- VxCOs(J> Bj sin(Y - 0 sin<|> )] (2.42) 

Fj- q [EgCOs(Y - H sin4» ) + Vx cos<^ By cos(Y - 0 sin0 ) 

- Vj^^sln^ Bx 3in(Y - H sin<j) )] (2.43) 

Note that Exf By » Ez , Bx » By , and Bz are the real 

magnitudes of the fields, with the phase differences taken separately 

slzi 

into account through ( Y 0 sin<|) ) terms. 

COS 

At this point we have three equations which can be used to 
describe the motion of particles in resonance with a whistler wave. 
However, it is desirable to reduce the number of required equations to 
simplify numerical simulations. In this case it is useful to combine 
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Che X and y couponencs of the Lorentz force In one perpendicular 
component. This Is done by taking the time derivative of the square 
of the perpendicular velocity vj ■ v^ + v^ 


vj - v^ + /4r 

•*- X y 'dt 


(2.44) 


dv. dVv UVv 

Vi— ~ ■ V ^ + V ^ 

^dt xdt ydt 


(2.44a) 


and multiplying it by m/v^ 


+ m is 

dt Vx dt v^ dt 


(2,45) 


v« Vv dvi 

However, ■ cos (J) , - sin d) , m -7-=- 

Vx V. ^ dt 

dVy 

m - Fy, and (2.45) reduces to 




and 


Fx ■ cosij) Fjj + sln(|i Fy 


(2.46) 


Now, combining Eqs. 2.46, 2.41, and 2.42 the perpendicular force term is 


Fx“ coscj) {q [ExSin(Y -risin(() ) + Vj^sincf) By sin(y -nsincj) ) 

- v„ Bz cos(y -nsiruj) )]} 

+ sin(<) {q [Eysin(Y -r)sln(J) ) + v„ B^ sin(Y -Tisin(|) ) 

- coscj) Bz sin(Y -Hsincj) )]} 


(2.47) 


The motion of a particle is now described in terms of the 
parallel and perpendicular forces, given respectively by Eqs. 2.43 and 
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s 

2.47. If the (Y~ ns In (^) terns In these equations are expanded 
cos 

(Appendix A) , the result Is an Infinite series of harmonics at 
frequencies ncu^ with amplitudes given by In a general 

formulation all terms must be kept and Eqs. 2.43 and 2.47 must be used 
as they stand. However, the equations can be considerably simplified 
when time averaged over one cyclotron period, Tg , because the higher 
order force terms (n s 2) vanish. Also, qualitatively, the v^Sy term 
should average out to zero since wave phase does not vary In the 
y<-dlrectlon. In the next section we present the necessary conditions 
for the averaging to be valid, along with the time averaged equations of 
motion. 


D. TIME AVERAGING OF EQUATIONS OF MOTION 

Before averaging Eqs. 2.43 and 2.47 over one gyroperlod we have 
to make sure that the wave phase variations, as seen by the particles 
during one gyroperlod, are negligible. For the small field case this 
condition can be stated as 

0 ) - k»v « o)g (2.48) 

which would certainly be the case for the Landau resonance described by 


U) - k«v 25 0 


(2.49) 
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Note Chat Eq. 2,49 Is Che equivalent of Eq. 2.8. 

Ue have stated condition (2.48) assuming small amplitude waves. 
This requires Chat the wave field be small enough chat It cannot move 
the particle by a substantial fraction of a wavelength during a 
gyroperlod. This condition can be stated as 

« ij <2-50) 

where Sp Is Che peak parallel acceleration » c Is the speed of light, n 

Is Che refractive index, ^ ^ wave frequency and fy ■ ~ is 

the electron gyrofrequency. The peak value of the parallel acceleration 

3tt 

ap during a gyroperlod can be taken to be that for (|> * 

7T 

Y - n sincj) ■ 2 . From Eq. 2.41 we have 

Kl - IS ( Ez-VyB^)l (2.51) 

In a order to express E 2 In terms of we have from Eq. 2.25 

Sy - Pj «2 (2.52) 

where p ■ ±~^ ^ ^ Substituting Eq. 2.52 in Eq. 2.26 

* n*-£ n^sinOcose Sj 

B Lcos0_ (2.53) 

X m f'z z 


or 
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Bx ^ 

Pz kcos6 


(2.54) 


Furthermore, for the near resonant particles ■ v_„ ■ v„ and 

kcosD P 

Eq. 2.54 yields 


Ez---^v„ (2.55) 

Pz 

Replacing the £2 in Eq. 2.51 with the above expression the peak 
acceleration | ap | Is 


where tan a ■ — . 

V,i 

The final step Is to substitute (2.56a) In (2.50) In order to 
get the condition on wave Intensity for which the averaging of equations 
(2.43) and (2.47) Is valid; 


(2.56) 

(2.56a) 


®x « ®u 


mfp c __ Pz tana 
qv^nf 1+ 1 P 2 1 tana 


(2.57) 


Thus Bu represents the upper limit on wave magnetic field intensity. 
Note that E^c Is equal to the total transverse for circularly 
polarized whistler waves. Assuming to have a value much higher 
( > 100 times ) than the typical field Intensities for whistler mode 
waves In the magnetosphere [Burtls and Helllwell, 1975], as shown later 
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In the text, we shall now time average Eqs. 2.43 and 2.47 over one 
gyroperlod. In doing so we use the Identities derived in Appendix A. 
The averaged equations of motion become 

dvz 

< m ■“ — > ■ < q iz >- < q vy 2x ^ (2.58) 

dt 

< m 4"^ > - < q 5y >- < q vz 8x > (2.59) 

n ^ • 


or 


dvii r vxcose Jx(n) T 

m q E^o(n) Ll - — Pz J sln( ait - k z cos6 ) (2.60) 

dvx r Vm kcos6 -i 

m TT"" - qPz Ez Jl(ri) Ll - 7 ; sln( ojt - k z cos0 ) (2.61) 

Qu CU 

dv„ dVx 

Since the brackets on the left hand sides are dropped, and 
should be understood to be the average rates of change of Vm and Vx , 
respectively. 

Finally, for an Inhomogeneous medium with Bq variable as In the 
magnetosphere, the adiabatic variations of Vn and Vx can be superposed 
on Che wave-induced perturbations as long as the variation of Bq In one 
wavelength Is negligible. Thus Che complete averaged nonlinear 
equations of motion become 


dvi, 

dt 


q _ . . V r- v^kcosB Ji(ri)T ... ..V Vx dB^ 

■^EjJo(n) [1 Pj-j^.lnCojC-kzcose ) - 


(2.62) 
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v„ kcos9 
0 ) 


■]sin(u)t-kzcos9 ) + 


Vii Vi dBp 
2Bq dz 


(2.63) 


We shall discuss 'Che relative importance of the different terms 
in Eqs. 2.62 and 2.63 in the next section. 


E. DISCUSSION OF FORCE EQUATIONS 


Two terms of the parallel force are: 


<q Sg> - q EjJpCn) siny 


<q Vy Sx> “ “ 9 2gJx(n) patan a sin y 


Also note that using (2.A9) 


k sin# to a k cosd 

ai a « tano — Vi 


0)H <*»H 


to 


n ■ — tan 9 tan oc 

WH 


for near-resonant particles. 


(2.64) 

(2.65) 


( 2 , 66 ) 


(2.6ba) 


The term in (2.64) proportional to qE 2 Jo(n) is similar to the 
qE^ term that would be prefsnt in the case of electrostatic waves. The 
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•Jo(n) represents the fact that the Ez field seen by the particle at 
different points In its transverse orbit is changing since has a 
transverse phase variation given by k x sin 0 . The term in (2.65) 
represents the effect of the q v x B force, and the fact that since the 
plane of rotation of the particle and the wave polarization ellipse are 
at an angle ( j - 6), there is a net longitudinal acceleration even 
after averaging over one gyroperlod. For cases in which (2.64) is the 
dominant term, the equations of motion for Interaction with whistler 
mode waves are much the same as those for electrostatic waves [Nunn 
1971, 1973]. 


Before comparing the relative magnitudes of (2.64) and (2.65) 

for the range of the parameters In the magnetosphere It should be noted 
dv^ 

that , given by Eq. 2.61, becomes very small for near-resonant 

.... ' V:ikcos9 v„ 

particles with Vi, ■ Vp„ . In this case 1 1 k 0, 


tii 


'Pn 


and the perpendicular motion of the particles is primarily governed by 
the adiabatic term of Eq. 2.63. In the following figures we present the 
magnitudes of (2.64) and (2.65), as well as the longitudinal 
polarization P 2 as a function of different parameters. 

Figure 2.4 shows a plot of the longitudinal polarization P 2 as a 
function of the wave normal angle 0 , for different values of normalized 

(I) 

frequency . The results are computed by using the cold plasma 

lAIJJ 

dispersion relation [Stlx, 1962]. The longitudinal polarization Is 
Pz ^ , as defined In (2.52). A plasma frequency fp ■ 180 kHz, 
corresponding to 400 el/cC at the magnetic equator at L ■ 4, along with 
the equatorial gyrofrequency fa ■ 13.65 kHz, were used In computing P 2 » 
For fp » fa the value of p 2 is not strongly dependent on fp. Note from 
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6 (degrees) 


FIGURE 2.4 MAGNITUDE OF THE WAVE LONGITUDINAL POLAR- 
IZATION (pzl - AS A FUNCTION OF 

WAVE NORMAL ANGLE 0. \ qz \ shown for three different 
normalized frequencies. 


f 
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Fig. 2.4 Chat P 2 Is In general higher at lower frequencies and decreases 
with Increasing 9 . Also recall that for longitudinal propagation, 
l.e. , 9 ■ 0*’, Ez " 0 and there Is no Interaction between the particles 
and the waves. 

Figures 2.5, 2.6 and 2.7 compare the peak magnitudes of the two 
terms as given by (2.64) and (2.65) for various parameters. Figure 2.5 
shows variation of both terms with pitch angle a , for various wave 
normal angles 3 and f ■ 0.5 fn. Ic can be seen that the ^qvy3x> term Is 
negligible for lower pitch angles, while it becomes equal to or larger 
than the <q5z> term for ot>30*. As long a<30®, the <q5z> term alone can 
be used to compute the motion of the Landau resonant particles with less 
than 10% error. 

Figure 2.6 shows the dependence on the wave normal angle for 
various pitch angles a and for f * 0.5 fg. The resonance cone angle for 
this frequency is -60° as shown. This result indicates that for any 
pitch angle a , the <qVyS^> term is more important at lower wave normal 
angles, but that there is a strong dependence on pitch angle as was also 
indicated in Figure 2.5. For 9 approaching zero Ji(n) goes to zero and 
Pz approaches infinity. As a result, the<qvy3]2> term will go to zero 

A 

and may be approximated by -qE 2 sinYtan a(l - f/f|j)/(2 + 2f/fjj) for small 
values of 9 (Appendix A). 

Finally, Figure 2.7 shows the variation of the terms with 
normalized frequency f/fg* The curves are for a ■ 40° and three 
different values of wave normal angle 9. It can be seen that the 
magnetic field term is more Important at lower frequencies, although the 
dependence on frequency is not as strong as that on 9 and a . 
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FIGURE 2.5 NORMALIZED PEAK MAGNITUDES OF THE <qvySx> 
AND <q«2> TERMS AS FUNCTIONS OF PITCH AN- 
GLE a. The results shown are for f « 0.5 and for 
three different wave normal angles 9 . 
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6 (degrees) 

FIGURE 2.6 NORMALIZED PEAK MAGNITUDES OF THE <qvy«jt 
AND <qS2> TERMS AS FUNCTIONS OF WAVE NOR 
MAL ANGLE 0. Both terms are calculated for three dif 
ferent pitch angles. The resonance cone angle for 
f “ 0.5 fu is *60® as shown. 
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We can also use Fig. 2.4 to show that the upper limit on wave 
magnetic field Intensity Is really satisfied, as It was assumed when 
averaging the eqiiatlons of motion. For the parameters of Fig. 2.4, and 
f ■ 5 kHz, a ■ 45® , and 9 ■ 30®, Bu ■ 1.3 x 10® pT, a value much 
larger than the typical field Intensities In the 0.1 to 100 pT range for 
whistler mode waves. Therefore, the required small wave condition for 
the averaging over one gyroperlod Is easily achieved In most cases. 

We have presented a simple set of equations describing cyclotron 
averaged motion of Landau resonant particles In a whistler mode wave 
propagating at an angle to the static magnetic field. We have argued 
that for the parameters of the earth's magnetosphere and for f < fH , as 
It Is the case for the whistler mode waves, this would be a very 
accurate description of the near resonant particles. The fact that the 
equations are compact and simple makes them suitable for analytical as 
well as test particle computer simulation studies presented In the next 
chapters. 


III. ANALYTICAL STUDY OF LONGITUDINAL RESONANCE INTERACTIONS 


A. INTRODUCTION 


In the preceding chapter we derived a set of equations of motion 
(Eqs.2.62, 2.63) for an electron Interacting with a whistler mode wave 
through a longitudinal resonance process. Before using those equations 
In numerical simulations It Is useful to have a seml-quantltatlve 
analysis of that Interaction process, the purpose of which Is to: 

a) Determine, qualitatively, the effects of different 
parameters on the resonance process, and to 

b) Provide a reference for the testing and explaining of 
numerical results. 

From the equations of motion and the resonance condition It Is 
evident that the most important factors that affect the Interaction 
process are: 

1) The magnitude of the wave parallel electric field En 

2) The magnitudes of Bessel terms In the equations of motion 

3) The wave phase velocity Vpi, 

4) The electron parallel velocity v„ 

The variations of Bessel terms have already been discussed In Section 
II. E. 

In the following text we discuss the remaining parameters 
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starting with calculations of expected magnitudes of E.i In the 
magnetosphere. Next we calculate the wave phase velocity Vpn and 
analyze the resonance condition Vpi,« Vm for a wide range of 
magnetospherlc parameters. We also stress the importance of the phase 
between a wave and the In.eracting electrons and examine Its variations. 
Finally, we discuss the energy exchange between the wave and electrons 
through the longitudinal resonance Interaction In an Inhomogeneous 
medium such as Che magnetosphere. 


B. RELATION OF £•• TO Ba AND MAGNITUDE OF Ei. FOR 
WHISTLER MODE WAVES 


Two equations of motion of an electron (Eqs.2.62, 2.63) are 
given In terms of the wave parallel electric field En (Ej^. However, It 
is useful to relate Eh to the wave perpendicular magnetic field Bx (B^ 
because most often wave amplitudes are given and referred to In terms of 
Bx. We proceed now with a derivation of the quantitative relationship 
between En and Bx • 

Using Che plasma dispersion relation (Eq. 2.22) It follows that 


n^ sin6 cos0 Ejj ”*■ • 0 (3.1) 


or 


U) 


n^ sln6 COS0 E ■ - (1 ^ n*sin^0 ) E, 


U)“ 


(3.1a) 
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Furthermore, from Maxwell's equation 7 x E 


we have 

9t 


k COS0 E„ - k sln0 E, ■ 

X- z y 


(3.2) 


Note that we use only amplitudes of $ and s , E and B , because 

z y z y 

both 5^ and 3^ vary as cos(ujt - k • r) . 

Now, substituting E from (3.2) in (3.1a) we have 

X 


n^slne C0.6 ( 

k COS0 



n^sln^e) 


(3.3) 


OI' 


n^sin© ojBy , ? 7 o 

^ + (n^sln^0 + 1 — n^sin^0 ) E, - 0 

k - a 


(3.4) 


w 


Finally, 


or 


n^sin© 0 ) 
0) 2 

k(-V - 1) 

0) 


(3.5) 


Ell 


c n sin9 

f ^/f^ - 1 
? 


B^ 


(3.6) 


Equation 3.6 relates Em to Bj. for whistler mode waves, and it 
can be further simplified if fp^» when it becomes possible to use 

the QL approximation for the refractive index. The refractive index is 
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then given as 


n 


2 


f 2 

P 

f(fjj cos6 - f) 


(3.7) 


and substituting (3.7) for n in Eq. 3.6 the final result is 


E ■ — sin9 

n (f^/f) COS0 - 1 


(3.8) 


Eq. 3.8 was also derived by Helllwell [l965]. It relates Eh to B* for 
whistler mode signals assuming that QL approximation for a refractive 
Index Is valid. 

Equation 3.6 can be applied to any whistler mode signal, 
although It Is possible to derive similar equations for some special 
cases of propagation. One such special case Is a whistler mode signal 
propagating In the Gendrln mode. This mode of propagation Is 
characterized by the Gendrln angle 6^ which can be found by setting 
^ (n COS0 ) - 0. The resulting wave normal angle 0^ is 



It clearly follows from Eq. 3.9 that the propagation In the 

and that 3^ varies from 0** to 
90** as f/fjj decreases from 0.5 to 0. The Interesting properties of 
propagation at the Gendrln angle are summarized as follows: 


Gendrln mode Is possible only If f < ^/2 


1) Substituting (3.9) in (3.7) the refractive index is 
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f 

n„(6„) - - 2 . 

“G'' G'' f 


(3.10) 


11) The phase velocity In the direction of Is 



V 

P 

. cfa. 

(3.11) 

^ G 

COS 0g 

2 f 

p 

111) The 

group refractive 

Index and velocity are 


■ "g<V ■ 

(3.12) 


f 

- V - C -p 
D 

(3.13) 


Iv) The group ray refractive Index and velocity are 


“8r ^®G^ " ‘"°®®G 

G G 


f 

2 

^ f 

H 


(3.14) 


'8r, 




£ £h 

2 f 


(3.15) 


Figure 3.1 Illustrates the shape of the refractive Index curve 

for f/f < 0.5, and also shows the Gendrln angle 0 . The second angle 

Indicated In Fig. 3. 1,6^ , Is the resonance cone angle where the 

refractive Index becomes Infinite. 

Thus, waves propagating at the Gendrln angle have their wave 

packets traveling In the direction of with the velocity Vg^ , which 

G 

Is Identical to the phase velocity In that direction Vp„ , and both 

G 
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velocities are Independent of the wave frequency. This property makes 
Gendrln mode waves rather interesting for longitudinal resonance 
interactions since electrons in resonance with those waves, l.e. 

Vii * Vpii * v^^ , do not drift through the wave packet during the 
interaction aa they do in the most general case when the wave phase and 
ray group velocities along the magnetic field line are different. 



FIGURES.!. REFRACTIVE INDEX SURFACE FOR f<fH/2. 0 r indicates the 
resonance cone where 6 g is the Gendrln angle, for 

which the ray is aligned with the static magnetic field. 


Returning to the derivation of the parallel electric field for 

the Gendrln mode waves we can substitute n(6_)» cos6_ and 
I ■ ' ■ G G 

sin6g " - cos^6^ for n, cos6 and sin6 in £q. 3.8 assuming that 

f_^/f*f„ « 1 is valid. The final result is then 
p n 



(3.16) 


Note that Eq. 3.8 represents the most general expression for En 
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(allowing for the QL approximation) and can also be used to compute 
Eiiq , whereas Eq. 3.16 is valid only for the Gendrln mode. At this 
point we can use Eqs. 3.8 and 3.16 to plot the magnitude of Che parallel 
electric field E«i as a function of frequency. Three curves shown In 
Figure 3.2 are calculated for different values of the wave normal angle 
(30”, 50” and 70”), while the wave perpendicular magnetic field Is 
taken to be 10 pT. This figure clearly shows the resonance cone effect; 
for a fixed wave frequency f Che parall il electric field E„ Increases as 
Che wave normal angle Increases and Ed approaches Infinity as 6 -»■ 6^ . 
The resonance cone angle 6^ can be found from Eq. 3.7 which yields (for 

f 

the QL approximation) cos6„ " and 6 as a function of frequency 

R fn R 

Is Illustrated by the dashed line In Fig. 3.2. At this point we recall 

Chat an upper limit on the magnitude of E„ was already set during the 

derivation of equations of motion when they were time-averaged. 

Although this limit Is not exceeded In most practical cases It Is 

possible that those equations become Invalid In a situation when 

0 - 6- <0.5”. In such a case It would be necessary to use the complete 
R 

eqiiations of motion (Eqs. 2.41, 2.42 and 2.43). 

Figure 3.3 shows the wave parallel electric field Sn as a 
function of frequency and parametric In (10, 20 and 30 pT), while the 
wave normal angle 6 for all curves is 30®. Figure 3.4 shows the wave 
parallel electric field Ehq for the Gendrln mode propagation as a 
function of frequency and parametric in Bj,. The dashed curves show 0 q 
and 9j^ as functions of frequency. By setting Ehq(9q) ■ 0 it can be 
shown that En^ reaches a maximum at the frequency f ■ 0.354 f^ at which 
6g ■ 45 ”. This result is Interesting in the light of data, on chorus 












The two dashed curves indicate the resonance cone and Gendrin angles ^ and 6 g 
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activity obtained by 3urtls [1974]. It was found that In the equatorial 
region there are often cjserved two narrow bands of chorus. The upper 
band Is commonly centered Just above half the electron gyrofrequency, 

0.5 f{{, while the lower band Is centered near 0.35 f[}. Therefore, It 
may be speculated that the chorus lower band Is made up of waves 
propagating In the Gendrln mode and that those waves are amplified 
through the strong longitudinal resonance due to their maximum . 

This wave growth could then account for the observed peak of chorus 
activity. 

Finally, Figure 3.5 shows E„ as a function of wave-normal angle 
6 ; different curves In that figure correspond to different wave 
frequencies, while the Ba Is 10 pT. Again we see the resonance cone 
effect where E„ as . 

AJ.l of the above calculations were done at the equator of the 
the magnetic field line given by L ■ 4 and assuming n„" 400 el/cc. 
Similar calculations can be carried out for different L values and 
corresponding values of Figure 3.6 shows the results of such 

calculations for a range of L values; corresponding values of n^^ used 
In those calculations are shown In Figure 3.7, with a plasmapause, 
cni.'rqcterlzed by the sharp decrease of electron density, located at 
L " 4. The wave parallel electric field Eh Is also normalized by Bj. and 
given In yV/m/pT. From this figure It Is evident that Eh for a given L 
value Increases as the frequency of the signal Increases, as already 
found before (see Fig. 3.2). Furthermore Eh Is larger outside than 
Inside the plasmapause, a fact whlr.h Is directly related to lower 
electron density outside the plasmapause. 





[id/ui/AT/] '^0/ 
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Summarizing, a stronger Eu (for a given Bi) can be achieved by 
increasing the wave frequency, or by raising the wave-normal angle, or 
both. 



FIGURE 3.7 EQUATORIAL ELECTRON DENSITY AS A FUNCTION OF L VALUE. 
The plasmapause is located at L ■ 4. 
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An additional Increase In Em Is also possible for waves 
propagating outside the plasmapause. However, waves with high 
wave-normal angles are usually associated with a non-ducted mode of wave 
propagation which In general Is not field aligned, whereas In the ducted 
mode the wave normals are very nearly aligned with the magnetic field 
[smith et al. I960]. In the latter case guiding Is based on the 
presence of linear field-aligned enhancement (or depression) of 
Ionization referred to as a duct. Therefore, the effects of the 
longitudinal resonance Involving ducted waves are limited by the low 
wave-normal angles of propagation at which magnitudes of the parallel 
electric field are low (see Fig. 3.2). There are other possibilities 
for wave guiding along the field line not limited to low wave-normal 
angle waves, such as when the plasmapause acts as a one-sided duct [inan 
and Bell, 1978]. Still another possibility Is to have a non-ducted wave 
which propagates In a field-aligned mode over a portion of the 
magnetospherlc path. Although those waves usually remain field aligned 
only for a short period of time, their large E„ may be sufficient to 
cause a strong longitudinal resonance Interaction. 

The Importance of field aligned propagation arises from the fact 
that electrons In the magnetosphere follow the earth's magnetic field as 
explained in Section II. A. Thus, if the ray path Is not field aligned, 
or Is only partially aligned, the Interaction may be relatively weak. 
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C. RESONANCE CONDITION v« - Vp„ 


Beside the equations of motion another Important factor to be 
considered Is the resonance condition Vn - Vpn (Eq.2.8). As discussed 
above, this condition requires that the wave phase velocity In the 
direction of B^ match the particle velocity In that direction. However, 
for an Inhomogeneous medium such Is the magnetosphere, both the phase 
velocity Vpii and the electron parallel velocity Vn are variable and 
their variations depend on the magnetospherlc model. Hence, In a case 
when the resonance condition Is satisfied for a given wave and electron 
at some location In the magnetosphere. It will not In general hold at 
some other location. For that reason It Is necessary to study how Vpn 
depends on different models used to represent electron density along the 
field line. It Is also essential to examine variations of both phase 
and parallel velocities with latitude and to study variations of Vn for 
different pitch angles. 

First, let us consider the phase velocity In the direction of B^ 
which is given as 


^P“ “ n • COS0 


(3.17) 


where n Is the refractive Index given by Eq. 3.7. Using Eq. 3.17 It Is 
a simple task to calculate the phase velocity of a whistler mode wave 
for a wide range of parameters. Figure 3.8 shows the equatorial phase 
velocity as a function of L value; values of n^^ used here are again 
those of Fig. 3.6. Figures 3.9a,b show the phase velocity as a function 
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(a) Vpii (km/sl 



(b) Vpii (km/;) 



FIGURE 3.9 PARALLEL PHASE VELOCITY AS A FUNCTION OF LATITUDE FOR 
DIFFERENT MODELS OF THE DISTRIBUTION OF ELECTRON DEN- 
SITY ALONG THE FIELD LINE. In (a) electron density along the field 
line Is represented by the diffusive equilibrium model DE-1, whereas 
In (b) the electron density Is calculated the colllslonless model R-4. 
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of latitude; Fig. 3,9a shows a typical shape of Vpu inside the 
plasmapause, while Fig. 3.9b shows Vpn outside the plasmapause. The 
difference between Figs. 3.9a and 3.9b reflects not only the assumed 
equatorial electron densities n^^ » but also the electron density 
distribution along the field line. Figure 3.9a is calculated using a 
diffusive equilibrium model [Park, 1972 ]» which is usually used inside 
the plasmapause. On the other hand, the electron density model of 
Fig. 3,9b is a collisionless model [Park, 1972] with the electron 
density along the field line approximated by 




1 

cos^X 


4 

) 


(3.18) 


where X is the latitude. 

Evidently, from Fig, 3,9, the phase velocity of whistler mode 
waves outside the plasmapause exceeds that found inside. Therefore, the 
parallel velocity of an electron, which has to match the phase velocity 
of the wave, is also larger outside the plasmapause. Since the 
electrons are moving faster when interactions take a place outside the 
plasmapause the corresponding interaction times are shorter compared to 
interaction times inside the plasmapause. Thus, the effects of a 
stronger wave parallel electric field Eu , related to propagation 
outside the plasmapause, tends to be offset by a reduced interaction 
time. 

The parallel velocity as well as the wave phase velocity varies 
with latitude, as already shown in Section II. A, but the two variations 
are generally different. By combining the first adiabatic invariant and 


a«cv*a 
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the law of energy conservation we find that the parallel velocity is 
given by 


Vii 


Vii 


®q 


1 + tan*a 

eq 


JT. 


3cos^X ^ 2 ^ 

- tan a 


cos 




eq 


(3.19) 


where Vugq is the electron equatorial parallel velocity, is the 
equatorial pitch angle and X is latitude. 



xn 

FIGURE 3.10 NORMALIZED ELECTRON PARALLEL VELOCITY AS A FUNCTION OF 
LATITUDE. Different curves correspond to different 
equatorial pitch angles. Note that the mirror point latitude, where 
Vii " 0, decreases as the equatorial pitch angle Increases. 

Figure 3.10 shows the normalized parallel velocity as a function 
of latitude for different values of the equatorial pitch angle. This 
figure also shows mirror point latitudes where Vh « 0. From Figs. 3.9 
and 3.10 It Is evident that the resonance condition for a given wave and 
electron may, or may not, be satisfied depending on the ratio of the 
equatorial phase and parallel velocities. Typical examples shown In 
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Fig. 3.11 are for three different ratios of the equatorial velocities. 
Note that the parallel velocities shown in Fig, 3,11 represent the 
unperturbed motion of electrons, l.e. Fig, 3,11 shows only adiabatic 
variations of Vh, Although the adiabatic motion of electrons is altered 
by the wave-particle Interaction, the electrons are identified in terms 
of their initial unperturbed equatorial parameters which simplifies the 
problem of comparing properties of different electrons. 

Those different variations of Vpii and Vu with latitude and their 
effects on the interaction process, along with effects of other factors 
are further discussed in the chapters on numerical results. 


D, PHASE BETWEEN WAVE AND ELECTRON IN LONGITUDINAL RESONANCE 

In Chapter II it was shown that the electrons trapped in the 
wave potential well execute an oscillatory motion around the bottom of 
the potential well. In general the analytical solution of the equation 
of motion for that case is very complex, but it is possible to derive an 
approximate solution if the maximum amplitude of the oscillation remains 
relatively small. From Eq, 2,12 the parallel electric field Eh , as 
seen by electrons in the wave frame, is given by 

Ell ■ Eiiq sin(kii*z) (3,20) 


Therefore, the force exerted on an electron is 
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FIGURE 3.11 RELATION BETWEEN v„ AND ALONG THE FIELD LINE. 

Depending on the ratio of /v ^ there may be 

one (a), none (b), or two (c) latitudes at whicK the longitudinal 

resonance condition v„ ■ v_ is satisfied. 

Pn 
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d^z 

dt^ 


q E„osin(k„.z) 


(3.21) 


which, for a small amplitude oscillation where sin(k„*z) - k„.z, can be 
written as 


d^z 

dt^ 


9. 

m 


EiiQkii »z 


(3.22) 


The solution of Eq. 3.22 is 


z ■ z_ sin(o) »t) (3.23) 

0 c 

where z is the position of the bottom of the potential well as shown in 
B 

Fig. 2.2, z is the position of the ele ctron and u)^ is the period of 

/ e E(ig kii 

oscillation given as - J ^ . It should be noted that although 

this oscillation period is computed for a homogeneous medium, this 
result can also be used in the case of a slowly varying medium such as 
the magnetosphere. Now, dividing Eq. 3.23 by the wavelength, we obtain 
the relative phase between the reference point at the potential well 
bottom and the electron. This relative phase is 

sin(oJ *t) (3.24a) 

Ik D C 

The relative phase between the wave and the trapped electron is 
also oscillatory in its nature and the phase variation is bounded such 
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chat i^_< 360“. It should also be noted that the smallest amplitude of 
o 

the phase oscillation- corresponds to the case of strongest trapping. On 
the other hand Che relative phase variation for unCrapped electrons is 
represented by constantly increasing (Vp„ > v.i) or constantly decreasing 
(vpii < vii) phase as those electrons drift backward or forward through 
the wave, respectively. 

All of Che above computations, as already pointed out, are 
carried out in the wave frame which moves in the z direction at Che 
phase velocity Vp„. In order to determine the total phase variation let 
us again assume propagation as exp l(u)*t - k*r). The Instantaneous 
frequency can be found by taking the time derivative (u)*t - k’r) 
which yields 

u>i - u) - k ^ (3.25) 

where Is actually the Doppler shifted frequency of the wave as seen 
by an electron placed at a location defined by radius vector r. It Is 
possible to rewrite Eq. 3.18 In the same form as that of Eq. 2.6 by 
using ^ ■ v„ and substituting for Equation 3.24 can now be 

used to examine a behavior of the total phase between a wave and an 
electron. First, rewriting (3.25) we have 

- u) - k„»v„ (3.26) 

If “ 0 Eq. 3.26 reduces to Eq. 2.7, or 
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‘'pH 


V„ 


(3.27) 


which Is Che original longitudinal resonance condition. Therefore, if 
Vii " Vp„ tlxe relative phase remains constant (Eq. 3.24a). 

However, if an electron has a parallel velocity which does not 
match the wave phase velocity exactly the Instantaneous (Doppler 
shifted) frequency has a non-zero value. In that case both the sign 
and the magnitude of u)^ depend on the difference between the parallel 
velocity and the phase velocity; when vn < Vp„, lo^ is positive and its 
magnitude Increases as v„ decreases assuming Chat Vp„ is constant; in a 
case when vh > Vp„ the instantaneous frequency is negative and its 
magnitude Increases as v„ increases, again assuming a constant Vp„. 

When is knoxm the total phase shift can be expressed as 


(p m 


t 


(3.28) 


or as 


f w 

<|) - ds (3.29) 

a P" 

ds 

where we have used the Identity dt ■ “ — . 

VpH 

Finally, Table 3,1 summarizes qualitatively the behavior of the 
total phase shift as a function of Vpii - vuc 

The phase between the wave and Che electron is a very Important 
factor in the trapping process. It is eventually the phase that 
determines if a given wave will trap any electrons, although all other 
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resonance conditions may already be met« l.e the parallel velocity Is 
clcse to the phase velocity and the parallel electric field Is strong 
enough to pull the electron Into the potential wel^. There Is no 
trapping If the phasing Is wrong, l.e. If electrons are accelerated when 
trapping would require deceleration or vice versa. The numerical 
results will show that a small difference In phase, less than 10**, can 
make a large difference In the behavior of electrons for which the 
resonance condition Vp„ ■ v„ Is satisfied. Furthermore, the phase 
directly translates Into the position of an electron within a wave 
packet (Eq. 3.24) and If there Is any space bunching of electrons there 
must exist a corresponding phase bunching. 


Velocity Conditions 

1 

_< 

V 

o 

0 

V 

1 

Magnitude of Total 

Positive and 

Negative and 

Phase Shift 

Increases with time 

decreases with time 

Rate of Phase 

Increases as Vp„ - v„ 

Increases as Vp„ - v„ 

Change with Time 

Increases 

decreases 


Table 3.1 PHASE SHIFT PROPERTIES OF LONGIDUTINALLY PESONANT 
ELECTRON AS A FUNCTION OF PARALLEL VELOCITY CONDITIONS. 



63 


E. ENERGY EXCHANGE 

In Chapter II we have discussed the energy exchange between the 
wave and trapped electrons in a homogeneous medium. For the case of an 
Inhomogeneous medium the energy exchanged during a longitudinal 
interaction can be computed in a similar fashion. However, we shall see 
later when presenting numerical results that the longitudinal resonance 
In the magnetosphere may, or may not, involve trapping of electrons. It 
will also be shown that electrons in both cases, whether they are 
trapped or not, exchange their energy with a wave. The energy exchange 
process is quite different In those two cases, but it is still possible 
to ise an equation similar to Eq.2.18 by using correct velocity limits 
for integration and an adequate value to represent the energy exchanged 
through the Interaction with a single electron. It is then also 
essential to compare contributions from both groups of electrons 
(trapped and untrapr.<.d), and to determine whether there are situations 
where the contribution from either group Is negligible. 

Here we recall that in the case of a homogeneous medium the 
energy Is exchanged only during the trapping process, i.e. only during 
the period when the electrons are accelerated/decelerated by the wave in 
order to match tha phase and parallel velocities, and there is no net 
energy exchange after that process is finished, or alternatively, an 
electron has to be trapped in order to exchange Its energy with a wave. 
There is still an instantaneous energy exchange after the trapping is 
completed because electrons oscillate at the bottom of the potential 
well, but when this instantaneous energy is averaged over one trapping 
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period there is no net effect. This Is so because the electron's 
oscillatory motion is perfectly symmetric around the bottom of the 
potential well, shown by Eq. 3.20, whereas In the magnetosphere or any 
other Inhomogeneous medium, the energy can also be exchanged after the 
electrons are trapped. This can be explained as follows; after an 
electron Is trapped Its parallel velocity is very close or equal to the 
wave phase velocity and It follows the phase velocity variations as long 
as that electron remains trapped. Thus, the perturbed parallel velocity 
Is different from the parallel velocity that a particular electron would 
have In the absence of the wave. This difference, Avh , Is directly 
proportional to the phase velocity changes [Brice, 19f ?] and it Is given 
as 

3v 3v 

^ "3^ ds + ( ^ df (3,30) 

where. In general, phase velocity depends on both frequency and 
position. For the positive sign of Avh the electron gains energy, while 
for the negative sign the wave gains energy. We shall discuss further 
various aspects of Eq. 3.30 later In the text. 

In the next chapters we present results of a test particle 
simulation of the wave-particle Interaction and Illustrate various 
aspects of the Interaction as they were discussed in the above analysis. 


I 

{ . 
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IV. DESCRIPTION OF THE NUMERICAL SIMULATION 


A. INTRODUCTION 


In this chapter we detail procedures used In numerical 
simulations of the time-averaged equations of motion. The method used 
In this report Is a test particle simulation. This approach uses a 
single particle to find wave Induced perturbations of the particle 
trajectory, and It Is feasible to test quantitatively the effects of 
various factors already considered In a qualitative analysis presented 
In Chapter III. The test particle approach can be further expanded to 
determine the perturbations of a full particle distribution by computing 
the effects of the wave on an adequate number of particles that are 
appropriately distributed In the phase-velocity space. However, there 
are restrictions Imposed on the full distribution simulations because 
there Is no feedback that should account for variations of the wave 
amplitude as particles and the wave exchange their energies. This 
feedback problem Is treated In more detail In a discussion of the 
numerical results. 

The actual listing of the particle code used In all simulations 
presented here Is given In Appendix B. Next we outline the basic 
operation of the program. 
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B. COMPUTATION OF PROPAGATION AND ADIABATIC MOTION PARAMETERS 


The representation of the static magnetic field along the field 

line Is based on a centered magnetic dipole model described by 

Eq. 2.1. Values of Bq obtained from that equation are then used to 

compute local values of the gyrofrequency fj{» as well as to compute a 

1 *^®o 

normalized gradient of the magnetic field ^ ^ . At the same time a 
cold plasma density variation along the field line can be calculated 
using two different models. One model assumes diffusive equilibrium 
[Angeraml and Thomas, 1964] with the electron density along the field 
line given as 





(4.1) 


where the 6 ^ are the relative concentrations of the Ionic species, n Is 
the number of species, G ■ rj^Cl - (r^/r)] , r^ Is the geocentric 
distance (In kilometers) to the base of the DE model, 

, i_i 

■ 1.506T(rjj/7370) (1/4 ), and T Is temperature at the base of the 

DE model (r « 1000 km). A second model Is a colllslonless model for which 

the density Is given by Eq. 3.18. The input parameters needed to uniquely 

define the field line and propagation properties are L value, the 

equatorial cold plasma density wave frequency f, and the 

wave-normal angle 6. Given those parameters the program divides the 

entire field line In spatial segments 10 kilometers long and than 

1 

computes, and stores, values of Vpi, (z), ki, (z), and - — ^ for each 
segment; z Is a distance between the equator and a particular 10 km 
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segment measured along the field line. The stored values of -g— , 

as seen from Eqs. 2.62. and 2.63, are used to compute adiabatic terms in 
the equations of motion. All of the above computations can be done 
either for a general whistler mode wave or for the Gendrln mode wave. 

In the latter case the program also computes, and stores, values of 
6q(z) and Emq(z) . In addition the program also computes, and stores, 
values of wave phase change given as jkiidz. In contrast to other 
parameters the values of Jkitdz are not symmetric about the equator and 
depend on the latitude where the particles are started. This starting 
latitude, l.e location where particles start their motion along the 
field line. Is also one of the Input parameters. 


C. NUMERICAL INTEGRATION OF THE EQUATIONS OF MOTION 


Before we start with simulations each particle must be uniquely 
defined by an appropriate set of parameters. Those parameters then 
describe the particle's position In phase-velocity space. For particles 
in the magnetosphere the velocity coordinate Is uniquely given by their 
equatorial parallel velocity and equatorial pitch angle, Ooeq* 

As particles move along the field line their corresponding equatorial 
parallel velocities can be computed with the help of Eq. 3.19. At the 
same time the local pitch angle Is related to the equatorial pitch angle 
through 



slna 

oeq 


slna 


(4.2) 
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where Bg(z) is the local value of the static magnetic field, and is 

the equatorial magnetic field. 

In this report a given particle is always identified in terms of 
the equatorial parameters which then simplifies the task of comparing 
properties of different particles. The conversion from local to 
equatorial values is made on the assumption of unperturbed particle 
motion. 

In addition to the velocity vn and pitch angle a there 

oeq oeq 

is a third parameter, the initial phase , which determines the 
position of a particle with respect to the wave packet at the beginning 
of the interaction (this is a local, as opposed to an equatorial, 
quantity) . In order to examine the dependence of the interaction 
results on the initial particle phase a simulation is actually done 
using twelve particles uniformly distributed in phase space; the 
parallel velocity and pitch angle are, however, identical for all twelve 
particles. This assembly of twelve particles uniformly distributed in 
phase is called a test sheet and is Illustrated in Figure 4.1. It 
should be recalled that, as already emphasized in Section III.D, uhe 
phase between a particle and a wave is directly related to the 
particle's position in the z~axls direction. This is Important because 
if particles are distributed in phase, l.e. space, the starting time t 
of the Integration must be increased by At ■ — from particle to 

''p II 

particle in order to maintain a correct phase separation between the 
particles in the sheet. This is especially important in particle phase 
(space) bunching calculations where particle positions determine the 
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extent of bunching. 

After particles are injected at a given latitude their motion Is 
altered due to the wave force which Is computed by numerical Integration 
of the equations of motion. A proper value of the startxng latitude, 
for interactions with a monochromatic CU signal as Illustrated In 
Fig. 4.2, was found experimentally by gradually Increasing the distance 
between the first resonance location and the location of particle 
Injection, and finding a latitude where further Increase of this 
distance caused no significant changes of the final results. The actual 
Integration of the equation of motions Is done using a simple 
predictor-corrector method using temporal steps with At • 0.001 msec. 
This time step size was also found experimentally, and for smaller size 
step there were only insignificant fluctuations of the final results In 
all of the examples presented later In the text. The Integration method 
Itself consists In predicting a position of a given particle after 
elapse of one time Increment using current values of force, l.e. using 
those forces acting on the particle at the beginning of the time 
Increment. However, after the particle reaches a new position forces 
acting on It are also different, and It Is necessary to recompute 
(correct) the particle's position by using the average force. This 
average force .U >und as a mean value of two forces, one at the 
beginning and one at the end of the time Interval At. This newly 
computed position of the particle Is then taken as a new starting point, 
and the whole process Is repeated. 

For a case of a monochromatic CW wave particles travel along the 
field line and reach the first resonance point (Fig. 4.2) where the wave 
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induced perturbations of particles trajectories become stronger and 
stronger. At this point further behavior of the particles is very 
dependent on the initial phase . Although all particles have their 
motion altered by the wave forces only a certain class of particles 
becomes trapped, l.e. only those with an appropriate phase, while other 
particles remain untrapped. However, in both cases the Integration is 
continued until all particles reach their second resonant point on the 
other side of the equator. After that moment the wave Induced 
perturbations become smaller and smaller as the difference between 
particles parallel velocities and the wave phase velocity Increases. 

The end point of the integration is then defined as the location where 
the absolute difference between the two velocities exceeds 10%. This 
value was determined experimentally, and the particular latitude where 
the above condition occurs is called the detrapping latitude. 

As the particle moves along the field line from the starting 
point toward a detrap point it has its adiabatic pitch angle variation 
modified by the wave. Finally, after the particle reaches its detrap 
point it will have certain Up and viip which are then transformed into 
the corresponding equatorial values o^Peq ''^"eq using (4.1) and 
(3.19). The difference otogq - Qtpeq gives the total pitch angle 
change, or scattering, while the difference Avh - vi^eq - vupgq gives 
the total energy exchange through 1/2 m Av.,^. The final scattering and 
the amount of transferred energy are given both for each individual 
particle and for a complete test particle sheet (mean value for 12 
particles) . 

In the next chapters we study the scattering of particles and 
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Che energy exchange process for different wave functions and a wide 
range of particle initial parameters. 


V. NUMERICAL ANALYSIS OF THE INTERACTION 


A. INTRODUCTION 


In the previous chapters we have derived a set of equations of 
motion for longitudinally resonant electrons, and we have studied 
analytically various aspects of the resonance process. Those analytical 
studies are now complemented by the results of the numerical simulation 
analysis. Numerical results should further Illuminate the physics of 
the interaction process, and enable us to compare the effects of various 
parameters on a quantitative basis, l.e. In terms of scattering and 
energy exchange efficiencies. The behavior of Individual electrons and 
sheets Is studied for a wide range of the parameters such as En , n^^, L, 

‘J’ » and for different wave functions, l.e. for different wave 
eq o 

amplitude variations along the field line. In our calculations we have 
used three different types of wave functions as they are described 
below: 

a) Monochromatic CW wave with a constant wave amplitude along 
the field line. 

b) One-sided wave function characterized by a very weak wave on 
one side of the equator and a strong wave on the other side. 
The transition region between the above regions is taken to 
be 1000 km long and starting at the equator. Such a wave 

t 

function can be created through a gyroresonance process. 
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c) Spatial amplitude pulse formed by a non-ducted wave 
when Its ray path Is partially field aligned. 

In the following discussion we present results of the numerical 
simulations . 


B. SCATTERING OF A SINGLE SHEET INTERACTING WITH CW SIGNAL 


For a case of monochromatic CW signal the interaction geometry 

o 

is already shown in Fig. 4.2, with electrons being injected at -15 
latitude. All electrons are identified in terms of their equatorial 
parameters, Vh and a , with the initial phase (p^ being a third 
parameter. First, we consider scattering of a single sheet (12 
electrons uniformly distributed in phase at the injection point) as a 
function of the initial equatorial pai^allel velocity Other 

parameters for this example are listed in Table 5.1 below. 


Field Line 

L - 4 

Equatorial Electron Density 

n ■ 400 el/cc 
eq 

Equatorial Gyrofraquency 

f jj » 13.65 kHz 

Equatorial Plasmafrequency 

f = 180 kHz 
P 

Wave Amplitude 

Bj. » 10 pT 

Wave Frequency 

f - 3 kHz 

Wave Normal Angle 

6 » 30° 

Equatorial Parallel Phase Velocity 

V - 9.924 10® m/s 
P"eq 


Table 5.1 PARAMETER VALUES FOR THE EXAMPLE CASE 





At this point we should note that we have used two 
approximations in numerical computations. First, it is assumed that the 
wave-normal angle is fixed, and second, the wave amplitude is also 
treated as though it has a constant value. However, it is well known 
that in the magnetosphere both wave-normal angle and wave amplitude 
change with location. The wave-normal angle changes as dictated by the 
guiding mechanisms [Helliwell, 1965] which is true for ducted waves, 
whereas wave-normals of nonducted waves can be found using ray-tracing 
analysis [Klmura, 1966, Burtls, 1974]. The wave amplitude variation 
arises from the inhomogeneity of the magnetosphere, and it is feasible 
to use a slowly-varying medium analysis to calculata those variations 
[Budden, 1961 ]. From ray-tracing and amplitude calculations it is 
obvious that both the wave-normal angle and the wave amplitude may 
change significantly along the field line, and affect the longitudinal 
resonance interaction. Nevertheless, if the interaction region is 
relatively small, the changes of wave properties are also small, and it 
is permissible to assume as a first order approximation that the 
wave-normal angle and wave amplitude are constant quantities. If there 
is a need for even more accurate analysis it is feasible to use 
ray-tracing along with WKB solution to derive exact solutions for both 6 
and Bi , and then Incorporate those results in the longitudinal 
resonance calculations. 

The mean scattering, <AOgq> (< > denotes averaging over the 
initial phases), of a single sheet of electrons as a function of sheet 
equatorial parallel velocity is Illustrated in Figure 5.1. The wave 
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Intensity ■ 10 pT corresponds to Eh" 15 yV/m. A solid curve shown In 

that figure Indicates the mean final scattering of a sheet at the end 

point of the Integration of equations of motion (as defined In Fig. 

4.2), while the dashed curve represents the mean scattering of a sheet 

computed at the equator. Comparing the equatorial, l.e. cumulative 

scattering when electrons reach the equator, and the final scattering It 

Is obvious that the final scattering Is, on average, one order of 

magnitude smaller than the equatorial scattering. It Is also clear from 

Fig. 5.1 that the equatorial scattering Is negative, l.e. the mean 

equatorial pitch angle of twelve electrons forming a sheet Is lowered. 

To explain those results shown In Fig. 5.1 It Is useful to study 

trajectories of Individual electrons. For example Figure 5.2 

illustrates typical electron trajectories and energy variatiors 

calculated for Interactions with a monochromatic C1<J signal. Four 

electrons shown In Fig. 5.1 belong to a test sheet specified by 

Vii ■ V , and a ■ 10®. A main difference between those electrons 
eqo Pi. eq eq 

are their Initial phases as Indicated In Fig. 5.1 and defined In 
Fig. 4.1. The left column of Fig. 5.2 shows energies of the four 
electrons as a function of Interaction time, while the right column of 
the same figure Illustrates variations of both parallel and phase 
velocities as a function of latitude. Note that the time scale and the 
latitude scale cover the same portion of the field line. Next consider 
Fig. 5.2a where, as the electron approaches the equator, the parallel 
velocity becomes better matched to the wave phase velocity, and the wave 
effects become more cumulative. Those wave effects cruse the 
oscillations of Vh and E, and as the electron comes closer to the 


FIGURE 5.2 SINGLE ELECTRON TRAJECTORIES FOR Bj. = 10 pT. The electron 
energy and parallel velocity are shown as a function of 
latitude as it interacts with CW wave. The initial parallel velocity 
'^•'eoq “ Vpneq * and a " 10® for all electrons. The initial phase <()o is 
30® in (a), 90® in (b), 120® in (c), and 270® in (d) . 
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equator the amplitudes of the oscillations increase. At the point 
t ■ 0.52 sec (X» -3.5* ) the parallel velocity of the electron equals 
the phase velocity, and that point is called the first resonance point. 
Electrons shown in Figs. 5.2b, 5.2c, and 5. 2d exhibit similar behavior 
before they reach the first resonance point. However, after electrons 
travel beyond the first resonance only the top three electrons shown in 
Fig. 5.2 are accelerated by the wave in such a manner that their 
parallel velocities become larger than the phase velocity. It is also 
clear from Figs. 5.2a, 5.2b, and 5.2c that this Increase of the parallel 
velocity is accompanied by an Increase of the total energy of the 
electrons. After those electrons have traveled beyond the first 
resonance their motion, as they travel across the equator, is still 
affected by the wave, but the parallel velocity remains larger than the 
phase velocity. However, on the other side of the equator the phase 
velocity again starts to Increase and the electrons approach their 
second resonance point. At this second resonance point the electrons 
are decelerated by the wave and consequently their energy is also 
decreased. Thus the electrons shown in Figs. 5.2a, 5.2b, and 5.2c are 
being accelerated at the first resonance point and then decelerated at 
the second resonance point. The amount of acceleration and deceleration 
in general depends on the actual phase between a given electron and the 
wave, and as a final result electron energy can be unchanged (Fig. 

5.2a), increased (Fig. 5.2b) or de^' < biased (Fig. 5.2c), Compared to 
those top three cases (Figs. 5.2a, 5.2b, 5.2c) a fourth electron 
trajectory Illustrated in Fig. 5. 2d is quite different. This electron 
became trapped after the first resonance Interaction and Its parallel 
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FIGURE 5.3 SINGLE ELECTRON TRAJECTORIES FOR = 10 pT. The electron 
parallel velocity Vn and phase (J> as a function of time. 

Time t » 0 indicates occurrence of the first resonance. Other parameters 
are the same as those in Fig. 5.2. 
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velocity, as well as the total energy, shows oscillatory behavior which 
is characteristic of the trapped electrons. 

Figure 5,3 is a time expanded view of the electron's behavior 
during a 400 msec window centered around the first resonance point at 
t ■ 0 msec. This figure shows both parallel velocity and electron phase 
behavior. From the phase diagrams it follows that the phase is 
increasing before the first resonance, with the rate of Increase 
decreasing as electrons approach the first resonance point. This type 
of phase variation is consistent with that found analytically in Chapter 
III. At the resonance point the phase does not change, l.e. it becomes 
constant, and the first derivative is equal to zero, as indicated in 
Fig. 5.3* After the first resonance untrapped and trapped electrons 
undergo different phase v«)rlatlons. Untrapped electrons are associated 
with a constantly decreasing phase as a result of v„ > Vp„ , while 
trapped electrons exhibit an oscillatory phase behavior as they 
oscillate at the bottom of the potential well. Note that an electron is 
considered to be trapped if it executes at least one complete phase 
oscillation. Figure 5.3 also clearly illustrates significance of the 
phase between electrons and a wave. By comparing the phase behavior of 
the electrons shown in Figs. 5.3c and 5.3d, we see that the difference 
in their phases at the resonance point (t * 0 msec) is less then 5 
degrees, but the electron of Fig. 5. 2c is not trapped, whereas the 
electron of Fig. 5. 2d is trapped. 

Those four sample trajectories are representative of typical 
perturbations of electron motion induced by the wave forces. Finally, 
to explain the results of Fig. 5.1 where the equatorial scattering is 








83 ORIGINAL PAGE IS 

OF POOR QUALITY 





X[»] 

FIGURE 5.4 NORMALIZED ENERGY OF TEST SHEET AS A FUNCTION OF LATITUDE. 

The normalized energy of a test sheet (12 electrons) in- 
creases about 2.5% around the equator when those electrons interact 

with a CW signal. The sheet initial parallel velocity is Vugg* ■'^nitea 
and a - 10®. h P 4 
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larger than the final scattering, the energies of all 12 electrons are 
added together and plotted as a function of latitude In Figure 5*4. 

From this figure It Immediately follows that there Is a region around 
the equator where the normalized total energy of the electron sheet Is 
Increased. This energy Increase Is on average about 2Z of the Initial 
total energy, £md it Is limited to latitudes between -A** and 4**. The 
Jump In the energy Is caused by the acceleration of untrapped electrons 
such as those shown In Figs. 5.2a, 5.2b, 5.2c, while the energy envelope 
oscillations are caused by trapped electrons such as that of Fig. 5. 2d. 

In the particular example there were 7 untrapped electrons and 5 trapped 
electrons. Beyond X 4** the total energy of the sheet returns almost 
to the Initial level. Here we recall that an Increase of the electron 
energy yields a decrease of the pitch angle, while a decrease of the 
electron energy yields an Increase of the pitch angle. Bearing this 
relation In mind It Is then easy to explain the results of Fig. 5.1 by 
translating energy variations shown in Fig. 5.4 Into pitch angle 
variations. This transformation Ijnmedlately reveals that the equatorial 
scattering Is negative and larger than the final scattering, again as 
Indicated In Fig. 5.1. It also explains why the final scattering can be 
both positive or negative because the final energy can be either larger 
or smaller than the Initial energy. The final scattering appears, due 
to Its randomness, as though It resulted from an Incoherent Interaction. 

On the other hand the equatorial scattering appears to be much less 
random Implying a larger degree of coherence. This Indicates that 
coherence of this particular type of longitudinal Interaction Is 

position dependent, and it Is necessary to examine electron trajectories | 

% 

„j 
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rather than to rely only on scattering results. 

The energy gained by the electrons is extracted from the wave 
which means that the wave amplitude must be reduced around the equator. 
For test particle studies involving only twelve particles this 
attenuation of the wave amplitude is negligible, but it should be 
considered in full distribution computations where significant loss of 
the wave energy will cause a strong wave attenuation and consequently 
weaken the Interaction process. 

From Fig. 5. 3d It follows that the trapping period Is about 82 
msec. Because the medium Inhomogenelty Is very small around the equator 
this trapping period can be also computed using a relation derived for 
the homogeneous medium 



Using (5.1) with k„ « 1,9 10”^ and E„ - 15 y V/m, the trapping period is 
computed to be 81.5 msec, which is In very good agreement with the 
numerical result. It Is also easy to check the oscillation period of v„ 
for untrapped electrons. For example consider the electron shown In 
Fig. 5.2b and Its parallel velocity at t 100 msec. The period of 
parallel velocity oscillation at that point is about 20 msec, which may 
also be found by computing the doppler shifted frequency of the wave 
( 1 ) - kiiVii, Taking co ■ 2ir* 3000 rad/sec, kii ■ 1.9 10 ^ , and Vn ■ 

10100 km/ sec yields 331 rad/sec; the equivalent oscillation period 
is of about 19 msec which is in a good agreement with numerical results. 

As mentioned earlier results shovm in Figs. 5.2, 5.3, and 5.4 
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FIGURE 5.5 SINGLE ELECTRON TRAJECTORIES FOR = 10 pT. Electron 
parallel velocities as functions of latitude for a case 
when the initial parallel velocity is Vng* 1.050 Vp„gq. The pitch 
angle and initial phases of electrons are the same as those in Fig. 

5.2. 
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are calculated for a sheet with initial equatorial parallel velocity 
v„eqg equal to the equatorial phase velocity of a wave. For 

purposes of comparison. Figure 5.5 shows the parallel velocity behavior 
of four electrons, from a sheet with v„„_ » 1.050 v_,,-^, and again as 

a function of latitude. The motion of the electrons is similar to that 
shovm in Fig, 5.2. The possibility of trapping, or not trapping, 
depends on the Initial phase ((i^ of nach individual electron, and the 
final scattering can be both positive and negative. 

The above results suggest that the longitudinal resonance 
interaction with a monochromatic CW signal is confined to a relatively 
small region around the equator. The controlling factor in the 
interaction is the variation of phase <p which determines if electrons 
become trapped or not, and affects the amount of exchanged energy. 


C. SCATTERING OF A SINGLE SHEET INTERACTING WITH CW WAVES 
AMPLIFIED AT THE EQUATOR THROUGH THE CYCLOTRON RESONANCE 

Next we consider the scattering of single electron sheet 
interacting with a monochromatic CW wave whose amplitude is Increased 
through the gyroresonance process. The amplification process of CW 
waves takes place close to the equator [Helliwell, 1967], and In our 
calculations the growth region is taken to be 1000 km long. The wave 
amplitude, before it reaches the equatorial growth region, is 0.1 pT. 

Figures 5.6 and 5.7 illustrate the scattering of a single sheet 
as a function of the initial parallel velocity Vugq^ . In all 
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FIGURE 5.6 MEAN SCATTERING AS A FUNCTION OF PARALLEL VELOCITY. The mean scattering is computed for 
longitudinal resonance Interactions with a CH signal which is amplified after it crosses 
the equator. The result^i indicate that the wave magnetic forces become important at larger pitch ang 
les; then it is necessary to Include Bessel terms in the equations of motion. 
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computations the wave uiplltude is Bj. ■ 10 pT, or E„ ■ 15 uV/rn, while 

the equatorial pitch angle is taken to be 10®, 30°, 50°, and 70°, The 

total sheet scattering is computed twice for each parallel velocity 

Increment; once it is computed using complete averaged equations of 

motion, and once using only the qE term of Eq.2.61 as though the wave is 

electrostatic, i.e, it is assumed that Jq(H) ■ 1 and Jx Cl) ■ 0, As 

discussed earlier, the effects of the Bessel terms, i.e. the effects of 

the wave magnetic field forces, should become significant at larger 

pitch angles, while at lower pitch angles the difference between the two 

computational methods is expected to be small. From Fig. 5.6a, which is 

calculated using a " 10°, it is evident that the two methods produce 
eq 

very similar results, as expected. On the other hand, as the pitch 
angle Increases the difference between the results becomes much larger 
and for ®eq " there is almost no scattering if v?e exclude the Bessel 
terms from the equations of motion (Fig. 5.7b), whereas che scattering 
calculated using the complete eqtiations is about *-6° at Vneq * Vp„eq* 

Those examples confirm the results of Chapter II, where it was found 
that the Bessel terms will be a very Important factor in governing the 
motion of electrons with high pitch angles. This is especially true for 
the Ji(n) term, which represents effects of the wave magnetic force, as 
already Indicated in Figs. 2.5, 2.6 and 2.7. 

As discussed earlier the longitudinal resonance interaction 
depends strongly on the wave amplitude. Tills wave amplitude dependence 
iff depicted in Figure 5.8. Three different curves shown in that figure 
represent scattering of sheets with three different initial parallel 
velocities Vneq^ • A sheet with Vneq^ * Vpiigq has the optimal parallel 

C I 
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velocity as required by the resonance condition. Two other sheets with 
'^"eqo " ^P"eq ^"eqo " Vp„eq are slightly off the 

resonance when they encounter the wave growth region at the equator; the 
first is slower and the second is faster than the phase front of the 
wave, respectively. The effects of different sheet velocities are best 
Illustrated by considering the amount of pitch angle scattering for a 
given wave amplitude. The particle sheet with v„gq - Vpi,gq is scattered 
about - 0.1° when interacting with a relatively weak wave with ■ 5 
pT. On the other hand, the other two sheet® require a wave with ■ 18 

pT to achieve the same amount of scattering. Below B^ ■ 18 pT 
scattering of the sheet with Viigq * 0.995 Vpneq is small and negative, 
whereas scattering of the sheet with Vii^q 1.005 Vpiigq is also small, 
but positive. We recall from Section III.E that the direction of energy 
exchange depends on the relative magnitudes of the parallel and phase 
velocities; if an electron is faster than a wave it is decelerated and 
loses Its kinetic energy; if an electron is slower than a wave it is 
accelerated and gains kinetic energy. An Increase, or decrease, of the 
kinetic energy is accomplished by changing the parallel velocity of the 
electron through the resonance process. If the parallel velocity of an 
electron is Incr. eased, its equatorial pitch angle becomes smaller, or 
equivalently, if the parallel velocity of an electron is decreased, its 
equatorial pitch angle becomes larger. It is this type of process that 
explains the behavior of the two sheets with Vi^q^ 0.995 Vp„eq and 
Viieqg“ 1.005 Vpiieq for Bj_ < 18 pT. It may be wondered why a sheet with 
viieq^ * Vpiieq does not show similar behavior, and what is happening when 
> 18 pX in the other two cases. The answers may be found by 
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examining trajectories of individual test electrons. From those results 
it was found that for weak waves all electrons remain untrapped 
regardless of their initial parallel velocities. As long as the electron 
is not trapped, l.e as long as the electron parallel velocity does not 
follow the phase velocity variation, the longitudinal interaction is 
generally limited to two relatively small regions around the two 
resonance points. In our case the interaction is further limited to 
only one side of the equator where the wave amplitude is sufficiently 
strong. Next, as the wave amplitude Increases beyond the equator the 
interaction becomes stronger, and from the trajectory calculations, it 
is evident that some electrons become trapped. This transition between 
the untrapped and trapped mode of the longitudinal interaction is 
characterized by a significant Increase in the scattering. The 
amplitude threshold at which the trapped mode scattering overtakes the 
untrapped mode scattering depends on the initial parallel velocity Vtigq^ , 
as shown in Fig. 5.8. The threshold amplitude for - Vp'ieq 

low as " 3 pT, with a relatively smooth transition between the two 
interaction regimes. The amplitude threshold in the two other cases is 
about Bj, 18 pT with a much sharper transition between two interaction 
regiaes. 

The individual particle trajectories are illustrated in Figures 

5,9, 5.10, and 5.11. Figure 5.9 shows parallel velocities and phases of 

four electrons with v„gq^ ■ Vpiigq, a ■ 10® , and different initial 

phases d) , as functions of latitude and time, respectively. The wave 
o 

amplitude is B^<* 10 pT. As in the case for a CW signal the parallel 
velocity variation of those electrons is controlled by the phase 
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lue. The threshold amplitude Increases as the absolute difference between the initial parallel 
and phase velocity becomes larger. 
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FIGURE 5.9 SINGLE ELECTRON TRAJECTORIES FOR Bi = 10 pT. Paral- 
lel velocity and phase behavior for electrons with 
v„eq"Vpeq and a » 10“ interacting with variable amplitude CW 
signal. The initial electron phase is 0* in (a), 120® in (b), 
150° in (c), and 300° in (d) . 
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FIGURE 5.10 SINGLE ELECTRON TRAJECTORIES FOR = 10 pT. Shown here 
are the parallel velocity and phase variations around the 
first resonance point at t»0. Other parameters are same as those in 
Fig. 5.9. 
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variation. For example, the electron trajectory of Fig. 5.9a Indicates 
absence of trapping because of an improper phase, whereas the number 
of oscillations for trapped electrons in the other three cases also 
depends on the phase at the moment when the parallel velocity equals the 
wave phase velocity. Figure 5.10 depicts a time expanded view of the 
electron trajectories around the first resonance point. Before 
analysing those trajectories we recall from section II. B that the 
variation of Z,i is described, in the vave frame, as cos kiiz, and that 
the bottom of the potential well is at as shown in Figure 2.2. In 
Figure 5.10 the time t » 0 indicates the first resonance where 
Vii ■ V . The phase at this point is a crucial factor governing the 

Pii 

further motion of a particular electron. For example, the phase of 
electron shown in Fig. 5.10a is such that it is strongly decelerated and 
by the time of phase reversal, l.e. electron acceleration, the parallel 
and wave phase velocity are too different for trapping to be possible. 
Observing the phase of the electron in Fig. 5.10b at t ■ 0 we find this 
phase to be significantly smaller than the phase in Fig. 5.10a. Due to 
this different phase the second electron is less decelerated, eventually 
becomes trapped, and executes one oscillation at the bottom of the 
potential well. For the next two electrons shown in Figs. 5.10c and 
5.10d the phases at t - 0 are even smaller resulting in an increasing 
number of oscillations. We note that the amplitudes of both velocity 
and phase oscillations decrease as the phase at t * 0 decreases. In the 
example shown in Fig. 5.10d the phase at t " 0 is vexy close to the 
optimal 90® which then results in the strongest trapping. As discussed 
earlier the 90® phase indicates that an electron is exactly at the 
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bottom of the poteiitxal well. To illustrate the effects of wave 
amplitude Figure 5.11 shows the same four electrons, but the wave 
amplitude Is Increased to ■ 30 pT. In this case even the first 
electron becomes trapped, and the other three electrons now remain 
trapped for longer periods of time. 

Figure 5.12 shows the scattering of Individual electrons as a 
function of their Initial phases for three different wave amplitudes. 
This figure confirms the importance of phase as a controlling factor In 
the longitudinal resonance Interaction. Figure 5.12 shows that It Is 
possible to achieve a significant Increase of the scattering efficiency 
by changing the Inltal phase (1)^ from 0** to 180°. We summarize the 
results of the above analysis In Figure 5.13 which shows the normalized 
total energy of a single sheet as a function of latitude. The Initial 
equatorial parallel velocity equals the equatorial phase velocity and 
wave amplitude Is Bj, ■ 10 pT. Before electrons reach the equator the 
wave amplitude Is very small and there are no significant changes of the 
sheet energy. After the equator crossing the wave amplitude starts to 
Increase and electrons become trapped. As long as those electrons 
remain trapped their parallel velocities Increase and so does the total 
energy of the electron sheet. As the electrons move away from the 
equator some of them become detrapped, but the energy Increase continues 
up to the point where the last electron becomes detrapped. At that 
point the energy of a sheet has reached Its maximum and remains 
constant. From Figure 5.13 we see that the particular sheet has gained 
about 4.6% over the Initial energy. The energy gain region Is between 
X " 1° and X • 7**. Recall that this energy Increase must be accompanied 



FIGURE 5.12 TOTAL SCATTERING, Aoeq. VERSUS INITIAL PHASE FOR DIFFERENT 
WAVE AMPLITUDES. The Initial pitch anyle is oieq" 10*. 
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FIGURE 5.13 NORMALIZED ENERGY OF TEST SHEET AS A FUNCTION OF LATITUDE. 

The energy of the test sheet Is Increased as It Interacts 
with the variable amplitude CW signal. 
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by wave attenuation which is not considered In the test particle 
studies, l.e there Is no feedback to account for wave amplitude changes. 
The feedback effects can be neglected In a test particle simulation 
where the number of electrons is sicall, but they must be considered In a 
full distribution analysis. 

Nejct we take Into account the scattering efficiency dependence 
on the wave'-normal angle. Figure 5. 14 shows <Aagq> vb. Q for 10 

pT,Qeq “ 10** Bxid v„Qq ■ '^p”eq* wave function corresponds to one 

given In Fig. 5.3. Also shown are the initial energy of the sheet, ^ ^ 
and the final energy < Eg >. We have found earlier that the main effect 
of the wave-normal angle Increase is seen through an Incraase of Eh . 
Thus, as the wave normal Increases the longitudinal Interactions become 
more effective, as Indicated In Fig. 5.14. Furthermore, when the 
wave-normal approaches the resonance cone electrons are scattered by as 
much as -5.5", and the sheet energy Is Increased about five times. For 
such a strong Interaction the wave amplitude would most likely be 
heavily attenuated, although to find the exact solution It Is necessary 
to include a previously discussed feedback term. The Inclusion of the 
feedback term tniuld than probably diminish the scattering effects as the 
wave amplitude becomes smaller with the Increasing scattering. 

In Chapter II we discussed the possibility of space bunching of 
electrons through the longitudinal resonance process. Figure 5.15 shows 
the phases of nine electri^ns from a sheet with “ ''p*'eq » Oeq “ 

and Interacting with a 30 pT wave. Three remaining electrons are 
omitted from this figure because they are very weakly trapped as already 
Illustrated In Fig. 5.12. Initially all electrons are uniformly 
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distributed in phase space and maintain this phase separation as they 
approach the equator. Ac the equator they reach a wave growth region 
and trapping takes place. As electrons become trapped around t 0.21 
sec their maxlmvim phase separation is reduced to about 150*’, and can be 
as small as 50*^ at the moment when all electrons reach the bottom of 
potential well nearly simultaneously at t>0.21, t>0.24, and t»0.27 sec. 
Thus the original spacing between the electrons Is reduced and we have 
a case of space bunching. In this particular example 9 out of 12 
electrons are bunched In about a half of the original separation. Thus, 
the density Increase Is roughly 9/12 ^ 360/130, or about 180% of the 
Initial density for Vneq^ ■■ Vpngq. For other velocities the density 
Increase Is smaller because the resonance condition Is not satisfied 
exactly at the equator. Note that after a few Initial oscillation 
periods electrons go out of phase and start to reach the bottom of the 
potential well at different times. It Is possible to have a new 
synchronization later In time, as occurs at t ■ 0.54 and t > 0.565 sec 
(Fig. 5.15). This problem may be understood as though we have 9 
harmonic oscillators with slightly different periods of oscillation 
caused by different phases at the moment those electrons entered the 
trap. 

Figure 5.16 shows vs. Vn ,and <AE> vs. v„ for 

Interactions taking place Inside and outside the plasmapause. Those 
results clearly show that Interaction? outside the plasmapause result in 
less scattering, but In more energy exchange, than those Interactions 
Inside the plasmapause. This Interesting result may be explained as 
follows; as n^^ drops outside the plasmapause the wave phase velocity 
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FIGURE 5.15 PHASE BUNCHING DUE TO THE LONGITUDINAL RESONANCE INTERACTIONS. 
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Increases and the parallel resonant energy becomes higher. Higher 
energy electrons move faster through the wave and hence have a shorter 
time to be scattered. Note that although the resonant energy Is about 
288 eV for n^^ ■ 400 el/cc It is 11529 eV for n^q "10 el/cc. Because 
of that difference In resonant energies even a relatively small 
scattering outside the plasmapause results In energy changes that are 
larger compared to those found Inside the plasmapause. 

This concludes our discussion of single sheet scattering 
Interacting with a one-sided wave function. In the next section we 
present results involving sheet scattering by a spatial pulse. 


D. SCATTERING OF A SINGLE SHEET INTERACTING WITH A SPATIAL PULSE 

In this section we examine the scattering of a single electron 
sheet as it moves through a spatial amplitude pulse formed by a 
non-ducted wave when Its ray direction stays field aligned for a certain 
portion of the wave path. As depicted In Figure 5.17a the ray direction 
Is field aligned between X • - lO** and A - - 7”, which Is equivalent 
to 1000 km In length. Other Interaction parameters are specified In the 
SL^e figure. The Interaction Is studied for a wide range of Initial 
parallel velocities, Av„, as Illustrated In Figure 5.17b. The minimum 
parallel velocity Is 1.012 Vpii^q, the maximum parallel velocity Is 1.106 
Vpiieq, and the parallel velocity Increment Is 0.001 Vpiieq. The wave 
amplitude is assumed to be zero everywhere except for - 10® < A < -7®. 

The scattering results are shown In Figure 5.17. To explain those 
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results we can use Figure 5.17b as follows; when the initial parallel 
velocity is small, for example v„pgq - 1.012 Vp„gq, the latitude of the 
first resonance point is also small, i.e. it is close to the equator. 
Hence, as those electrons travel up the field line toward the equator 
they encounter the spatial amplitude pulse but parallel and phase 
velocities are rather different resulting in a very weak Interaction. 

As the Initial parallel velocity of a sheet is Increased the first 
resonance point moves away from the equator and closer to the amplitude 
pulse, and the two velocities become better matched. This better 
velocity match results In a stronger Interaction and a negative 
scattering . A negative sign of <Ao^q> means that electrons are 

accelerated. This acceleration Is consistent with the relative ratio of 
two velocities; namely, before electrons reach the first resonance point 
their velocity Is less than the wave phase velocity In which case 
electrons are accelerated In order to match the phase velocity. 

However, further Increase of the parallel velocity beyond 1.082 Vp„gq 
results In a change of sign of the effective scattering. This occurs 
when the first resonance point falls within approximately ±0.5® of the 
pulse front edge at - 10®. The principal difference Is that electrons 
become trapped as they Interact with the pulse, whereas for lower 
parallel velocities there were no trapped electrons. When trapping 
takes place the parallel velocity follows the phase velocity, which 
decreases as electrons approach the equator, and this results In a 
positive sign of scattering <Attgq> in Fig. 5.18. Furthermore, as the 
parallel velocity is increased beyond 1.094 Vp„gq the first resonance 
moves even further down the field line and interactions become small 
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FIGURE 5.17 INTERACTION WITH SPATIAL AMPLITUDE PULSE 
EXTENDING BETWEEN A=-10® AND X=-7®. Shown 
In (a) is the position of spatial pulse on the field line. 
The range of affected initial parallel velocities is shown 
in (b) . 


ORIGINAL PAGE IS 
OF POOR QUALITY 



3 V 
h- M 
Ul Q 
OQ 

(U 

CS 00 
Z C 
•— <0 
o u 


>; > T3 

Ul < 0) 

a 

Ut 01 o. 
(/) X 10 
-J U h 
5 4J 

a. e e 
•H a 

Ui £ 

o 4J e 


_j CO qj 
Q. (U h 

5 

< 4J 0) 
•H 01 

^ U <H 

< O U 

iH "H 

h- 0> O 

SC > O 

a» rH 
</> iH 0» 

X a ii 

K- B H 


10 

<0 

a 


LO 4J 
Z -H 
O 5 


>-> <n (0 

u c -H 

< O 4J 

Ul <J c 
O -H 


s 01 
>-* iH 

Cl 

U 

01 

ac 

JS 

c 

<J 

u • 

o 

sfl r>« 

JS 

Z 1 

4J 

II 

•H 

cac r< 

» 

Ul 

i=i 

CO 

s 

«c < 

o 

c_> 

h 

CO 

u 

o 

o 

Z r— 

01 

< 1 


Ul II 

01 

z ^ 

CO 

01 

rH 

•H 

.e 


9 

U) 

m 

Ul 

•o 

0) 

oc 

a 

= 

a 

CP 

(0 

>— < 

Vi 

u 

4J 



ORIGINAL PAGE 19 
OF POOR QUALITY 




FIGURE 5.19 INTERACTION WITH SPATIAL AMPLITUDE PULSE EXTENDING 
BETWEEN AND X=10°. The format is the same as 
that of Fig. 5.X8. 
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again. The shaded area in Fig. 5.17a Indicates the trapping velocity 
bandwidth A which la also Indicated In Fig. 18. When comparing 
areas of positive and negative scattering In Fig. 5.18 they turn out to 
be approximately the same which means that the energy exchange Is small. 

This example Is a good Illustration of the different features of 
the longitudinal resonance Interaction. We see that the electron 
behavior Is very dissimilar In cases with and without trapping. 

Untrapped electrons change their velocity depending on the relative 
ratio of phase and parallel velocities, while trapped electrons become 
space bunched and their parallel velocity follows the wave phase 
velocity. 

Figure 5.19 Illustrates a similar type of interaction as the one 
discussed above, only the spatial amplitude pulse is on the other side 
of the equator. The corresponding scattering results are shown in 
Figure 5.20. Those results may be explained using the same analysis as 
the one used in the previous example. The trapping occurs when the 
first resonance point is close to the pulse front edge at X • 7^, 
although the trapped electron scattering is now negative as the phase 
velocity Increases. The untrapped particle scattering is positive 
because the phase velocity is smaller than the parallel velocity before 
the resonance point is reached. 


VI. FULL DISTRIBUTION SIMULATIONS 


A. INTRODUCTION 


In Chapter V we have presented results of single sheet 
simulations. The purpose of that analysis was to clarify various 
aspects of the longitudinal resonance process. In this chapter we carry 
those calculations one step further by Increasing the number of test 
electrons In order to simulate a full distribution. Such calculations 
are Interesting for two reasons: 

1) It Is possible to calculate a precipitated flux, and 

2) It Is feasible to estimate wave amplitude changes due 
to the energy exchange. 

In the following examples of full distribution calculations 
electrons are assumed to interact with a one-sided wave function. As It 
was already shown In Chapter V, this type of wave function may produce a 
significant amount of scattering, whereas Interactions with narrowband 
signals (not amplified through gyroresonance) may result In a very small 
final scattering. Therefore, based on those results. It appears that 
the constant amplitude CW signals represent a very weak source of 
precipitation, although those CW waves still may have some amplitude 
variations around the equator as a consequence of the Interaction with 
electrons. 

The energetic electron population Is readily described In terms 


113 


114 


of an equatorial distribution function f^q (v,ieq ,oieq) • From this point 

on we drop the subscript 'eq', and all quantities represent equatorial 

values unless specified otherwise. The distribution function xs given 

in Vii - a space because It Is a convenient representation which directly 

shows the pitch angle scattering, M , and It is easy to determine a 

normalized velocity vn/vpii which Is one of the prime factors affec 'lng 

the Interaction process. The velocity space volume element Is then 

gi'^enas v* dadv„d(j) [inan, 1977]. 

cos^a 

Now we recall results of Figures 5.6 and 5.7 showing the mean 
scattering of a single sheet as a func.:lon of the sheet Initial parallel 
velocity. From those figures It is evident that the trapping velocity 
range considered Is limited to a narrow strip around v„ ■ vpn , while 
the pitch angle range extends from to . The value of may 

be as large as 90**, and specifically In our calculations It may be 
limited to a slightly lower value due to time averaging In the equations 
of motion. The angle 5.5*^ Is the nominal loss cone angle for the 

dipole field lire L " 4, l.e. all electrons with pitch angles lower than 
5.5® have mirror points at Ionospheric heights (h < 200 km) and are 
assumed to be lost through precipitation. As already shown In Figs. 5.6 
and 5.7, the trapping velocity bandwidth Increases with Increasing pitch 
angle due to the effects of the wave magnetic field forces. This 
trapping velocity bandwidth Av„(. Is about 0.4% of Vp„gqfor g - 10®, and 
about 1% of Vpiigqfor g " 70®. Again, It should be noted that this 
velocity bandwidth refers to the trapped electrons only. The untrapped 
electrons have a quite different behavior; if the Initial parallel 
velocity Is smaller than the lower trapping velocity limit the 
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scattering Is negligible because the wave phase velocity and the 
parallel velocity of the electron are never matched along the field 
line. On the other hand. If the Initial parallel velocity of an 
untrapped electron Is larger than the upper trapping velocity limit 
there are always two resonances; at the first resonance scattering Is 
negligible because the wave amplitude Is very small, whereas at the 
second resonance point, where the electron parallel velocity exceeds the 
wave phase velocity, the untrapped electrons are decelerated. All of 
the above mentioned classes of electrons are Illustrated In Figure 6.1. 
The scattering of untrapped electrons is much smaller than It Is for the 
trapped electrons, but the Interaction velocity range for untrapped 
electrons Is larger than the trapping velocity bandwidth. The effects 
of trapped and untrapped electrons on the wave amplitude are exactly 
opposite; the trapped electrons are accelerated and the wave loses 
energy, whereas the untrapped electrons are decelerated and the wave 
gains energy. This dissimilar behavior of trapped and untrapped 
electrons indicates that. In order to' calculate a net transfer. It Is 
necessary to consider a wide range of Initial parallel velocities of 
electrons which then requires a very large number of test electrons. 
While the wave amplitude variation calculations require a large number 
of test electrons the precipitation calculations may be carried out by 
considering a significantly smaller number of electrons, because only a 
certain class of electrons can be scattered Into the loss cone, l.e. 
only trapped electrons with sufficiently small Initial pitch angles are 
precipitated In the Ionosphere. 

From Fig. 6.1 It is obvious that there is always an a /2 

IDSX 
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FIGURE 6.1 GENERAL DISTRIBUTION FUNCTION. Differently shaded areas 
Indicate the various behavior of electrons as they Inter- 
act with the variable amplitude wave. 
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such that electrons witha>ot cannot be scattered into the loss cone. 

max 

As noted above those scattered electrons must have been trapped, l.e 
only trapped electrons may have their pitch angles decreased by the 
amount required for precipitation. Based on the above limits for v„ 
and a it is feasible to define a region in V|, - a space (cross-shaded in 
Fig. 6.1) containing electrons that can be scattered into the loss cone. 
This region in the v„ - a space is further divided into a number of mesh 
points identified by their Vh and a , and this mesh then represents the 
initial distribution. The number of electrons at each mesh point is 
equal to twelve, reflec :lng the fact that electrons are uniformly 
distributed in phase. Figure 6.2a Illustrates the unperturbed 
distribution function; note that we use the number density of electrons 
Ng rather than f(vu ,cx). The number density and f(vn ,a) are related 
through [inan, 1977]: 

Np « 27T f(v,^a) vf, Av„ Aa (6.1) 

^ cos 'a 

Using Eq. 6.1 It Is also possible to find the actual number of 
electrons represented by a single test electron. 

During the Interactions the Initial distribution of electrons 
(Fig. 6.2a) Is perturbed by the wave, and the final distribution Is 
shown in Figure 6.2b. Note that the velocltj^ mesh size is different in 
Figs. 6.2a and 6.2b, since the energy of the electrons tends to be 
signlflcancj.y Increased through the interaction process. Beside an 
overall Increase in electron energies, three electrons are scattered 
Into the loss cone. In the next section precipitation fluxes are 
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FIGURE 6.2 SIMULATION OF THE DISTRIBUTION FUNCTION. (A) The 

unperturbed dlstrlbutlor. (B) Perturbed distribution. 
The numbers In each Individual cell Indicate the number density of 
electrons. 
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computed for three particular cases. 


B. PRECIPITATED ELECTRON FLUX 

Here we compute the precipitated electron fluxes Involving a 
one-sided wave function, and for three different maximum wave 
intensities (Em ■ 50, 150 and 250 yV/m) . The maximum Initial pitch 
angle considered In these calculations Is lO** , since there are no 
electrons with a > 10° scattered Into the loss cone even when the 
electrons Interact with a very strong wave, l.e. En * 250 yV/m. The 
Initial unperturbed number density function Is the same In all three 
examples, and was already shown In Fig. 6.2a. Furthermore, the 
distribution function Is taken as 

f(v,a) • ^ g(a) (6.2) 

where A Is a constant and g(a) Is some function of pitch angle. In our 
calculations g(a) Is assumed to be an Isotvoplc function given by 

g(a) - gi(a) “ 1 a > (6.3) 

- 0 a < Oic 

The following analysis Is similar to that presented by Inan 
Cl 977], although In his work electron scattering was due to 
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gyroresonance Interactions. First, before computing the precipitation. 
It Is feasible to compute the wave Induced pitch angle perturbations 
given by f(a) which Is obtained by Integrating f(vn,a) over the velocity 
range of Interest. In our examples. Involving a 5 kHz wave. It Is found 
that the maximum parallel velocity after the Interaction Is 
^“max" l*8vp(i, whereas the minimum parallel velocity Is 0.98 Vp„. 

The equatorial phase velocity Vpci for a 5 kHz wave Is 11.23 10 m/sec. 
Thus the pitch angle distribution Is given by 

v„ - 1.8 Vp„ 

f(a) " I vfidv,, (6.4) 

v„ - 0.98vp„ 

remembering that electrons are uniformly distributed In Initial phase, 
which results In the factor 2 ^ In Eq. 6.4. 

Figure 6.3 shows the normalized pitch angle distribution f(a) as 
a fuTiCtlon of a for different wave intensities. The dashed curves show 
the Initial unperturbed distributions, whereas the solid curves Indicate 
the final distributions. These results show that the longitudinal 
resonance Interaction requires rather strong waves In order to scatter 
electrons Into the loss cone. For a wave with E» " 50 UV/m 14 pT) 

the perturbations are very small, and only a few electrons are scattered 
below When the wave amplitude Is Increased the loss cone starts to 

fill with electrons, and also electrons with higher pitch angles are 
scattered down to lower pitch angles. This process Is best Illustrated 
In the case of a 250 yV/m wave, where the loss cone Is filled with 
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FIGURE 6.3 NORf^ALIZED ELECTRON DISTRIBUTION f(a). The dashed 
lines represent the unperturbed distribution. The 
solid curves represent the perturbed distribution. 
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electrons having a wider range of Initial pitch angles than the 
electrons reaching the loss cone In the two other cases. 

The total number density of electrons precipitated In the 
velocity range 0,98 Vp„ to 1,8 Vp„ Is given by 


“ic 


N ■ 2tt| I f(vii,a) v« dviida L®(l+3sin^X) 


Vz 


(6.5) 


cos a 


0 0.98vp„ 


where the factor L^(l-f3sln^ X ) accounts for the convergence of the 
field line going from the equator to ionospheric heights. The 


precipitated energy deposition rate Is computed In similar fashion by 

1 


Including the energy weighting factor j m 
yields 


Vn 


cos^a 


In (6.5) which then 


1.8vp„ 

Q - 2ttI I f(v„,a) — l/2m-^- dv„da (l+3sln^ X )^^ (6.6) 

J J cos^a cos*a 

0 0.98Vpn 


The Integrals In Eqs. 6.5 and 6.6 are easily evaluated by a numerical 
integration. For the three examples considered the normalized energy 
deposition rate, defined as Q/A where A Is defined In Eq. 6.2, are: 


Eh 

- 50 uv/m 

Qjj - 0.9652* 10'*“ 

erg/cm^ / sec 

•J 

Eh 

- 150 yv/m 

Qjj - 0.8129* 10’*^ 

erg/cm^/sec 


Eh 

- 250 yV/m 

Qjj - 0.3565*10*“ 

erg/cm^/sec 

- j 
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To compute Che total energy deposition It Is necessary to 
evaluate the constant A. This can be done by computing the total number 
density In el/cc In the specific velocity range 0.98-1.8 Vpii . In 
this case 


TT 1.8Vp„ 


- 


Si 


27T \ \ ^ sina dv da 


0 0.98vp„ 


The above Integral yields 


(6.7) 


A - 2x10® N_ (6.8) 

£i 

Finally, to compute A It Is necessary to estimate from the 

reported measurements. From Schleld and Frank [1970] we find Chat N * 1 

el/cc In the 1-2 KeV range and that the number density varies as v~*^ with 

velocity (E“* with energy) . In our case the electron energies are 

300-1000 eV which results In N " 10 el/cc, since the number density 

£ 

Increases with decreasing electron energy. Substituting N ■ 10 el/cc 
In Eq. 6.8 we find Chat A ■ 2x10®. 

The next step Is to compute the absolute energy deposition rates 
by multiplying the normalized rates by the constant A. The results 
are shown below: 


E„ - 


50 yV/m 


Q ■ 1.94xl0"®erg/cm^/sec 
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E„ ■ 150 yv/m Q ■ 1.66xiO”^erg/cm^/sec 

Eh - 250 yv/m Q ■ 7.40’<10"^erg/cm^/aec 

The above values Indicate that the fluxes precipitated by a 5 
kHz wave, which Is amplified at the equator through the gyroresonance 
Interaction, are rather small, especially when compared to those 
computed for gyroresonance Interactions. Results for the gyroresonance 
process calculated by Inan [1977] Indicate flux levels of 0.01-0.2 
erg/cm^/sec for a 10 pT wave. Note that 10 pT corresponds to Eh * 30 
yV/m for 6 ■ 30” and f ■ 5 kHz. Thus, the scattering efficiency Is 
considerably higher for the gyroresonance than It Is for the 
longitudinal resonance. 


C. ENERGY EXCHANGE AND BALANCE 

From the analytical and numerical studies It Is evident that the 
scattering of electrons Is always associated with energy transfer, l.e. 
If electrons gain energy then the wave Is attenuated, or If electrons 
lose energy then the wave Is amplified. Also, a large scattering Is 
always associated with a large energy exchange. Such behavior 
constitutes another major difference between the longitudinal and the 
gyroresonance processes; namely, electrons can be scattered 
significantly through the gyroresonance Interactions with a very small 
amount of energy transfer. This Is explained by the fact that In 
gyroresonance It Is the momentum transfer that causes pitch angle 
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changes, whereas Che energy remains almost unchanged [Inan, 1977]. 

The cocal energy balance calculaClons for Che longitudinal 
resonance process are extremely complicated as they involve a large 
number of electrons. As indicated in Fig. 6.1 Che electrons with 
parallel velocities close to the wave phase velocity become trapped 
which then results In scattering from -0.2** up to -6^ for pitch angles 
from 10** to 70°, respectively. The scattering of untrapped electrons Is 
smaller and positive, about 0.05-0.1° on the average. However, only a 
fraction of the electron population becomes Crapped, while Che number of 
untrapped electrons Is much larger. From the sample calculations It 
was estimated that the upper velocity limit for untrapped electrons can 
be as high as 1.30 Vp„ , l.e. even If the Initial parallel velocity of 
the electron Is v„ ■ 1.30 Vp„ the electron is still scattered more than 
±0.005°. The scattering of ±0.005° represents a practical threshold of 
resolution for the numerical Integration method used In our simulations. 
This resolution limit was found by setting Em ■ 0 UV/m, l.e. computing 
only the adiabatic motion of the electrons and comparing the Initial and 
the final pitch angles. Theoretically, the difference between these two 
pitch angles should be zero, whereas Che numerical results have shown 
10.005 fluctuations, which are than used as the limit of accuracy 
(resolution) . These fluctuations are primarily due to the Integration 
scheme, which uses linear interpolation. Returning to the energy 
exchange problem. It Is evident that both trapped and untrapped 
electrons play important roles, and it Is rather difficult to find an 
exact solution to this problem as the muaber of electrons Involved Is 
very large. 
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However, Ic Is possible to estimate the energy transfer as 
follows; let us consider the example of Fig. 5.7a (solid curve) showing 
scattering as a function of the Initial parallel velocity for a fixed 
Initial pitch angle a ■ 10°. This curve can be replutted substituting 
energy changes for pitch angle changes and also expanding the velocity 
range. Note that these results must be weighted by an appropriate 
function to account for different number densities at different 
velocities. This weighting function Is assumed to have a v~^ 
characteristic (Eq. 6.1). Figure 6.4 shows both unweighted and weighted 
energy transfer, l.e. the average energy gain (loss) per electron with a 
given Initial parallel velocity, as well as the weighting function 
(dashed curve). Now It Is possible to use a numerical Integration to 
estimate the total energy balance for this particular case. 

The total energy exch>\nge Is given as 

£ 2 ( 1 .03Vpii) 

Ae » I dE (6.9) 

Ei(0.99vp„) 

where AE represents the total energy exchanged through the longitudinal 
Interaction with electrons whose Initial parallel velocities are In 
0.99-1.03 Vpii range, and all those electrons have the same pitch angle 
a ■ 10°. The quantity dE gives the weighted amount of energy exchanged 
per electron at a particular parallel velocity, and It Is shown In Fig. 
6.4. The final result of the above Integration Is AE " 0.03 eV. Though 
this number Is obtained using only twelve electrons It Is evident that 
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FIGURE 6.4 ENERGY TRANSFER FROM WAVE TO ELECTRONS FOR THE CW SIGNAL INTERACTION. 

The energy transfer is shown as a function of electron parallel velocity, 
and positive <AE> Indicates that electrons gain the energy. 


128 


the totel energy exchange at the particular pitch angle Is very small 
even when the actual number of electrons Is much larger. 

To compute the overall energy balance similar calculations 
should be done for other pitch angles. A rough estimate using Figs. 5.6 
and 5.7 Indicates that the total energy transfer Is very small, since 
the positive and negative scattering cancels out, l.e. the total area 
underneath <Aa^q> curve is approximately zero. 

Summarizing, It appears that both the precipitation and wave 
amplitude amplification (attenuation) for our sample case are small. 
Thus, It may be very difficult to observe the presence of this type of 
longitudinal Interaction using ground observations. Another possibility 
for detection would be to use satellite borne particle detectors and 
to measure a relatively sharp depletion or electron density around 
Vii ■ Vpii . However, the problem Is that particle detectors measure 
energies and pitch angles rather then parallel velocities and pitch 
angles. Note that the problem arises from the fact that the narrow 
range of parallel velocities which are affected (and wide range of pitch 
angles) maps Into a wide range of energies (and pitch angles). 

For example, if the parallel velocity equals the phase velocity, 

Vii ■ Vp„, and pitch angles vary from 5“ to 70®, the corresponding 

electron energies vary from E to E (l+tan^70®)/(l+tan^5*) • 8.48 E , 

o 0 o 

where E Is the energy of the electrons with 5** pitch angle. Beside the 
o 

above mentioned spreading effect, which tends to dilute the effects of 
the longitudinal resonance when measured on an energy basis, the 
particle detector resolution Itself may pose a problem. The typical 
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resolution of particle detectors Is about 2.5*’- 5*’ In pitch angle, and 
about 15% In where Is the energy of Interest. For example, if we 
want to measure the density of electrons with energy Ejjj - 2E^, and pitch 
angle a ■ 45“, the corresponding resolution cell would be as shown in 
Figure 6.5. On the other hand, the longitudinal resonance will tend to 
remove electrons from a narrow velocity band around Vpii, leaving a 
depletion region In the distribution (Fig. 6.5). The width of the 
depletion region is very email, so tlu t it occupies only about 30% of 
the resolution cell, s indicated In Fig. 6.5. Therefore, even if we 
remove all of the electrons from this depletion region, the particle 
detector would see only a 30% decrease in the number of electrons within 
th^ resolution cell. We recall from Chapter V that longitudinal 
resonance interactions. Involving moderate amplitude waves, result in 
trapping of only about 30% of the electrons that satisfy the resonance 
condition (we considered only the trapped electrons, because only those 
electrons undergo sufficient change in vn to be moved from one 
resolution cell to another) . Thus the maximum total depletion factor 
for the resolution cell is estimated to be about 10%. On the other 
hand, typical particle detector measurements (e.g. Kimura at al., 1982) 
indicate large temporal variations of the electron flux, approaching r.n 
order of magnitude in intervals as short as 50 sec. Fo‘ that reason the 
particle detector sensitivity is reduced, because it becomes very 
difficult to distinguish between variations due to spatial changes in 
particle distribution and wave Induced variations. Thus, present 
particle detectors are probably not capable of detecting perturbations 
of the electron distribution due to longitudinal resonance interactions. 


VII. APPLICATIu:i3 TO MAGNETOSPHERIC PHENOtffiNA 


Although It was found that the scattering efficiency of the 
longitudinal resonance process Is small. It Is possible that the 
bunching effects of the process may have Important magnetospherlc 
application <. In this chapter we consider applications of the 
longitudinal resonance to the generation of whistler precursors and to 
the generation of broadband VLF hiss. We also discuss the size of the 
Internal electric field created in the bunching process. 


A. GENERATION OF WHISTLER PRECURSORS 


Whistler precursors are discrete rising tone emissions that 
precede two-hop whistlers, starting shortly (0.1-0. 3 sec) after the 
one-hop delay. The precursor may consist of one or more discrete 
emissions. For the particular measurements of August 2, 1973, the 
number of emissions varied from one to seven. Figure 7.1 Illustrates 
three typical cases of precursors showing both one-hop whistlers 
(recorded at Slple, Antarctica) , and precursors with corresponding 
two-hop whistlers (recorded at Roberval, Canada). There Is no precursor 
In Fig. 7.1b, Illustrating the fact that not all whistlers propagating 
on the same path trigger a precursor. Figure 7. Id depicts a single 
emission precursor, while Fig. 7. If shows a multi-emission precursor. 
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FIGURE 7.1 SPECTROGRAMS OF WHISTLER PRECURSOR EVENTS RECORDED AT 

SIPLE/ROBERVAL CONJUGATE STATIONS. The causative spheric 
Is marked vlth an arrow, and the whistler component which triggers the 
precursor is marked by a W. (b) shows no precursor, (d) shows a single 
emission precursor, and (f) shows a multi-emission precursor event. 
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FIGURE 7.2 EXPANDED SPECTROGRAM OF THE PRECURSOR AT 1400 UT FROM 
FIGURE 7.1. (b) shows the corresponding amplitude 

variation In a 300 Hz bandwidth, and (c) Indicates the rate of 
frequency change of the frequency- tracking filter used. 
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These particular data were analyzed by Park and Helllwell [1977], and it 
was found that the precursors were triggered only by the whistlers 
propagating in one particular duct, and that the precursors themselves 
propagated In the same duct. The duct parameters were L " 3.6 and 
equatorial electron density n^^ 440 el/cc. The center of the 

plasmapause was located at about L ■ 4.2 where the equatorial electron 
density dropped by factor of ten. Figure 7.2 shows an expanded 
frequency-time spectrogram of the precursor at 1400:03 UT, along with 
amplitude and frequency changes measured using a frequency- tracking 
filter. The growth rate deduced from that figure Is about 105 dB/sec, 
and the rate of frequency change Is about 6.5 kHz/sec. 

Park and Helllwell [l977] have reviewed different proposed 
generating mechanisms for precursors. Including the hybrid mechanism 
suggested by Helllwell [l965] and Dowden [1972]. This Is based on the 
presence of hybrid whistlers, which first propagate In the earth- 
ionosphere waveguide to the conjugate hemisphere and than return through 
the magnetosphere and trigger precursor emissions. Other mechanisms 
Include one pro\>osed by Reeve and Rycroft [l976] In which the nonducted 
whistler is reflected In the conjugate hemisphere at the lower hybrid 
resonance (LHR) frequency, and is then deflected by the plasmapause such 
that It enters the duct near the equator, triggers the precursor through 
the gyroresonance, and then leaves the duct. A third mechanism 
Involving a nonlinear multiple wave Interaction known as parametric 
decay has been suggested by Reeve and Boswell [l976]. 

Considering various precursor mechanisms for the Aug. 2, 1973 
case, the hybrld-whlstler hypothesis can be Immediately excluded because 


-«■- . 


135 


there was no evidence of hybrid whistlers. The mechanism suggested by 
Reeve and Rycroft [l976] requires special propagation conditions which 
are difficult to apply to multi-component precursors with a wide range 
of starting frequencies (-1 kHz for the example shown in Figure 7. If). 
Furthermore, the L-shell values of the duct and the plasmapause differed 
by more than the 0.15 required by their model. Finally, the parametric 
decay mechanism cannot explain the multicomponent precursors; hence Park 
and Helllwell Cl977] have suggested a new mechanism. 

The new mechanism Is Illustrated In Figure 7.3 and Its time 
sequence Is described below: 

a) A lightning Impulse In the northern hemisphere produces a 
whistler propagating toward the equator. 

b) The whistler wave train signal and the energetic electrons 
streaming toward the equator Interact with one another through the 
longitudinal resonance process. 

c) Due to the longitudinal Interaction, electrons become space 
bunched, which then temporarily Increases the electron flux within a 
certain range of parallel velocities. 

d) This enhanced electron flux reaches the equator while the 
whistler signal that caused the bunching continues to travel toward the 
southern hemisphere. 

f) After crossing the equator the enhanced electron flux 
Interacts with northward traveling power line harmonic (PLH) waves 
through the gyroresonance process. The enhancement of the electron flux 
is sufficient to lower the threshold of this interaction below the level 
required for triggering of an emission by one or more lines of PLH 
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FIGURE 7.3 SCHEMATIC ILLUSTRATION OF THE WHISTLER PRECURSOR 
GENERATION MECHANISM. See the text for details. 
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waves. These emissions travel toward the northern hemisphere. 

g) While the triggered emission (precursor) travels toward the 
northern hemisphere^ the one-hop whistler reaches the conjugate point In 
the southern hemisphere, where It Is reflected. It then travels back to 
the northern hemisphere. 

h) The precursor reaches the northern hemisphere followed by a 
two-hop whistler, resulting In a frequency- time spectrograms similar to 
those depicted In Fig. 7.1. 

The detailed timing of this process was worked out by Park and 
Helllwell L^977] and It was shown that this mechanism can explain 
different properties of the Aug. 2, 1973 precursors such as variable 
starting frequency, multicomponent emissions and variable starting time. 
However, there are some special requirements that have to be met In 
order for this mechanism to work. First, the enhancement of the 
electron flux achieved through longitudinal resonance must be large 
enough and should last about 200 ms, so as to provide both the threshold 
for triggering through gjnroresonance as well as the temporal growth time 
required for emission generation. Second, the PLH waves (which 
obviously must be present for this mechanism to work) must have 
amplitudes such that they approach the triggering threshold level. 

PLH activity appeared from time to time In the August 2, 1973 
case; during some Intervals It dominated the VLF spectrum. Park and 
Helllwell [1977] found that the PLH propagated In the same duct with the 
precursor; this suggests that PLH waves were present at the time of the 
precursor observations and, when not detected, were probably close to 
the threshold for triggering emissions through cyclotron resonance. 
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As already stated the gyroresonance triggering mechanism will 
work only If the electron density perturbations achieved through the 
longitudinal resonance result In an electron flux Increase which lasts 
at least **200 ms. The 200 msec requirement is associated with a typical 
temporal growth time [Stlles and Helllwell. 1977], l.e. a typical delay 
from onset of temporal growth to emission triggering. This flux 
Increase can be achieved. In principle, through electron bunching. We 
have shown in Chapter V that the longitudinal resonance Interaction 
results In significant space bunching, which In our particular case of a 
monochromatic signal was about 180%, l.e. the electron density was 
enhanced roughly by factor of two at v» ■ with the density 

enhancement decreasing for other parallel velocities. 

However, In order to explain multi-component precursors It Is 
necessary to Increase the electron flux over a relatively wide range of 
parallel velocities. At each velocity the flux increase should last for 
about 200 msec. To Illustrate this process we consider a multi- 
component precursor consisting of two emissions with starting 
frequencies fj * 2 kHz and f 2 * 3 kHz, and assume that those emissions 
are triggered at the equator, although the triggering location must be 
slightly off the equator to account for the rising frequency-time 
characteristics. From the gyroresonance condition at the equator 
f(l+ Viigq/Vp^^gq) » f|| the parallel velocities at which the flux must be 
Increased are ■ 76.6 10® m/s and Vi,gq 2 " 57.1 10® m/s, where we 

used fgeq* 18.7 kHz and fpgq * 188.8 kHz. Thus the whistler Interacting 
with the energetic electrons must be able to produce an Increased flux 
at those two velocities for 200 msec. We also recall from Chapter III 


that the parallel velocities vnj and Vn 2 vary along the field line as 
indicated in Fig. 3.11, and that the electrons with higher pitch angles 
mirror closer to the equator. 

Next we recall that the longitudinal resonance condition is 

V V 

given aa vn ■ Vp^^ ■ c f - f) » which yields the resonance 

frequency f ■ 1/2 ^fn±Cf| ” 4(vtifp/c)^ (the plus sign gives f 
where the waves become unducted» so we can disregard that solution). 

The resonance frequency changes as we change the parallel velocity. 

For example, if we consider electrons with V'leq^ ^"eq2 assume 

a ■ 10°, their parallel velocities at 50° latitude are Vu^* 0.30 
22.9 10^ m/s and vh^" 0.30 corresponding 

resonant frequencies are fj^ ■ 2.65 kHz and f 2 ■■ 2 kHz. Thus a whistler 
train of appropriate frequency range can Interact with electrons with 
different parallel velocities, such that when those velocities are 
mapped back to the equator they satisfy the gyroresonance condition at 
different frequencies. If the perturbations of the electron flux at 
those different velocities are large enough and last long enough ("200 
msec), they could result in emission triggering at those frequencies. 
This would then provide a basis for explaining the generation of 
multi-emission precursors. 

We want first to Illustrate that the flux perturbation at a 
given parallel velocity (actually in a narrow range of about 1% around 
that velocity) can last longer than 200 msec. In order to do that we 
recall the results for the interaction with a spatial pulse from Chapter 
V. From Figs. 5.17, 5.18, 5.19, and 5.20 we see that a 1000-km-long 
spatial pulse can trap electrons in a narrow band of velocities (-21), 
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and chac those electrons beside being trapped, l.e. space bunched, 
undergo pitch angle scattering on the order of a few tenths of a degree. 
Although this spatial pulse Is stationary and monochromatic, the results 
from that analysis can be related to the whistler train If we consider 
the whistler train to be composed of segments of approximately constant 
frequency, ffe consider one of those segments with frequency f • 2 kHz; 
the group velocity of that segment at SO" latitude (L ■ 3.6) Is about 
30 10‘ m/s , and If It Interacts with electrons for about 2000 km (this 
Is comparable to the length of the spatial pulse considered In Chapter 
V) the total Interaction time Is about 70 msec* On the other hand, as 
long as an electron Is trapped It does not matter If the trapping signal 
Is a stationary amplitude pulse (not moving along the field line) or a 
moving segment of a whistler. If the length of the Interaction region 
In the two cases Is comparable, the trapping and scattering effects 
should also be comparable. 

This segment of the whistler Is therefore capable of Increasing 
the flux in a narrow band of parallel velocities, but this Increased 
flux should last at least 200 msec at the equator In order to provide 
the basis for emission triggering. The total duration of the flux 
perturbation depends on the latitude at which the resonance takes place, 
and on the pitch angle of the electrons Involved. For example. If we 
want the triggered emission to start at 3 kHz It Is necessary to 
Increase the electron flux In a narrow band of velocities around 
v„ - as noted above. However, electrons with Vi,« will have 

different pitch angles at the equator, and will thus mirror at different 
latitudes (see Fig. 3.10). For a " 10® the mirror point is at 53® 
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latitude, while for a ■ 50® the mirror point is at 20® latitude. Thus 
our whistler segment at 2 kHz, as it travels toward the equator (from 
higher latitudes toward lower latitudes), first encounters electrons 
with a 10® at about SO® latitude (the time of this encounter is the 
reference time t <■ 0). As noted earlier, if the interaction lasts for 
about 70 msec, it should be sufficiently long time to bunch the 
electrons. During those 70 msec both wave and electrons move from about 
50® to about 48® latitude* After the interaction is over it takes about 
0.43 sec for the bunched electrons to reach the equator, or essentially 
the travel time from 48® latitude to the equator. When the electrons 
arrive at the equator they have Vh ■ C''® have neglected the 

parallel velocity changes due to the interaction, as it is assumed that 
the scattering is small). Furthermore, as our whistler segment gets 
closer to the equator it Interacts with electrons with progressively 
higher pitch angles. The arrival time at the equator for those 
electrons with higher pitch angles can be calculated using the above 
described method. For a 50® the interaction occurs at 20^ latitude, 
and those electrons arrive at the equator at t = 0.69 sec (0.5 sec for 
whistler travel time from 50® to 20** latitude, ‘'0.1 sec for the 
interaction, and 0.18 for particle transit from 20® to the equator). 

Thus the perturbation at the equator irould last about t ■ 0.69 - 0.43 ■> 
0.26 sec, which is sufficient for the development of emission 
triggering. Computations for the whistler segment with f " 2.65 kHz 
indicate that the corresponding flux perturbation lasts about 210 msec. 
Therefore it is found that the electron flux perturbation may last long 
enough and may cover a sufficiently wide range of parallel frequencies. 
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Noce chat similar computaClons were done by Park and Helllwell [l977], 
but wlChour consideration of the Interaction time. 

As noted earlier In Chapter V, this perturbation (space 
bunching) Is associated with an amplitude threshold of the waves driving 
Che longitudinal resonance. This suggests that one could measure the 
amplitudes (on the ground) of whistlers with and without precursors, and 
therefore test for the presence of the threshold. Such amplitude 
measurements were made on one-hop whistler » recorded at Slple, 
Antarctica, and propagating at L ■ 3.6 on August 2, 1973. The data were 
taken at two frequencies, 4000 Hz and 4600 Hz, using a bandpass filter 
with Af ■ 300 Hz. This provided the temporal resolution needed to 
distinguish a particular whistler component connected with precursor 
generation from other multipath components. The results of those 
measurements are shown in Figure 7.4 as amplitude vs. time diagrams. 

The whistlers without precursors are indicated by crosses, the whistlers 
with single emission precursors are indicated by circles, and the 
whistlers with multicomponent precursors are Indicated by squares, where 
the numbers above the squares represent the number of individual 
emissions forming a single precursor event. 

Figure 7.4 shows that the amplitudes of the one-hop whistlers 
decreased, on average, from -15 dB (0 dB level corresponds to 100 pV/m) 
to about -22 dB for f ■ 4600 Hz. For f - 4000 Hz the average amplitude 
decreased from -13 dB to about -17 dB in the same period of time between 
1335 UT and 1415 UT. This overall decrease of the whistler amplitudes 
is most likely a result of increased absorption in the ionosphere 
because of transition from nighttime to daytime conditions (sunrise time 
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was around 1400 UT) * Helllwell Cl96S] has shown ChaC there Is a 
significant Increase in the ionospheric absorption at VLF for the 
night-day transition, aad that the amount of the absorption increases 
rapidly with increasing frequency. This prediction is consistent with 
the above observations; the amplitude level at 4000 Hz drops about 4 dB, 
whereas the amplitude level at 4600 Hz drops about 7 dB. If we further 
assume that the luaxlmum ionospheric absorption occurs in the D region at 
about 100 km altitude CHelllwell, 1965] it is possible to estimate the 
duct exit point using the path L value as one coordinate and sunrise 
time at 100 km altitude as the recond coordinate. From Fig. 7.4 we see 
that the amplitudes of the whistlers start to decrease around 1355 UT 
which is then assumed to indicate the beginning of sunrise effects. On 
the other hand calculations show that for sunrise times of 1355 UT and 
1405 UT at 100 Ion altl . ude, the terminator reaches the latitudes of 71**S 
and 72'*S, respectively. This period of time (1355-1405 UT) is the time 
when the whistler amplitudes are rapidly decreasing (Fig. 7.4), 
suggesting that the latitude of the whistler duct exit point was between 
7l”S and 72°S. Because the whlstl^ir duct was on L ■ 3.6, we can find 
where this line Intercepts the above latitudes; the result is shown in 
Figure 7.5. The estimated location of the duct exit point lies in the 
north-west direction from Slple Station, at a distance of about 490 km 
for 71*’S latitude, and about 830 km for 72°S latitude. 

A more Important feature of Fig. 7.4 is the presence of a 
threshold level that a whiscler amplitude must exceed in order to 
trigger a precursor. This amplitude threshold is most clearly seen 
between 1335 and 1350 UT. As found earlier in Chapter V, such behavior 
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FIGURE 7 5 ESTIMATED LOCATION OF THE WHISTLER DUCT blT POINT FOR THE PRECURSOR EVENTS 
OF AUGUST 2. 1973. 
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Is one of Che characteristics of the longitudinal resonance interaction, 
which then supports the precursor generation mechanism suggested by Park 
and Helliwell [l977]. We note Chat the apparent gap in the precursor 
activity between 1350 and 1400 UT is artificial. At least five 
precursor events were observed at Roberval, but it was not possible to 
measure the corresponding amplitudes of the one-hop whistlers due to the 
operation of a VLF transmitter at Slple (receiver preamplifier muted). 

In the next period of time, between 1400 and 1415 UT, Che 
precursor activity still exhibited a threshold, although not as clearly 
as before. The presence of many multicomponent precursors indicates 
favorable triggering conditions for the gyroresonance interaction 
between electrons and PLH waves. This is supported by the level of 
spontaneous magnetospherlc emissions, which increased sharply around 
1400 UT, and strong PLR (power line radiation) which was observed for a 
period of a few minutes. 

The data show that the precursor generation was associated with 
an amplitude threshold in the driving whistler, but the model suggested 
by Park and Helliwell [ 1977] also requires that the space bunching 
produced by the one-hop whistler be sufficient for triggering emissions. 
As it was found earlier, the space bunching process can roughly double 
the electron density (flux). According to Helliwell and Inan [l982] who 
proposed a feedback model to explain VLF growth and discrete emission 
triggering in the magnetosphere (through gyroresonance), a doubling of 
the electron flux is usually sufficient to result in the triggering of 
emissions. In their model the loop gain G is directly proportional to 
the electron flux. For G<1 the system acts like an amplifier, while for 
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G >1 Che system becomes unstable and can generate emissions. Therefore, 
a doubling of Che flux could easily boost Che loop gain G to a value 
larger than unity and thus result in triggering. 

Thus the precursor generating mechanism suggested by Park and 
Helllwell [1977] appears to be supported by the results found for Che 
longitudinal resonance. Including both the amplitude threshold and the 
level of Che density bunching. 

In the next section we discuss some other aspects of the 
longitudinal resonance interaction that may be important In other 
magnetospherlc processes. 


B. VLF HISS 


One of many magnetospherlc processes for which the generating 
mechanism Is not certain Is VLF hiss, most often observed on the ground 
as relatively broad band (several kilohertz) noise. VLF hiss often 
shows no discrete structure, having the appearance on a spectrogram of 
band-limited white noise. This type of spectrum is characteristic of 
auroral and plasmaspheric hiss, whereas mid-latitude hiss usually shows 
some kind of discrete structure. Therefore, Che hiss generating 
mechanism must be such that it can explain the generation of relatively 
wideband signals, and also account for the observed amplitudes of such 
signals. 

An electron propagating in a dielectric medium does not radiate 
as long as its velocity remains less than the phase velocity in that 
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medium; If the electron velocity Is larger than the phase velocity we 
have a case of Cerenkov radiation. The two situation are depicted in 
Figure 7.6, and we note that the electron radiates at only one angle 
when v„ > c/J^, However, In the case of a dispersive medium different 
frequencies are radiated In different directions, as shown In Figure 
7.7. In the magnetosphere the radiated frequencies are within the VLF 
range. Thus If the amplitude of the Cerenkov radiation Is large enough 
It could account for the hiss generation. It should be noted that the 
condition for Cerenkov radiation Is exactly the same as the condition 
for longitudinal resonance, l.e. the electron velocity must match the 
phase velocity (In the direction of electron travel) In a particular 
medium. 

In the magnetospherlc case It can be shown that there are In 
general two Cerenkov frequencies radiated at each angle, and that the 
radiation condition Is not met when the parallel velocity exceeds the 
critical velocity v„ [Brice, 1964]. The critical velocity corresponds 
to propagation In the Gendrln mode, which was defined In Section III.B. 
As noted earlier, the broadband nature of Cerenkov radiation makes It 
Interesting as a possible source of VLF hiss, and It was considered by 
many authors [Ellis, 1959,1960; Dowden, 1960; McKenzie, 1963; Llemohn, 
1965; Mansfield, 1967; Seshadrl, 1967; Jorgenson, 1968; Llm and Laaspere, 
1972; Taylor and Shawhan, 1974]. However, all of the power density 
calculations fell short of explaining the observed power density of VLF 
hiss. Indicating that Incoherent Cerenkov radiation Is not sufficiently 
strong to account for VLF hiss. For this reason other mechanisms were 
suggested which are still based on the Cerenkov radiation, but In which 
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radiation is either coherent [Taylor and Shawhan, 1974], or amplified 
through interaction with an electron beam [Swift and Kan, 1975; Maggs, 
1976] . In the case of the coherent radiation it is assumed that the 
radiation from n electrons is in phase, resulting in F > n^P, where P is 
the power radiated by each electron. On the other hand, if all n 
electrons radiate incoherently (random phase) the total radiated power 
is given by P ■ nP. 

Due to the n^ dependence, a relatively small number of electrons 
radiating coherently could produce power levels which are in agreement 
with the meaudi'ements. Thus the problem is to identify a process that 
could result in electron bunching such that the bunch dimensions are 
much less than a wavelength (smaller dimensions mean greater coherence). 
As already shown, the longitudinal resonance interactions may produce 
such bunches of electrons, and It may be speculated that the radiation 
coherence needed to explain VLF hiss Is created In the following way: 

(1) first a strong signal bunches a significant number of electrons 
(stronger waves would produce better coherence), and (11) the bunched 
electrons become detached from the bunching wave. The detachment may be 
due to difference In phase and group velocity, as Is the case for the 
whistler mode where the phase and the group velocity are always 
different (except for f > fg/2) . For example, consider a pulse with 
f < %/2 so that Vg > Vp,, . Electrons trapped by this pulse will have 
^11 ■ ^pii » because the wave energy propagates with Vg > Vi, , those 
electrons slide backwards through the pulse, and eventually emerge from 
the tall end of the wave packet. Such a blob of electrons could radiate 
coherent Cerenkov radiation. 
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However, It remains to be seen how long this blob of electrons 
remains bunched, because It may contain electrons with different pitch 
angles and different parallel velocities. For the moment let us assume 
that all electrons have the same parallel velocity, but different pitch 
angles which means that they have different variations of parallel 
velocity as required by their adiabatic motion. Thus, for a given 
spread In pitch angle It may be determined how long It takes the 
separation between the low and high pitch angle electrons to become 
larger than the wavelength, which than destroys the radiation coherence. 
The sample calculations have shown that the coherence time for a given 
Initial spread In pitch angles depends strongly on the latitude where 
the electrons become detached from the bunching wave, l.e. on the 
latitude at which their motion begins to be entirely governed by the 
static magnetic field. For example, assuming the Initial range of pitch 
angles to be from a " 10** to a ■ 20**, and detachment at 20** latitude 
(electrons are moving toward the equator). It takes only about 1 msec 
before the separation between 10** and 20** electrons becomes larger than 
one wavelength. On the other hand. If the detachment occurs at 
latitude (for the same Initial range of pitch angles) It takes about 0.2 
sec for the same process to occur. Note that after 0.2 sec the 
electrons reach 4^ latitude, but on the other side of the equator. 

A blob of electrons created through the longitudinal resonance 
Interaction and with a spread In pitch angle only could radiate 
coherently for a substantial period of time (few tenths of second). 
However, the electrons within a blob have slightly different parallel 
velocities, e.g. a t 3 rplcal spread In parallel velocity Is about 400 
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km/ sec (Figs. 5.10 and 5.15). Thus It will take only about t - 2/400 ■ 

5 msec for those electrons to become separated more than a wavelength at 
the equator, assuming the wavelength to be 2 km at the equator. From 
this result It Is evident that spreading due to the finite range of 
parallel velocities occurs much faster than the spreading due to a 
finite range of pitch angles, and that the life time of the blob Is 
about one hundredth of a second. We also note that the blob of 
electrons could further be dispersed due to Interaction with other 
waves. 

Thus it is possible that the short life time during vdilch the 
blob can radiate coherently, together with the fact that there may not 
be many electrons within a single blob, makes the radiated power level 
Insufficient to account for the observations. However, there could be 
more than one blob formed through the above described process, which 
could further enhance the radiation (as long as the radiation from 
different blobs does not Interfere). Even stronger radiation effects 
could probably be achieved If the velocity of the electron blob equals 
the critical velocity, because In that case all radiated frequencies 
satisfy the Gendrln condition given In Chapter II. The enhancement of 
radiation Is expected because for the Gendrln mode the ray direction Is 
field aligned for all radiated frequencies, and the group velocity Is 
Independent of the wave frequency so that wave packets radiated at 
different frequencies travel together [Helllwell, private communication]. 

Another explanation for VLF hiss generation Is based on 
amplification of Incoherent Cerenkov radiation through the wave-beam 
Interaction where the beam provides for the 'bump-on-tall* distribution. 
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As mentioned earlier, a distribution function which has a positive 
slope, as Is the case for the bump-on-tall distribution, may result in 
Landau growth* 


C. COMMENTS ON THE INTERNAL FIELDS OF THE BUNCH 

At this point we should note that space bunching always gives 
rise to an Internal electric field through the Poisson equation. This 
electric field will then act to debunch the electrons, as It opposes the 
wave bunching field. Although this effect can be neglected In test 
particle simulations where the number of electrons Is small, it may 
become Important depending on the actual flux of particles. We have 
shown that significant bunching occurs for a parallel electric field 
around 50 yV/m and higher, so that we choose 5 yV/m as the limit for the 
Internal field, l.e. we ass\ime that Internal fields up to 5 yV/m do not 
significantly affect the bunching process. Using the 5 yV/m field we 
can find an electron density N that Is needed to produce that field. In 
Chapter IV we showed how the twelve test electrons are uniformly 
distributed In phase before the Interaction, and In Chapter V (Fig. 

5.15) we showed that the same electrons are compressed In phase space, 
l.e. space bunched. The typical compression Is about 90^ In phase, or 
500 m assuming X > 2 km. 

At the same time each single test electron actually represents a 
large number of electrons In the real distribution, l.e. each test 
electron represents a sheet of electrons* Thus the question Is, If we 
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have twelve Initially equidistant sheets of electrons, and we displace 
those sheets so that the total displacement Is 500 m, what Is the 
maximum electron density for which the Internal field (due to the 
compression of the sheets) does not exceed 5 yV/m? It turns out that 
this computation Is rather simple, and the electron density is given as 
[Buneman, 1980] 


N 


e As 


(7.1) 


where E is our maximum allowable internal field (negative), and As the 
total displacement of the sheets. Using E ■> 5 yV/m, As ■ 500 m, and 
Cg" 8.854 10 we find N * 0.55 el/m’ which is the maximum allowable 
density, i.e. densities larger than this produce internal fields 
stronger than 5 yV/m, which can reduce the bunching effects. When the 
density of the electrons is knotm we can relate it to the electron flux 
as discussed below. 

It was shovm that trapping occurs in a narrow range of parallel 
velocities centered around the wave phase velocity, so we use 1% as a 

typical value. The next step is to compute the actual number of 

electrons in that velocity range, and then to compare with the 
previously computed N ■ 0.55 el/m’ . The electrons are assumed to have 
an Initial energy of 300 eV and a ■ 10® , so that the corresponding 

parallel velocity is v„ « 9.654 10® m/s. In that case the total number 

of electrons, within 1% velocity range around v,„ is given as (assuming 
an isotropic distribution in pitch angle) 
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N 


ir V 2 



0 vj 


sin a 


dv da 


(7.2) 


where A Is a constant that can be deduced from the flux. It can be 
shown [Inan, 1977] that for E ■ IkeV and a - 90®, A ■ 2 <|), where <j) is 
the differential energy spectrum for 1 keV electrons with a ■ 90° . Note 
that this relationship between and A holds only for a v''^ distribution, 
and It is necessary to use a different relation for other distributions. 
Thus, substituting for A In Eq. 7.2, and Integrating we have (n ^ 3) 


N - 4ttA 


y-n+2+1 

n-2-1 


V2 


Vl 


(7.3) 


whereas for n ■ 3 we have 


N 


" 4 TTA In V 


(7.3a) 


and Table 7.1 shows the results for various values of the differential 

flux (j) (1 keV, a ■ 90**) and various values of n (the constant A Is given 

as where m Is the electron mass), 

m 

Thus, from Table 7.1 we can find the values of n and for which 
the electron density Is lower than 0.55 el/m^, l.e. we see when It Is 
possible to have bunching without creating a strong Internal electric 
field which may significantly decrease the bunching effects. Also note 
that only the trapped electrons contribute to the Internal field. 
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Flux 

cm~^ sr“* s~* kev) 

n 

A 

i 

N (el/m*) 

10® 

3 

1.3 

10“^ 

1.3 10"* 

10* 

4 

2 

10* 

21658 

10* 

5 

2.9 

10” 

3.1 10* 

10* 

6 

4.4 

10* » 

6.8 lO” 

10“ 

3 

1.3 

10"” 

1.3 10"^ 

10“ 

4 

2 

10“ 

2.17 

10“ 

5 

2.9 

lo” 

31 

io“ 

6 

4.4 

10*“ 

6.7 10® 

10^ 

3 

1.3 

10"” 

1.3 10"* 

10^ 

4 

2 

10* 

0.02 

10^ 

5 

2.9 

10” 

0.31 

10^ 

6 

4.4 

10** 

6.7 10** 


TABLE 7*1 Total number of electrons within 1% velocity 
bandwidth for 300 eV electrons as a function of flux and 
various distribution functions. 

Because most of the flux measurements are made at higher energies the 
exact fluxes and distributions at lower energies are uncertain, but as 
those data become available Table 7.1 can be used as a guide to 
determine If the bunching of the electrons Is affected by the Internal 
fields. Present measurements Indicate that the flux can be on the order 
of 10^ to 10^ , and the exponent n can vary between 3 and 5 [Klmura, 
1982; Shield and Frank, 1970]. 

We have presented two examples In which longitudinal resonance 
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interac''i.ons may play an Important role, along with an analysis of the 
limiting electron flux for the bunching. We conclude our discussion 
with a suimna..:y and suggestions for future work. 


VIII, CONCLUSION AND SUGGESTIONS FOR FUTURE WORK 


A. SUMMARY 


We have analyzed the nonlinear longitudinal resonance 
interactions between energetic electrons and coherent VLF waves in the 
magnetosphere. The longitudinal resonance, which may result either in 
wave growth or wave damping, and also causes space bunching of energetic 
electrons, was numerically simulated using time averaged nonlinear 
equations of motion. The simulations were done for single electrons, 
sheets of electrons, and a full distribution of electrons. Those 
studies, done for different types of wave functions, have shown how the 
the wave forces modify the electron trajectories, and that the 
trajectory perturbations result in nonlinear pitch angle scattering. 

The nonlinear pitch angle scattering variations have been studied for a 
wide range of the initial pitch angles, wave amplitudes, cold plasma 
densities and wave normal angles. It was found that there are two basic 
groups of electrons, trapped and untrapped, where the trapped electrons, 
in contrast to the untrapped electrons, are trapped in the potential 
well formed by the *wave. The trapped electrons cause the space bunching 
which increases the electron flux at certain parallel velocities. 

The nonlinear scattering for the longitudinal resonance is found 
to be much smaller compared to that for the gyroresonance interactions, 
indicating a higher efficiency for the gyroresonance process. This is 
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so because the scattering for gyroresonance Is achlevei. - .‘i. ugh the 
conversion of perpendicular momentum of the electron Into parallel 
momentum with very small energy exchange between the wave and electrons, 
while the scattering for the longitudinal resonance Is solely based on 
the energy exchange. Due to the smaller scattering efficiency a full 
distribution simulation produced only small precipitated fluxes, l.e. 
for moderate strength VLF waves the precipitation due to the 
longitudinal Interactions Is below the detectable level of about 0.01 
ergs/cm /sec. 

In a study of magnetopsherlc applications we found support for a 
mechanism proposed by Park and Helllwell [l977] to explain whistler 
precursors. We conclude that the longitudinal resonance Is a likely 
candidate to drive a process In which a whistler wave perturbs the 
particles along a field line through longitudinal resonant bunching. 

This bunching has the effect of creating an enhancement, near the 
equator, of particle flux in a particular parallel velocity range. The 
enhancement Is of sufflcl'vat amplitude and duration to permit a 
gyroresonance Interaction with wave activity such as power line 
harmonics. We fxnd that the longitudinal resonance Is not at first look 
a likely process for creating coherence In Cerenkov process of hiss 
generation, but that features of the longitudinal resonance may merit 
further study In this direction. Also presented was an analysis of the 
limiting electron flux for the bunching, l.e. we estimated the electron 
densltj' at which the Internal fields of the bunch may become large 
enough to affect the bunching process. 
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B. SUGGESTIONS FOR FUTURE WORK 

In our presentation we have shown the results of computer 
simulation of the nonlinear longitudinal resonance Interactions with 
constant frequency whistler mode waves In the magnetosphere. This work 
could be further extended as described below: 

I) We have Indicated In Chapter V that both the wave amplitude 
(Ell) and the wave normal angle are treated as though they are constant 
quantities. It was said that this approximation will be valid as long 
as the Interaction region Is small, but there may be cases where it Is 
necessary to Include effects due to the variation of those quantities. 
The wave amplitude can be computed as a function of position using a 
standard WKB approach, while the wave normal angle variations can be 
calculated using a ray tracing analysis. Those additional computations 
could either be done separately and entered as data, or they could be 
added to the existing code. 

II) Another extension of the present work could deal with CW 
pulse signals propagating along the field line. In this case It should 
be realized that the wave group and parallel velocities have In general 
different values (except for the Gendrln mode^ which poses additional 
problems. It can be easily visualized that an electron trapped In the 
wave potential well, l.e. an electron whose parallel velocity Is very 
close to the wave phase velocity, has to slide either backward or 
forward through the wave packet when the group velocity Is either 
smaller or larger than the phase velocity, respectively; for f < f^/2 
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the whistler mode group velocity always exceeds the phase velocity. In 
the case of a CW pulse signal It would also be possible for electrons to 
enter the wave packet from both ends, depending on the ratio of their 
parallel velocities, and the group velocity of the pulse. 

From the above discussion It Is obvious that this problem would 
require significant changes In the present program, but could also 
reveal some additional features of the longitudinal resonance. 

Ill) Another extension of the work presented here would be to 
Investigate the longitudinal Interaction for the case of variable 
frequency pulse signals. In this case the calculations would have to 
take Into the account the fact that different frequencies of the signal 
Interact with different electrons, and also at different locations along 
the field line. It should be feasible to Investigate the behavior of 
whistlers Interacting with energetic electrons by approximating the 
whistlers with an appropriate number of segments of linearly changing 
frequency, as was done In the discussion of the precursor. 

Iv) It was noted earlier that the wave amplitude may be 
significantly changed due to the Interaction, especially In a full 
distribution simulation.. Although In our particular case In Chapter VI 
It was found that the total energy exchange Is small. It will change for 
other distribution functions. For example, If we assumed a v"^ Instead 
of a v”*^ dependence, there would be many fewer electrons at higher 
parallel velocities, as the weighting function would change from v to 
v”** (see Fig. 6.4). In this case there would be more energy transferred 
from the wave to the trapped electrons compared to the energy 
transferred from the untrapped electrons to the wave, and the final 


result would be wave attenuation. Thus In cases like this It may become 
necessary to Include an energy feedback term that accounts for the 
amplitude changes. However, for a single particle simulation this 
feedback effect Is very small and 


can be omitted 


APPENDIX A; USEFUL IDENTITIES 


Below Is the list of Identities used In the derivation of time averaged 

equations of motion, as well as the derivation of an approximation for 

the <qv 3 > term for small 9. 

y X 

cos(y ” |i8in(j)) ■ cos Ycos(nsln<|)) + slnY slrCnsln^)) 

sin(Y - nsln(^) ■ slnY cos(nsln(j)) - cosy sln(nsln(|)) 

cos(nsin4i) ■ Jo(n) + 2 JaCn) cos(2<)») + 2 J4(n) cos(4<()) + ... 

sln(nsin4>) ■ 2 Ji(n) sin(<p) + 2 JsCn) sln(3<j)) + 2 Js(n) sin (5<j>) + ... 
27T 

I cos(Y - nsln4i) d<P - Jo(n) cosY 

0 

1 slni^ cos(y “ nslncj)) d(j> ■ Ji(n) slnY 
o 

r 

I cos4»cos(y ” nsin(<)y d<|> ■ 0 

S 2n 

sin(Y- r(Sln<<>) d<|> ■ Jo(H) slnY 
o 


•2ir 


f 


sln<j> sin(Y - nsin4>) d<(» 


cos<J> sln(Y - nsln<J>) d<j> 


Jx(n) COSY 


0 
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The <qVy2jj^ term (Eq, 2.62) is given as 


k cos 9 


S > - -qE«Ji(n) slny Pz 

7 * 0) 


For small 9 sin 9=9, cos9 = 1, and tan9=9 so that n - — tane tana. 

ujh 

as already found in Section XI. C. Furthermore, we note that 


Vj, k cos 9 


nu> 


U) 


u> 


cos9 


Vx 

V 


and that near the resonance Vp„ = v,, so that 


k cos9 V* 


U) 


• tana. 


'pit 


Therefore, <qv 3 > - -qE„slnYtanap Ji(n)» and for small 9 

y y z 


n _ u)0 tana 


Also, for small 9, p_ls given as 

z 


Pz 


, 1 — U)/o) 

1 H 1 

w/wu 1 + aj/(u„ 9 

H 


Substituting for J; (n) and p 

Z 


in the expression for <qv S > 

y X 


9 


result is 


the final 


<qv 3 > 
^ y X 


-qEiisinYtan^a 


1 - W/Wjj 

2(1 + cj/ujjj) 
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APPENDIX B: PROGRAM LISTING 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
46 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 
67 


C ANGLES ARE IN RADIANS EXCEPT IN INPUT AND OUTPUT 

DIMENSION Z( 3000 > , PH 1 ( 3000 ) . BZC 3000 ) , WN< 3000 ) . VP ( 3000 ) 

DIMENSION FDAT( 20) ,EN0AT( 10) ,ALPDAT( 10 > , KTEMP ( 40 ) 

DIMENSION BESEL( 2) ,ETA( 3000) ,BMULT( 3000) 

DIMENSION RKDZ( 3000) .RKDZK 3000) ,CTHG( 3000) 

DIMENSION AMPLOW( 3000 ) , AMPLH I < 3000 ) 

COMMON OVPA,EQALD.ALGRO,VPA,FVPA( 400 ) , SO I ST. ALEQ , A , SVPA , FD I ST( 18 

1 0,400).EQAL.FPOIST( 1 80 ) , P I , EM . EL . RPH I . VPE . E . EV , KMAX , VMI N . VPMAX . 

2 ALMIN.ALMAX,ALDC( 1 2 ) . R , RO . VPAEQ, EPA . EV0C( 12 ) . IG , EPAG< 3000) 
COMMON/BLOCKI/ KFDISTC 1 80 . 400 ) , IFD I ST( 180.20) 

COMMON/BLOCK2/ SFDISTC 1 80 ) , I IAS , I lAF . NVG . ALFALO . ALFAH I 

1 .ALFA(35).JLO.JHI 

COMMON/BLOCK3/ TC( 400.12) .CARGUC 400 . 1 2 ) . VPHAC 400 . 12 ) 

1 ,VPARA( 400, 12 ) .ENER( 850. 12 ) .PBCARGUC 505 . 12 ) .P3VPH( 505 . 12 ) , 

2 PBVPA< 505,12) ,TMIN.TMAX,TR{ 12),TTRACE( 12 ) . INDEX! 12 ) , 

3 MLO,MHI ,MSTEP ,TEN<850) ,TPB< 505 ) ,DI STAN ( 850) ,DISTAN1< 505 ) 

C 

C Z-ARC LENGTH, PHI-INVARIANT LATITUDE, 

C BZ-< 1/B)*DB/DZ, WN-WAVE NUMBER K, VP- PHASE VELOCITY 

C 

C IREAO IN ALL NECESSARY DATA 

C 

C CHOOSE GENORIN MODE OR NOT < IG-1 OR IG-0) 

READ(5,350) IG 

350 FORMAT! 12) 

READ!5,3S0) ICONT99 
READ! 5, 350) ICONT88 

C 

C COLLISIONLESS MODEL OR DIFUSSION MODEL ! ICLM-1 OR ICLM-0) 

C GET MODEL PARAMETERS TE , X 1 0 , X IH , X I HE , ENEQ 

C 

READ!5,351 ) ICLM, POWER, TEMP ,XIO,XIH ,XIHE , ENEQ 

351 FORMAT! 12, 6F 10. 5) 

C 

C 

C 

C 

C THE WAVE AMPLITUDE IS DIFFERENT IN THE CASE OF GENDRIN MODE 

C THAN IT IS IN NON-GENDRIN CASE. GENDRIN MODE WAVE INTENSITY IS 

C BW V/HILE NON-GENORIN MODE WAVE INTENSITY IS LABELED EPA. 

C THE PROPER SETTING OF WAVE INTENSITIES IS DONE IN FOLLOWING WAY : 

C 1 ) GENDRIN MODE 

C WAVE INTENSITY IS 8W-CONST»2** IBW WHERE I3WLO< IBW< I BWH I . 

C IBWLO AND IBWHI ARE READ FROM INPUT CARO DECK. AT THE SAME TIME 

C EPA IS NOT USED WHICH IS ACCOMPLISHED SETTING IELO-IEHI-1 

C 

C 2) NON-GENDRIN MODE 

C WAVE INTENSITY IS EPA-CONST*2**IE WHERE IELO< IE< lEHI . 

C AT THE SAME TIME GENORIN MODE IS SUPPRESEO USING IBWLO-IBWHI-1 

C 

READ! 5, 352) IBWLO. IBWHI , lELO, lEHI 

352 FORMAT! 412) 

C 

C FREQUENCY ITERATION 

C ENTER THE NUMBER OF DIFFERENT WAVE FREQUENCIES AND THEY WILL BE READ 

C FROM INPUT CARO DECK 

C 

REAO!S,350) INFREQ 
00 1013 ICNT-1 . INFREQ 
REA0!5.3Q3) FOAT! ICNT) 

1013 CONTINUE 

353 FORMAT! F 10.5) 

C 

C READ L VALUE AND ANGLE BETV/EEN K&30 ! THETA) 

C 

READ! 5. 354) EL. THETA 

354 FORMAT! 2F 10.5) 
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68 C 

69 C DEFUtE DIRECTION OF PROPAGATION 

73 C IWO-I — > POSITIVE DIRECTION 

71 C IW0«-l --> NEGATIVE DIRECTION 

72 C 

73 READ(S.3S0) IWD 

74 C 

75 C PARAMETERS ALONG FIELD LIME PRINTED IF IC0NT1>1 

76 C 

77 READ(S,350) ICONTl 

78 C 

79 C FULL DISTRIBUTION USED IF IFULL-1. ADIABATIC APPROXIMATION USED 

80 C BEYOND RESONANCE POINT IF IA0IA«1, DIFFUSION COEFFICIENTS 

81 C COMPUTED IF IOIFF-1 

82 C 

83 C PROGRAM CAN TRACE EITHER A SINGLE PARTICLE OR GIVEN DISTRIBUTION 

84 C GIVEN BY THE DISTRIBUTION FUNCTION FDIST( VPARALEL .ALFAEQ ) . 

85 C 

86 C 1) SINGLE PARTICLES TRACING 

87 C TO DO SINGLE PARTICLE TRACING IT IS NECCESSARY TO SPECIFY ITS 

88 C PARALLEL VELOCITY AND EQUATORIAL PITCH ANGLE. 

89 C THIS IS DONE DEFINING TWO PARAMETERS: IV (LOOP 206) 

90 C AND lA (LOOP 204). 

91 C GIVEN RANGE IVtlVS.IVFI PARTICLE VELOCITY IS GIVEN AS 

92 C VPAI-VMIN*( !♦< IV-1 )/10)*l .05 AND PITCH ANGLE IS READ FROM 

93 C INPUT CARO DECK USING lA AS POINTER WITH RANGE CIAS.IAFl. 

94 C 

95 C 2) FULL DISTRIBUTION TRACING 

96 C IN THE CASE OF FULL DISTRIBUTION ALL DATA CONCERNIG SINGLE PARTICLE 

97 C WILL BE NEGLECTED. THE INITIAL DISTRIBUTION IS GIVEN BY THE NUMBER 

93 C OF BINS IN VELOCITY AND PITCH ANGLE RANGE. 

99 C NUMBER OF VELOCITY BINS IS READ AS INPUT DATA ( NVG > SAME AS PITCH AN 

100 C RANGE CIIAS.IIAF]. 

101 C 

102 REA0<5.355) lAOIA, IFULL . lOIFF 

103 355 FORMAT! 312) 

104 C 

105 IF( IFULL. EQ.l ) GO TO 1015 

106 REA0<5,352)IVS.IVF.IAS.1AF 

107 DO 1014 ICNTl-IAS, lAF 

108 REA0(5.3S3) ALPOAT< ICNTl ) 

109 1014 CONTINUE 

110 GO TO 1016 

111 1015 REA0(S,341) NVG . I IAS , I lAF , VRANGE , VIN ITL 

112 341 FORMAT(3I2,2F10.5) 

113 1016 CONTINUE 

114 IF< IFULL .EQ. 1 ) IVS-1 

115 IF( IFULL .EQ. 1 ) IVF-1 

116 1F( IFULL. EQ.l ) IAS-1 

117 1F( IFULL .EQ. 1 ) IAF-1 

118 C 

119 C PRINT PHASE ANGLE YES-1. NO-0 

120 REA0<5.363) ICOMT2 .MLO.MHl .MSTEP .TMIN ,TMAX 

121 363 FORMAT(4I2,2F10.5) 

122 C 

123 C READ THE STARTING LATITUDE WHERE TRACING SHOULD BEGIN 

124 C 

125 READ(5,353) SRPHIO 

126 C 

127 C READ WAVE AMPLITUDE INFORMATION 

128 READ<5,350) IGROW 

129 REA0(5.3S8> XPHIOO .XLEN , XAMPL 

130 3E8 FORMAT! 3F10. 6 ) 

131 READ(5.3S0) IC0NT5 

132 READ(5,350) 1C0NT25 

133 READ(S.3S3) XMAX 

134 READ(5.353> VOELTA 
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135 

136 

137 
136 
139 
U6T 

141 

142 

143 

144 

145 
14S 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 
187 
1G8 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 


C 

C 

C ITERATE FOR DW IN GENDRIN MODE. IF GM IS NOT USED SET IBW-I.l 

WRITE( S.70J7 > IG, ICLM.TEflP ,;nO,XIH ,XIHE .ENEO 
7007 FORMAT(2I3.5F10.5) 

WRITE( 6. 7008 ) IBWLO. IBV/HI , lELO, lEHI , INFREQ,FOAT( I ) 

7009 FORflAT(SI3,Fl0.5) 

WRITE (6,7009) EL, THETA. IWD, ICONTl .lAOIA.IFULL.IOIFF 
7.009 FORMAT(2F10.5.5I3> 

WRITE! 6,7010) IVS, IVF, IAS. lAF, ICONT2,SRPHIO,ALPOAT( 1 ) 

7010 FORMAT! 5 1 3. 2F 10.5 ) 

00 210 IBW-IBWLO. IBWHI 
BW»3.75E-12*2’*«IBW 
DO 208 IEF«1 .INFREQ 
F-FDAT< lEF ) 

C 



C DEFINE ALL NEEDED CONSTANTS 

C 

E-1 .6021E-19 
C-2.9978E8 
PI-3.1416 
R0-6.37E6 

PHI0-ATAN(SQRT(EL-1 . )) 

A-3. 1416927/180. 

EM-9. 1066E-31 
DZ-1 .E4 
R1-7.37E6 
CTII-COS(THETA*A) 

STh-SIN<THETA*A) 

OM-E.^PI'P 
BOLTZ-1 .3805E-16 
EMI-9. 1066E-28*1837. 

G1-980.67*RO’*RO/R1/R1 
OMS-(PI/12./3600. )**2 

C 

C 

C TEST PROGRAM FOR FULL DISTRIBUTION 

IF( ICONT88.EQ.0) GOTO 713 
WRITE<6,3958) 

3958 FORMAT(////’TEST BESSEL FUNCTION COMPUTATIONS'//) 

ARG-0. 

CALL BESJR(ARG.1,BESEL.IER) 

WRITE! 6, 3956) ARG,BESEL( 1 ),BESEL(2) 

3956 FORMAT! 3F 12. 4) 

ARG-1 . 

CALL BES3R!ARG, 1 .BESEL. lER) 

WRITE! 6.3956) ARG.BESEL! 1 ),8ESEL!2) 

713 CONTINUE 
C 

c 

C DENSITY MODEL DATA ARE READ FROM INPUT CARD DECK 

C COMPUTE PF!PLASMA FREQUENCY), FH !GYROFREQ.) AND R IND! REFRACTIVE 

C INDEX) ALONG GIVEN FIELD LINE USING QL APPROXIMATION. 

C ALS( COMPUTE WN!WAVE NUMBER IN MAG. FIELD DIRECTION) AND VP!PHASE 

C VFi.^' rv IN MAG. FIELD DIRECTION). 

C '.!N) AND PHI!N) GIVE POSITION ALONG THE LINE. 

C WN AND VP ARE DIFFERENT FOR GENORIN AND NON-GENDRIN MOOES. 

C 

C 

HH-B0LTZ«TEMP/EMI/G1*1 .E-2 
C SCALE HEIGHTS ARE CONVERTED TO METERS 

HHE-HH/4. 

HO-HH/16. 

GPHEQ-Rl-Rl*Rl/R0/EL-0MS/2./Gl/R0/EL*! ! RO-EL )**3-Rl**3) 
ENFAC»XIH*EXP! -GPHEQ/HH )+XIHE*EXP! -GPHEQ/HHE I+XIO^EXP! -GPHEO/HO ) 
EHFAC-ENEQ/SQRT! ENFAC ) 

N-1 


I 
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203 



Z(N )-0. 

204 



PHKN )-0, 

205 



R-RO“EL 

205 



COEL-1 . 

207 



BZ< N )-0. 

208 



FP-SQ«T( 80.6"ENEQ"1 .E6> 

209 



FH-8. ,'36E5/EL*"3 

210 



RIND-FP/SQRT{ F«{ FH*CTH-F > > 

211 

C 


VP(N) IS PHASE VELOCITY IN MAG FIELD DIRECTION 

212 

C 


IF 6ENDRIN MODE USED THAN NEXT LINES EXECUTED, OTHERWISE 

213 

C 


GO TO 11 

214 

C 



215 



IF< IG .NE. 1 ) GO TO 11 

216 



VP(N)-C/2.*FH/FP 

217 



CTHG<N)-2.*F/FH 

218 



WN<N)-2.*PI*F/VP(N)/CTHG(N) 

219 



EPAG(N)-C*BW*F/FP«SQRT» 1 . -4 . »F*F/FH/FH ) 

220 



GO TO 12 

221 


11 

VP<N)-C/RIND/CTH 

222 



WN<N)-RIND/C»2.*PI*F 

223 


12 

RKDZ(N)-0. 

224 

C 



225 

C 


NEXT LOOP (LABEL 10> COMPUTES ALL MEDIUM PARAMETERS ALONG GIVEN 

226 

c 


FIELD LINE 

227 

c 



228 


10 

N-N>1 

229 



Z(N)-Z(N-1 >*DZ 

230 



PHI<N)-PHI(N-1 )-*.0Z*C0EL/R 

231 



CPHI-COS(PHI(N>) 

232 



SPHI«SIN(PHI(N) ) 

233 



R-PX*EL*CPHI*»2 

234 



SRF-SORT< 1 .+3*SPHI**2> 

235 



CDEL-CPHI/SRF 

236 



SDEL-2.*SPH1/SRF 

237 



B2(N)-3./R*( SPHI*CPHI»COEL/SRF/SRF+SDEL ) 

238 

c 


BZ IS DELTA B OVER DELTA Z DIVIDED BY B 

239 



Z( N)-RO/2./SQRT(3. ) /C03( PH 10 > **2*( ALOG ( SQRT( 3 . >*SPHI^SRF ) 

240 


1 

♦ S0RT<3. )*SPHI*5RF ) 

241 



GPH-Rl-Rl*Rl/R-OMS/2./Gl/RO/EL*(R**3-Rl**3> 

242 



EN»XIH*EXP( -CPH/HHl+XIHE^EXPI -GPH/HME ) +X 10-EXP ( -GPH/HO > 

243 



EN-SQRT(EN)-ENFAC 

244 



IF( ICLM.EQ.l ) EN-ENEQ*(RO*EL/R)**POWER 

245 



FP-SQRT< 80.6-EN-l .E6 ) 

246 



FH-8.736E5-(RO/R)*-3-SRF 

247 



RINO-FP/SORT( F-<FH-CTH-F ) > 

248 



FACTl-l-( FP/F )--2 

249 



FACT2«l-FP-*2/{F*-2-FH«-2) 

250 



FACT3«(FH/F )*FP— 2/( F--2-FH — 2 ) 

251 



IF( IG.NE . 1 ) GO TO 14 

252 



VP< N >»C/2.-FH/FP 

253 



CTH6(N)-2.*F/FH 

254 



V/M(N)-2.*PI*F/VP<N>/CTHG(fl) 

255 



RIN02-<FP/F><*-2 

256 



STHG2-1-CTHG<N>— 2 

257 



STHG»SQRT(STHG2) 

238 



EPAG(N l-C-avZ-F/KP-SORTC 1 . -4 . -F -F/FH/FH ) 

259 



RKD2(N)»(V/N( N)*WN( N-1 ) ) /2 . -OZ-CTKGl N ) + R:<DZ( N- 1 > 

260 



GO TO 15 

261 


14 

WN( N )»RIN0/C-2."PI-F 

262 



VP< N )-C/RIMO/CTH 

263 



RKDZ(N)-<V/N(N)+WN<N-1 ) ) /2 . -OZ-CTH*R!<DZ( N- 1 ) 

264 



CTHG( N >-CTH 

265 



RIND2«RIND"*2 

266 



STHG2-STH--2 

267 



STHG-STH 
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263 

15 

IF ( VP( ri ) .GT.VPI N-1 > > VPMAX-VP(, 

269 


BMULTl N >-FACT3/( RIN02-FACT2 )*( Rl, j2*STHG2-FACT 1 )/RIN02 

27,7 

1 

/STHG/CTHG< N ) 

271 


ETA(;i)-WN<N)-3THG/FH 

272 


MMAX-N 

273 


IF (R.GT.RO) GO TO 10 

274 

C 

ALL PARAMETERS COMPUTED 

275 


N-0 

276 


RPH!"SRPHIO*A 

277 

46 

M"N^1 

278 


IF(ABS(RPHn.GT.PHI(NM GOTO 46 

279 


INDMAX-N 

280 


RKOZL(N)-0. 

281 

47 

N‘N-1 

282 


RKI32L(N)»<WN<N+1 >+VN< N ) >/2 . *02*CTH+RKDZL( N+1 ) 

283 


IF(N,GT. 1 ) GOTO 47 

284 


00 48 N-l.NMAX 

285 

286 

48 

RKDZ(N)-RKDZ(N)^RKOZL( 1 > 



287 

c 


288 

c 

TO PRINT PARAMETERS ALONG FIELD LINE ICONTl-1 

289 

c 


290 


IF< ICONTl .NE. 1 ) GO TO 6000 

291 


I«0 

292 

6002 

H- I "10+1 

293 


IF(N.GT.NMAX) GO TO 6000 

294 


PHID-PHI(N>/A 

295 


WRITE (6. 6001 > PHIO.Z(N) ,EPAG(N},VP(N>.CTHG(N),UN(N) 

296 

297 

298 

299 

6001 

FORMAT! F 10. 2, SE 12. 4) 
1"I*1 

GO TO 6002 

u 


300 

c 


301 

6000 

CONTINUE 

302 

C 


303 

C 

THIS CODE WILL BE EXECUTED IF VARIABLE AMPLITUDE WAVE 

304 

C 

IS USED PROGRAM 

305 


IF( IGROW.NE. 1 ) GO TO 8061 

306 


XPHIO-XPHIOO*A 

307 


XSTART»R0/2./SQRT( 3 . )/COS( PH 10 ) • »2*( ALOG( SORT! 3 . >*SIN 

308 

1 

(XPHIO> + SQRT( 1 .+3.*SIN< XPHIO >*<»2 > > + SQRT{ 3 . )“S I N( XPH 10 ) • 

309 

2 

SORT! 1 .f3.*SIN(XPHIO)**2 ) ) 

310 


XEND-XSTART+XLEN*1000. 

311 


DO 8032 I >1,3000 

312 


AMP LOW! I )»0. 

313 

8032 

CONTINUE 

314 


DO 8033 1-1,3000 

315 


AMPLHI! I >-0. 

316 


IF! ! PHI! I ).GT. 0.1 221 7) .AND.! PH I! I ) . LT .0 . 1 7453 > ) AMPLH I ! I )-45 . E -6 

317 

8033 

CONTINUE 

318 

8061 

CONTINUE 

319 

C 


320 

C 

AMPLITUDE DATA STORED 

321 

C 


322 

C 
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323 

C 


324 

C 

INITIALIZE FINAL DISTRIBUTION FUCTION TO 0 IF FULL DISTRIBUTION 

325 

C 

IS USED IN PROGRAM. 

326 

C 

THE INITIAL DISTRIBUTION IS SET UP ACCORIDINGLV TO MVG FOR VELC 

327 

C 

BIN AND HAS AND I lAF FOR PITCH ANGLE 3 1 ft . 

328 

C 

THE FINAL DISTRIBUTION BINS ARE COMPUTED FOR VELOCITY TO GIVcI 

329 

C 

THE BEST RESOLUTION AND FIXED FOR PITCH ANGLE (0.5 DEGREE IN 

338T 

C 

0-90 RANGE ) 

331 

C 


332 


IF { IFULL .EQ.0) (20 TO 43 

333 

C 

J FOR ALPHA GOES FROM 1-180 

334 

C 

NVG IS NUMBER OF GRIDS IN VPARALLEL IN INITIAL DIST FUNCT 

335 


DVPA»VP( 1 )«VRANGE/{ NVG+1 ) 

336 


K-1 

337 


FVPA< 1 >«0.25*VP< 1 ) 

338 

40 

K-K*l 

339 


FVPA(K>-FVPA<K-1 )+DVPA*10 

340 


IF<FVPA{K).LT.(VP( 1 )*3.24>) GOTO 40 

341 


KMAX-K 

342 


DO 42 K-1 .KMAX 

343 


CO 41 J-1.180 

344 


IF(K.LT.21> IFOISTI J,K)-0 

345 


KFOIST( J,K)-0 

346 

41 

FDISTCO .K >-0. 

347 

42 

CONTINUE 

348 

43 

CONTINUE 

349 

C 


350 

C 


351 

C 

'PARTICLE TRACING STARTS 

352 

C 

ITERATE FOR WAVE INTENSITY 

353 

C 

FOR GENDRIN MODE WAVE INTENSITY IS SPECIFIED BY 

354 

C 

MAGNETIC COMPONENT NEAR BEGINNING OF PROGRAM. 

355 

C 

NEXT DO LOOP SHOULD HAVE ONLY ONE LOOP ( IBWLO-IBWHI-1 > 

356 


DO 207 lE-IELO.IEHI 

357 


EPA-4S.E-6 

358 


IF < ICONT25.EQ.0) GOTO 4080 

359 


EPA-1 .E-6*XMAX 

330 

4080 

IF(EPA.EQ.0) IMAX-1 

361 


IF (EPA.NE.0) IMAX-12 

362 

C 

FOR GENDRIN MODE EPA IS REPLACED BY EPAG< 1 ) FOR OUTPUT PRINTING 

363 


IF < IG.EQ. 1 > EPA-EPAG( 1 ) 

364 


VFMIN-1 .E16 

363 

C 


366 

C 

INITIALIZATION OF PLOTTING DATA ARRAYS 

367 

C 


363 


if: ICONT2.EO.0) GOTO 1721 

369 


00 1716 1-1,12 

370 

1716 

TR( I >-100. 

371 


00 1717 1-1,850 

372 


00 1761 J-1,12 

373 


ENER< I . J )--l . 

374 

1761 

CONTINUE 

375 

1717 

CONTINUE 

376 


00 1713 1-1,12 

377 


00 1719 J-1,505 

373 


ir<J.GT.400> GOTO 1720 

379 


TC(J,I )-l.E36 

380 


CAPvGU( J, I )-l .E36 

381 


VPHA( J, I )-l .E36 

382 


VPARA(0 , I )-l .E36 

383 

1720 

PBCARGUl J , I >-: .E36 

334 


PBVPH(J. I )-l .E35 

385 


PBVPA( J. I >-l .E36 

336 

1719 

CONTINUE 

387 

1718 

CONTINUE 


. ‘r 


388 

289 

39.9 

391 

392 

393 

394 

395 

396 

397 

398 

399 

400 

401 

402 

403 

404 

405 

406 

407 

408 

409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 

433 

434 

435 

436 

437 
433 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 

451 

452 

453 

454 

455 
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DO 1762 1-1,350 
IF( I .GT.505 ) GOTO 1733 
TPB( I )-l .E36 
D I STANK I )«1 .E33 
1763 TENC I >-l .E36 

DISTAN( I )“1 .E36 
1762 CONTINUE 
1721 CONTINUE 
VFMAX-0. 

JCOUNT-0 

EOTOT-0. 

tFTOT-0. 

ALFALO-1 .E10 
ALFA!! 1-0. 

IF( IFULL.EO.0) VINITL-l. 

VMI"l-UINITL*VP( 1 > 

C ITERATE FOR PARTICLF VELOCITY 

C IVS AND IVF ARE VELOCITY RANGE DATA FOR SINGLE PARTICLE TRACING 

IF (IFULL.EQ.l) IVF-IVS 
DO 206 IV-IVS, IVF 
VPAI-VMIN»< 1 . 102+IV-0.001 > 

IF ( ICONT25.EO.0) GOTO 4081 
VPAI»VP( 1 l^VDELTA 
4081 IIVS-1 

IF( IFULL.EO.0) NVG-0 
IIVF-NVG+1 

IF <IFULL.FQ.0) IIVF-IIVS 
DO 205 riV-IIVS, IIVF 
VPAI I-VMIN + DVPA"( I IV-1 ) 

IF( { IIV.EQ. nVS) .AN0.( IFULL.EQ.l)) VSTART-VPAII 
IF< < I IV.EQ. IIVF ) .AN0.( IFULL.EQ.l)) VEND-VPAII 
IF < IFULL.EO. 1 ) SVPA-VPAI I 
IF (IFULL.EO.0) SVPA-VPAI 
C ITERATE FOR EQUATORIAL PITCH ANGLE 

C lAF AND IAS ARE PITCH ANGLE RANGE DATA FOR SINGLE PARTICLE TRACING 

IF ( IFULL.EO. 1 ) lAF-IAS 
DO 204 lA-IAS.IAF 
ALEQI-ALPDAT< lA) 

C HAS AND IIAF ARE PITCH ANGLE RANGE FOR FULL DISTRIBUTION 

IF< IFULL .EQ.0) IIAF-1 
IF( IFULL.EO. 0) IIAS-1 
IF < IFULL.EO.0) HAF-IIAS 
ALMIN-5.25+0.5*IIAS 
ALMAX-5.25+0.”*‘HAF 
DO 203 IIA-IIAS.IIAF 
ALEQII-5. 25-^0. 5»IIA 
IF (IFULL.EQ.l) ALEQ-ALEOn 
IF (IFULL.EO.0) ALEQ^ALEQl 
If (IFULL.EO.0) WRITE (6,998) ALEQ 
C ALEQ IS IN DE6,REES 

998 FORMATdHl,' F.Q P ITCH ANGLE- ‘ ,F7 , 3/ ) 

C ITERATE FOR BETA 

DO 2.02 I -1,1 MAX 
6ETAD-30.-I-30. 

8ETA-CETAD*A 

C STARTING LATITUDE IS INPUT DATA 

RPHI-SP.PHID-A 
SPHI-SIN(AfS(RPHI )) 

CPHI-COS(AflS( RPHI ) ) 

3RF-SQRT( 1.+3.«SPHI**'2) 

S-RO/2./SQRT(3. )/COS( PH 10 )**2»( ALOG( SORT! 3 . ) 

1 -SPHI+SRF )>SQRT( 3. )*SPHI*SRF ) 

IF ( PvPHI .LT.0;- S-0.-S 
TANS«TAN(ALEU»A >*•2 
FHRaT®SURT( 1 .+3 .-SPHH-2 )/CPHI**6 
VPA-SVPA«*SQRT( 1 . +TANS-FHRAT’*TANS ) 

SVPE»SVPA>»TAN( ALEQ*A) 

VPE-SVPE-SGRT(FHRAT) 
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456 

457 
453 

459 

460 

461 

462 

463 

464 
463 

466 

467 

468 
4c9 

470 

471 

472 

473 

474 

475 

476 

477 

478 

479 

480 

481 

482 

483 

484 

485 
436 

487 

488 

489 

490 

491 

492 

493 

494 

495 
436 

497 

498 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 

513 

514 
sl5 

516 

517 

518 

519 


EO-EM/2 . •< Vr:'*VPE+V?A*VPA > 

EVO-EO/E 

IF( IFULL .EQ. 1 ) GOTO 135 
IF( I .ME . 1 ) GOTO 135 
IF< ICONT25.EQ.I ) RATIO-VDELTA 
IF( ICONT25.EQ.0> RATIO-VPAI/VMIN 
135 COMTinUE 

IF( < I .EQ. 1 ).AN0. (IFULL.EQ.0) > WRITE ( 6 . 705 1 > SRPHID 
7051 FORMAT! • TRACING STARTS AT '.Fo.2,' DEGREES LATITUDE* > 
IFIIGROW.EQ. 1 ) EPA-XAMPL 

IF(< IFULL.EQ.0>.AND.< I.EQ.l > ) WRITE! 6,999 ) EL , ENEQ, F , SVPA , EVO, EPA 
1 ,VPA, VP! 1 ) .RATIO 

999 FORMAT !* EL - ' ,F5 . 2 , 3X . ‘ EQ DEN- ' ,F6 . 1 , ' CM-3 ' , 3X . ‘ FREQ- * , -3PF6 . 3 , 

1 'KHZ’,3X,‘EQ PAR VEL - ' ,0PE 10 • 3 , * M/SEC 3X . ‘ INIT ENERGV- ' , E 1 2 . 6 , 

2 ■ EV’,3X. EPA-',E10.4,‘V/I<r/' VPA- ’ , E 1 1 . 4 , • M/S ’ , 3X , 

3 * EQ PHASE VEL»' ,E11 .4, *M/S' ,3X, 'RATIO! VPAR/VPHASE >-• ,F7.5 > 
IROONE-0 

lROON-0 

ITl-0 

IT2-0 

IC-0 

INCi>0 

IC2-0 

IMOONE-0 

IMIRR-0 

T-0. 

DT-0.0001 

IT-0 

N-1 

100 N-N+1 

IF (ABS!S).6T.Z!N>) GO TO 100 

NU-N 

ML -N-1 

IF! I.EQ.l) WRITE!6,49) INDMAX , NL . NU 
49 F0RMAT!//3I5//) 

VPHASE-IWD*!VP!NL )♦! VP! NU )-VP! NL ) )*! ABS! S )-Z! NL ) )/! Z! NU )-Z! NL ) ) ) 
IF !VPA.GF.!VPHASE»IWD)) ITEST-1 
IF !VPA.LT.!VPHASE*IWD) ) ITEST--1 
1 10 BZF-< BZ! NU )-0Z< NL ) )*! ADS! S )-Z! >.L ) )/!Z! NU )-Z! NL > )+BZ! NL ) 
IFIS.LT.0.) RKDZ1-RKDZL!NL ) 

IFIS.GE.0.) RKDZl-RKDZ! NL ) 

IFIS.LT.0.) RKDZ2-RKDZL!NU) 

IF!S.GE.0.) RKDZ2-RKDZ! NU ) 

RKF-IWD*! RKDZ1+! RKDZ2-RKD21 )*! ABS! S )-2! NL ) )/! Z! NU )-2! NL ) ) 

1 ) 

IF !S.LT.0. ) BZF — 1*BZF 
CARG-OM*T-RKF+BETA 

IF!! IGP.OW.EQ.l ) .AND . ! S . LE .0. )) EPA-AMPLOW! NU ) 

IF!! IGROW.EO.l ) . AND . ! S . GE .0. )) EPA-AMPLH I ! NU ) 

IF! ICONT99.EQ.0) GOTO 3703 
COSINE-CTH 

IF! IG.EQ. 1 ) COSINE-!CTHG!NU )+CTHG! NL ) )/2. 

TE RM 1 -VP E • ! WN ! NU ) +WN ! NL ) ) / 2 . *COS I ME /F / 2 . /P I 
ARG-!ETA!NU) + ETA!NL ) )/2.»VPE 
CALL BESJR!ARG.1,8ESEL,IER) 

TERM3-BESEL! 1 )•! l-TERMl"! BMULT! NU ) ♦3MULT! NL ) )/2 . *BESEL ! 2 ) 

1 /BESEL!!)) 

GOTO 3709 

3708 TERM3-1. 

3709 CONTINUE 

IF! IG .NE. 1 ) GO TO 5000 

EPAF“EPAG!NL)+!EPAG!MU)-EPAG!NL ) )*! ABS! S )-Z! NL ) )/!Z! NU )-Z! NL ) ) 
VPAT-VPA-VPE«'«2/2.’*G2F‘’OT-E/EM’»£PAF“TERM3’*COS!CARG)-DT 
GO TO 5001 


4 



520 

52' 

522 

523 

524 
523 
523 

527 

528 

529 

530 

531 

532 

533 
G34 

535 

536 

537 

538 

539 

540 

541 

542 

543 

544 

545 

546 

547 

548 

549 

550 

551 

552 

553 

554 

555 

556 

557 

558 
553 
EC0 

561 

562 

563 

564 

565 

566 
557 

568 

569 

570 

571 

572 
572 

574 

575 

576 

577 

578 

579 

530 

531 

582 

583 

584 
535 

586 

587 
533 
589 
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5000 V?A7-VPA-VPE*’*2/2.»BZr'*DT-E/EM'*EPA»TE2M3'*C0S( CA?vG )-'DT 

5001 ST*S*(VPAT+VPA)/'2."0T 
NUO-NU 

IF (ABS(ST).LE.ABS(S) ) GO TO 101 
MU«NU-1 

102 NU'i'IU + 1 

IF <ABS(ST'.GT.Z(NU) ) GO TO 102 
M.-NU-1 
GO TO 104 
101 NL-NL+1 

103 NL-NL-1 

IF (ABS(ST>.LT.Z(NL>) GO TO 103 
NU-NL+1 

104 CONTINUE 
BZS-<BZ(NU>-BZ(NL> 

IF(S.LT.0.) RKDZ1-RKDZL( NL ) 

IF<S.GE.0.) P.KDZl-aKDZ( NL ) 

IF<S.LT.0.) RKDZ2»P.K0ZL(NU) 

IF<S,GE.0.) RKDZ2«RKDZ( NU > 

RKS-U/0«(RKDZ1-XRKDZ2-RKDZ1 ) »< ABS( ST ) -Z< NL ) )/( Z( NU )-Z( NL > 

1 ) ) 

IF< ST.LT.0. > BZS--1*BZ5 
CARG-OM*T-0.5"( RKF+RKS )+BETA 
IF< ICONT99.EQ.0) GOTO 3715 
COSINE-CTH 

IF< IG.EQ.l ) COSINE-<CTHG<NU>+CTHG(N’.)>/2. 

TERM1-VPE*< WN( NU )+WN< NL > ) /2 . *COS I NE /F /2 . /P I 
ARG-( ETA( NU l+ETAC NL ) )/2. -*VPE 
CALL BESJR(ARG,1 .BESEL.IEP) 

TERM3-BESEL< 1 )>•( 1 -TERMl “I BMULT( NL ) +BMULT{ NU ) ) /2 . *BE SEL ( 2 ) 

1 /BESELCD) 

GOTO 3716 

3715 TERM3-1. 

3716 CONTINUE 

IF< IG.ME . 1 ) GO TO 5005 

EPAS»EPAG< NL > + < EPAG( NU )-EPAG( NL > )>*( ASS ( ST )-Z( NL ) )/(Z( NU >-Z( NL ) > 
VPAT«VPA-VPF*»2/4 . aZF+BZS )’*DT-E/EM*( EPAF + EPAS ) 

1 /2.»TERH3»COS<CAaG>'*DT 
GO TO 2400 

5005 V?AT-VPA-VPE»*2/4 .*< 8ZS + BZF )«OT-E/EM*EPA»T2RM3'‘COS(CARG >'»0T 

2400 IF( IG.EQ.l) EPATEM»< EPAFh-EPAS )/2. 

IF( IG.NE . 1 ) EPATEM-EPA 

IF< ICONT99.EQ.0) GOTO 3726 

TERM2»VPAT*< WN( NL >+WN<NU) ) /2 . *COS I NE /F /2 . /P I 
TERM4-BESEL( 2 )•< BMULT< NU ) + BMULT< NL > >/2.*< 1-TERM2 ) 
VPE»VPE*VPAT*VPE/4.'»( BZS + BZF > *DT+E/EM*TERM4*COS< CARG ) '*OT 
1 -EPATEM 
GOTO 3727 

3728 VPE-VPE+VPE'»VPAT/4.*(BZSfBZF >*0T 
3727 CONTINUE 

SC-S+( VPA+VPAT)/2.»DT 
C CHECK FOR EQUATOR CROSSING 

IF ( ( SC*S ) .GT.0) GO TO 2401 
CALL EQCOMV 

IF( IFULL.EQ.0) WRITE(6,2402> T . EV , EOALD . ’ RDONE 
2402 FORMAT (' EQUATOR Id NG ' , 3X , ' T- ' . F7 . 4 . 3X . ' E NERGY» ' , E8 . 3 , ' E V . 

1 3X,'EQ PITCH ANGLE-* ,F6.3,3X, 'NO OF RE SONANCES- ' , 1 3 > 

IRDONE-0 

2401 CONTINUE 

C FIND MIRROR POINT 

IF ( IMOOME.EQ. I > GO TO 2500 
IF < ( VPA-VPAT) .LT.0) IMIRR-1 
IF ( IMIRR.NE. 1 ) GO TO 2500 
CALL EQCOMV 

IF (IFULL.EQ.0) WRITE (6,253) H , RPH I D , S , T , EV . EQALD 
253 FORMAT (' MIRROR POI NT * . 3X . ' H» * . E 1 2 . 5 . * KM * , 3X . ' P H I » * , F7 . 3 . 3X . 

1 *S-* .E12.5,3X, 'T»* .F7.4.3X, •ENERGY-' ,E8.3, *EV ,3X, 

2 *EQ PITCH ANGLE-' ,F6.3 ) 

IMDONE-1 

GO TO 312 
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S90 


2500 

VPA-VPAT 

591 



ENGY-EM/2./E*(VPE*VPE*VPA»VPA) 

592 



ERROR-ENGY-EO/E 

593 



AL«ATAN( VPE/YPA ) 

594 



IF( ABSt SC ) .LE .A3S( ST ) > GO TO 105 

595 



NU«iNU-l 

596 


106 

NU-nU+1 

597 



IF < ABS( SC ) .GT.Z( NU ) > GO TO 106 

598 



NL-NU-1 

599 



GO TO 108 

600 


105 

NL'NL^l 

601 


107 

NL-NL-1 

602 



IF <A3S(SC>.LT.7{NL)> 60 TO 107 

603 



NU-NL+1 

604 


108 

CONTINUE 

605 



NUO=NU 

60S 



S-SC 

607 



RPHI-(PH1<NU)-PHI(NL ) >*< ABS( S )-Z( NL > )/( Z( NU >-2( NL ) ) + PHHNL > 

608 



IF < S.LT.0. ) RPHI-0.-RPHI 

609 



RPHID-RPHI/A 

610 



R-RO^EL-COSIRPHI >**2 

611 



H-(R-RO)/1000. 

612 



VPHASE-IWD*(VP<NL) + <VP{NU>-VP(NL)>*(ABS(S)-Z(IIL>)/{Z<NU)-Z{NL,)) 

613 

C 


FIND RESONANCE POINT 

614 



IF < <VPA*IWD).LT.0) GO TO 250 

615 



IF (<(VPA-VPHASE)*ITEST).LE.0) GO TO 251 

616 



GO TO 250 

617 


251 

CONTINUE 

618 



IF( ( IFUU.EQ.0).ANO.( IRDON.EQ.0)) TR( I )>T 

619 



IRD0N-IRD0N4-1 

620 



IF( IFULL .EQ. 1 > GOTO 137 

621 



CARGO-CARG/A 

622 


139 

IF(AES(CARGD).LT.360. > GOTO 138 

623 



IF(CARGD.GT.0. ) CARGO-CARGO-360. 

624 



IF(CARGD.LT.0. ) CARGO-CARGD>360. 

623 



GOTO 139 

626 


138 

IF(CARGO.LT.0. ) CARGO-CARGO-«-360. 

627 


137 

CONTINUE 

628 



IF(( IFULL.EQ.0).ANO.( I ROONE . EQ .0) ) WRITE! 6,252 ) VPMASE . R . RPH 1 0 , S . 

629 


1 

.CARGO 

630 


252 

FORMAT {■ RESONANCE VEL » ' , E 1 2 . 5 . 5X , ' AT R- ' , E 1 2 . 5 . 5X , • PH I ■ ’ , F7 . 3 , 

631 


1 

5X. ’S-' ,E12.5,5X, 'T-' ,F7.4,3X. ‘BETA-’ ,F7.2) 

632 



IROONE-IROONE4-1 

633 



ITEST-0-ITEST 

634 


250 

CONTINUE 

635 



T-T+DT 

636 

C 


THE NEXT CARD .GO TO 300, BYPASSES WRITING OF PHASE ANGLI- 

637 

c 



638 

c 


SAMPLING OF PLOT OATA 

639 

c 



540 



IF((ICONT2.EQ.0).OR.(IFULL.EQ.in GOTO 1732 

64' 

c 



642 

c 


RESONANCE POINT SAMPLING 

643 

c 



644 



IF(T.GT.5.0) GOTO 1732 

645 



IT-IT+l 

646 



IF< IT.LT.20) GOTO 1726 

647 



IF(ABS(<VPA-VPHASE)/VPA).GT.0.10) GOTO 1729 

648 



IF((T-TR( I )>.GT.0.20' GOTO 1729 

649 



CARGO-CARG/A 

650 


1727 

IF <ABS<CARGO ) .LT.360. ) GOTO 1728 

651 



IF<CARGD.GT.0. ) CARGO-CARGO-360. 

652 



IF<CARGO.LT.0. ) CARGD-CARGO^360 

653 



GOTO 1727 
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654 

1728 

655 

656 

1741 

657 

658 

659 

660 
661 
662 

1729 

663 

C 

664 

C 

665 

C 

666 

1726 

667 

668 

669 

670 

1742 

671 

672 

673 

674 

675 

676 

1730 

677 

C 

678 

C 

679 

C 

680 

681 

682 

683 

684 

685 

1743 

636 

687 

688 

1778 

689 

690 

691 

692 

1779 

693 

694 

695 

696 

697 

698 

699 

700 

701 

702 

1732 

703 

704 

C 

705 

c 

706 

c 

707 

708 

709 

710 

310 

711 

712 

201 

713 

714 

715 

4000 

716 

717 

718 

719 

311 


IC-IC+l 

IF( < IC.LT. 1 >.0R.{ IC.GT.40J) ) WRITE< 6 , 1 741 ) I.T.IC 
FORMAK/- FIRST RESONANCE ERROR (BAD I NOE K) ' , 1 5 , F 1 J . 5 . 1 5 ) 

IF((iC LT. 1 ) .OR.( IC.GT.40J) ) GOTO 1726 
TC< IC , I )-T 
C^K’C'.!( IC . I > -CARGO 

vpi:; ( ‘c, I )-vPHASE/i 0 j 0 . 

VF7 RA' IC, I )»VPA/1000. 

IT-0 

ENERGY SAMPLING (EVERY 8 MSEC) 

ITI-ITl+I 

IF< ITl .LT.60) GOTO 1730 
IND-INT( T*I000/6 ) + l 

IF<< IND.LT.l ).OR.{ INO.GT.850) ) WRITE! 6 . 1742 ) I.T.INO 
FORMAT!/’ TOTAL ENERGY ERROR (BAD INDEX ) ' , 1 5 , F 10. 5 , 15 ) 

IF( < IND.LT. 1 ) .OR.( INO.GT.350) ) GOTO 1730 
ENER< INO, I )-ENGY 

IFd.EQ.l) DISTAN( IND)-PHI(NL )/A 

IF< ( I .EQ.l ).AND.( S.LT.0. ) ) DISTAN( IND )»-l . *PHI ( NL )/A 

ITl-0 

CONTINUE 

PHASE BUNdiNG DETECTION ( TMI N<T<TMAX ) 

IF< (T.LT.TMIN).OR.(T.GT.TMAX) ) GOTO 1732 
IT2-IT2+1 

IF< IT2.lt. 20) GOTO 1732 
IC2-IC2>1 

IF<( IC2.lt. 1). OR. (IC2.6T.505>) WRITE ( 6 , 1 743 ) I,T,IC2 
FOr.MAT(/’ PHASE DATA ERROR (BAD I NDEX ) ' . 1 5 . F 10 . 5 , 1 5 ) 

IF< ( IC2.lt. 1 ) .OR. ( IC2.GT. 505 ) ) GOTO 1732 
CAR60-CARG/A 

IF(ABS(CARGD).LE.360. ) GOTO 1779 
IF(CARGO.LT.0. ) CARGO-CARGO-^3S0. 

IF(CARGD.GT.350. ) CARGO-CARGO-360. 

GOTO 1778 
CONTINUE 

PBCARGU( IC2, 1 )-CARGD 
PBVPHC IC2, I )-VPHASE/1000. 

PBVPA< IC2.I )-VPA/1000. 

IF( I .EQ. 1 ) TPB( IC2)-T 

IF( I .EQ. 1 ) DISTANK IC2 )-PHI(NL )/A 

IF( < I .EQ. 1 ).AND.(S.LT.0. ) ) D I STAN 1 ( IC2 ) — 1 . *PH I ( NL )/A 
IT2-0 

IF( IROONE.GT. 10) INOEXd)-! 

IFdRCONE.LE. 10) INDEXd )-0 
CONTINUE 

IF (T.GT.10) GO TO 209 

TEST FOR DETRAPPING, IF PARTICLE VEL DIFFERS FROM WAVE VEL BY 
MORE THAN SPECIFIED AMOUNT, NO INTERACTION IS ASSUMED AND ALL 
PARTICLE PARAMETERS CALC FROM ADIABATIC THEORY 
IF ( IAOIA.EQ.0) GO TO 310 

IF < ( VPA*IW0 ) .GT.0.ANO. IRDONE.GT.0.AND.(ABS( V?HASE-VPA)/VPHASE ). 

GE.0.2) GO TO 311 

IF (R.LT.(RO+l .E5) ) GO TO 201 

SO TO 110 

CONTINUE 

CALL EQCONV 

IF (IFULL.EQ.0) WRITE (6,4000) H , RPH I D . S , T . EV , EQALO 

FORMAT ( 'LANDING POI NT ’ ; 3)( , ’ H- ’ , E 1 2 . 5 , ’ '<M ’ , 3X , ’ PH I - ’ , F7 . 3 , 3X , 

•S“' .E12.5,3X. ’T-' ,F7.4,3X, 'ENERGY-' ,E8.3, ’EV ,3X, 

■EQ PITCH ANGLE- ' ,F6.3 ) 

GO TO 312 
CALL EQCONV 
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720 


721 

313 

722 

1 

723 

2 

724 

312 

725 

C 

726 

c 

727 


728 

314 

729 


730 


731 


732 


733 


734 


735 


736 


737 


733 


739 

202 

740 


741 


742 

203 

743 

204 

744 

205 

745 

206 

746 


747 


748 

3200 

7!3 

1 

1 -. ' 

2 

7o 1 


V 3 !. 


■’33 


' <. 

3300 


1 

756 


757 


758 


759 


*■60 


761 


762 

3505 

763 

3510 

764 


765 


766 


767 


763 

3502 

769 

3501 

770 

35.03 

771 

617 

772 


773 

3610 

774 


775 


776 


777 


778 


779 

3604 

780 

3603 

731 

3605 

782 

36;r6 

783 


784 


785 

3640 


IF (IFULL.EQ.0) WRITE (6,313) H , RPHI 0 , S , T. EV , EQALD 

FORMAT <• OETRAP PO I NT ' , 3X , • H- ' , E 1 2 . 5 , ' KM ' , 3X , ' PH I • ' , F 7 . 3 , 3X . 

■S-‘ .E12.5.3X. 'T-- ,F7.4,3X, 'ENERGY-' .E3.3.'EV.3X. 

'EQ PITCH ANGLE-' ,F6.3 > 

IF < IFULL .EQ. 1 > CALL DP UNC 

IF PARTICLE CROSSES EQUATOR, IRDONE PRINTED HERE IS COUNTED 
FROM EQUATOR CROSSING. 

IF< IFULL .EQ.0) WRITE(6.314) BETAD , I RDOtlE 
FORMAT!' BETA-' ,F7.2,5X, 'NO OF RESONANCES- ', 13/ ) 

ALDC( I )-EQAL 
EV0C( I )-EV 
EOTOT-EOTOT*EVO 
EFTOT-EFTOT+EV 

IF< VPAEQ.IE.''FMIN) VFNIM-VPAEQ 
IF< VPAEQ.CE.VFMAX) VFMAX-VPAEQ 
IF(EQALO.GT.ALFAHI ) ALF AH 1 -EQALD 
IF(EQALD.LT.ALFALO) ALFALO-EQALD 
JCOUNT-JCOUNT*! 

TTRACE( 1 >-T 
CONTINUE 

IF ( IFULL.EQ.0. AND. lOIFF. SO. 1 ) CALL OIFCO 

IF( < ICOKT2.EQ. 1 >. AND. ( IFULL.EQ.0) ) CALL PLOTTING 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

IF ( IFULL. EQ.l) CALL SUMARY 

1F( IFULL. EQ.! ) WR I TE ( 6 . 3200 ) VSTART , VEND , VFMI N , VFMAX 
FORMAT! ////• DISTRIBUTION FUNCTION PARAMETERS'/// 

• SVPAMIH-' ,E10.4, ■ SVPAMAX-' ,E10.4, • FVPAMI N- ' , E 10. 4 , 

' FVPAMAX-' ,E10.4// ) 

IF! IFULL .EQ. 1 ) DVPAl-OVPA'*10 
IF! IFULL. NE.l ) GOTO 3504 
WRITE! 6,3300) OVPA,DVPAl 

FORMAT!/' INITIAL VEL. B IN- ' , E 10. 4 , ' FINAL VEL. BIN-' 

.E10.4 ) 

Kl-INT! ! VFMIN-F'/PA! 1 ) )/DVPAl ) »1 
K2-INT! ! VFMAX-FVPA! 1 ) )/0VPAl )+l 
01»INT!ALMAX*2)+2 
VIRITE! 6,3510) JCOUNT 
IF! JHI .LT.35) GOTO 617 
WRITE! 6,3505) 

FORMAT!/' FINAL DISTRIBUTION !# OF PARTICLES PER CELL)') 
FORMAT!////' TOTAL NUMBER OF TRACED PARTICLES WAS-', 16//) 

DO 3501 K-K1,K2 
DO 3502 J-1,J1 
PITCH-0‘0.5-0.25 

WRITE (6, 3503) P ITCH , K , KFD 1ST! 0 , K ) 

CONTINUE 

CONTINUE 

FORMAT! F10.4, 14, ' # OF PARTICLES- ', 1 4 ) 

CONTINUE 

WRITE(6,3610) 

FORMAT!//' INITIAL DISTRIBUTION AFTER SCATTERING'/) 

DO 3603 K-1,20 

DO 36,04 J-1,01 

IF! K.GT. ! NVG+1 ) ) GOTO 3603 

PITCHl-O‘0.S-0.25 

WRITE! 6. 3605) PITCH1,K, I FOIST! ).K» 

CONTINUE 

CONTINUE 

FORMAT! F10.4, 14, ' NUMBER OF PAR^ .CLES- ' , 1 4 > 

CONTINUE 

OIFEM-EFTOT-EOTOT 
WRITE! 6,3640) DIFEN 

FORMAT!/' TOTAL ENERGY EXCHANGE ! EV )- ' , E 10. 4 ) 
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786 

C 


787 

c 

FULL OISTIUBUTION TABLE 

768 

c 


789 


IF< ( Jhl-.1L0).GT.32> GOTO 601 

790 


DO 602 J-1 .32,2 

791 

602 

ALFA( J )-J'^0.5-0.25 

792 


WRITE(6,603) <ALFA<J>, J-1 .32,2) 

793 

603 

FORMAT! IHl EQUATORIAL DISTRIBUTION FUNCTION OF 7ARTICLES > ' 

794 

1 

/■ VPARALEL <KM/SEC)' ,50X, • PITCH ANGLE (DEG)' 

795 

2 

/8X. 16F6.2/9X,32( • I’)> 

796 


IF(K1 .GT. 1 > Kl-Kl-1 

797 


IF(K2.LT.400) K2«K2*1 

798 


DO 604 K-Kl ,K2 

799 


VEL-FVPA( K )/l000. 

800 


DO 606 J-I.33 

801 

606 

KTEMP< J )-KFDIST( J.K) 

SJ2 


WRITE (6, 605) VEL,<KTEMP(J),J»1,33) 

803 

605 

FORMAT! IX. F8.0, ',33!I2,' ')) 

804 

604 

CONTINUE 

S05 

601 

CONTINUE 

806 

3504 

CONTINUE 

807 

207 

CONTINUE 

808 


IF! IGROW.NE.l ) GOTO 9003 

809 


WRITE!6,8201 > 

810 

8201 

FORMAT!/' WAVE AMPLITUDE DATA') 

811 


WRITE! 6,8081 ) XSTART , XENO . XLEN . XAMPL 

812 

8081 

FORMAT!/' START-' ,E12.4, ' END- ' , E 12 . 4 , ' LENGTH- ' ,F 10. 3, ' AMPL- 

813 

1 

• E12.4) 

814 


PHI1-PHI!NT0P )/A 

815 


PHI2-PHI!NBOT)/A 

816 


WRITE!6,9020) PHIl ,PHI2 

817 

9020 

FORMAT!/' ABSOLUTE VALUES OF STARTING NO ENDING LATITUDE ARE: 

818 

1 

' ,F10.5,3X,F10.5) 

819 


IF! IC0NT5.E0.1 ) GO TO 9002 

820 


GO TO 9003 

C21 

9002 

CONTINUE 

822 


WR I TE ! 6 . 9004 ) 

823 

9004 

FORMAT!/' WAVE AMPLITUDE DATA') 

824 


CO 9005 11-1,3000,10 

825 


WRITE (6,9006) II,Z! II ),AMPLOW! U ),AMPLHI( 11 ) 

826 

9006 

FORMAT! I5.3X,3!E12.4,3X)) 

827 

9005 

CONTINUE 

828 

9003 

CONTINUE 

829 

208 

CONTINUE 

830 


GO TO 210 

831 

209 

WRITE (6,3001) 

832 

3001 

FORMAT !///' INTEGRATION TI.1E EXCEEDS 10 SEC LIMIT') 

833 

210 

CONTINUE 

834 


STOP 

835 


END 
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836 

C 



837 



SUBROUTINE PIOTTING 

838 



COMMON/BLOC'<3/ TC< 400 , 1 2 > . CARGU ( 400 , 1 2 ) . VPHA< 400,12) . 

839 


1 

VPARA( 400, 12 ) ,Et;ER< 850, 12 > , P3CARGU( 505 , 12 ) , PBVPH( 505 , 1 2 ) 

840 


2 

,PBVPA( 505. 12 ),TMIN,TMA;<,TR< 12 ) ,TTRACE( 12 ), INDEX! 12 ) 

841 


3 

, Ml.0 ,MH I .MSTEP , TE M< 850 > . T?8( 505 ) . D I STAN! 850 > , D I STAN 1 ( 503 ) 

842 



DIMENSION SAVE! 350) .XXI < 850) ,XX2( 400 ) . XX3! 400 , 2 ) , XX4! 400 ) . 

Q43 


1 

TLO! 12) ,THI! 12) 

844 



TMAXI-0. 

84S 



DO 1 J-1 . 12 

846 



IF!TTRACE!0 ).GT,TMAXI ) TMAX I -TTRACE ! J ) 

847 


1 

CONTINUE 

648 



DO 2 0-1 . 12 

849 



DO 3 1-1,400 

850 



IF! ( INTCTC! 1 . 0 )*10000 )- INT! TR! 0 ) * 10000 ) ) . EQ . 0 ) INDEX!0 )-I 

851 


3 

CONTINUE 

852 


2 

CONTINUE 

853 



WRITE<6,6) TMAXI 

854 


6 

FORMAT! ///• PLOTTING ROUTINE STARTED'//’ MAXIMUM TRACING TIME- 

855 


1 

.F10.5) 

856 



DO 10 0-1,12 

857 


10 

WRITE!6,11) 0,TR!0),TTRACE!0) 

853 


11 

FORMAT!' PARTICLE#' ,12, • F IRST RES F 10.5 , ' ENO-',F10.5) 

859 

C 



860 

C 


FILL UP ENERGY ARRAY 

861 



DO 20 0-1,12 

862 



DO 21 1-1,850 

863 



IF!ENER! 1 ,0 ).LT.0. ) ENER! 1 ,0 )-ENER( ( 1-1 ) ,0 ) 

864 


21 

CONTINUE 

865 


20 

CONTINUE 

866 

C 



867 

C 


SUM ENERGIES FOR ALL PARTICLES 

868 

C 



869 



DO 22 1-1,850 

870 



TEMP-0. 

871 



DO 23 0-1,12 

872 


23 

TEMP-TEMP+ENER! 1,0) 

873 



ENER! 1 , 1 )-TEMP/1000. 

874 


22 

ENER! I ,2)-ENER! I . 1 )/EN£R! 1 , 1 ) 

875 



WRITE! 6,24) ENER! 1 ,1 ), ENER! 850,1 ) 

876 


24 

FORMAT!' TOTAL ENERGY DATA'//' INITAL ENERGY !EV)’,E12.4/ 

877 


1 

' FINAL ENERGY ! EV )- ' , E 1 2 . 4 ) 

878 

C 



879 

C 



880 

C 


SET UP TIME ARRAY 

881 

C 



882 



1 1- 1 NT! TMAXI* 1000/6 ) + l0 

883 



DO 60 1-1,850 

884 



IF! I.LE. II ) GOTO 61 

885 



TEN! I )-l .E36 

386 



ENER! 1 , 1 )-l.E36 

887 



ENER! I,2)-1.E36 

888 



GOTO 60 

889 


61 

TEN! I )-l*0.006 

890 


60 

CONTINUE 

891 

c 



892 

c 


PLOT ENERGY VS. TIME ! DISTANCE) 

893 

c 



894 

c 



895 

c 


DEFINE CURVE WINDOW 

896 

c 



897 



KK-l 

893 


50 

FORMAT! ' THIS IS STEP' ,13) 

899 

c 



900 



CALL AGSETF! 'GRID/LEFT. ' ,0.10) 

901 



CALL AGSETF! 'GRID/RIGHT. ' ,0.90) 

902 



CALL AGSETF! ' GR I O/BOTTOM . ' ,0. 10) 

903 



CALL AGSETF! 'GRID/TOP. ■ .0.85) 


904 

905 

906 

907 

908 

909 

910 

911 

912 

913 

914 

915 

916 

917 

918 

919 

920 

921 

922 

923 

924 

925 

926 

927 
923 

929 

930 

931 

932 

933 

934 

935 

936 

937 

938 

939 

940 

941 

942 

943 

944 

945 

946 

947 

948 

949 

950 

951 

952 

953 

954 

955 

956 

957 

958 

959 

960 

961 

962 

963 

964 

965 
965 
967 
963 
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c 

c 

c 

c 

c 

c 

c 


c 


c 

c 


63 

C 

C 


64 

C 


C 

C 

C 

C 

C 


C 

C 

C 


C 

c 


66 

67 

65 


DEFINE BACKGROUND 

CALL AGSETK 'BACKGROUND. ' ,3 ) 

TURN ON WINDOWING 

CALL AGScTK 'WINDOWING. ' ,1 > 

CALL AGSETFC 'LABEL/NAME. '. 'L' ) 

CALL AGSETK 'LlNE/MUMBEft. ' , 100) 

CALL AGSETPI 'LINE/TEXT. ', ' ENERGY (KEV)S'.l) 

CALL AGSETFI 'X/MINIMUM. ’ ,0.0) 

CALL AGSETFI 'X/MAX. ' ,TMA)n ) 

CALL AGSETK ' BOTTOM/MAOOR/TVPE . ’ ,1 ) 

CALL AGSETFI 'BOTTOM/MAJOR/BASE. ' .0.5) 

CALL AGSETK 'BOTTOM/MINOR/SPACING. ' ,4) 

CALL AGSETFI 1 1HLA3EL/MAME . , UIB ) 

CALL AGSETK 'LINE/NUMBER. ' ,-100) 

CALL AGSETPI 10HLINE/TEXT. ,1 IHTIME ISEOS.l) 

LOAD TEMP ARRAYS WITH DATA 
DO 63 1-1,850 
XXII I )-ENERI I .1 ) 

CALL EZXYITEN.XXl ,8S0,22HTOTAL ENERGY VS. TIMES) 

DO 64 1-1,850 
XXII I )-ENERI 1,2) 

CALL AGSETFI ‘LABEL/NAME.' ,'L' ) 

CALL AGSETK 'LIME/NUMBER. ' , 100) 

CALL AGSETPI ' LINE/TEXT. ' , 'E/E0S* , 1 ) 

CALL EZXYITEN.XXl ,850, 27HNORMALIZEO ENERGY VS. TIMES) 


RESET X AND REDEFINE 'NICE' 

CALL AGSETFI 'X/MAX. ■ .1.E36) 

CALL AGSETK 'X/NI . ‘ .0) 

PLOT ENERGY VERSUS LATITUDE 

CALL AGSETFI 'X/MIN. ' ,1 .E36) 

CALL AGSETFI *X/MAX. • .1.E36) 

CALL AGSETFI IIHLABEL/NAME. ,1HB) 

CALL AGSETK 'LINE/NUMBER. ' ,-100) 

CALL AGSETPI 10HLINE/TEXT. , 19HLATITUDE I DEGREES )S, 1 ) 

CALL EZXYIOISTAN.XXI .8S0.31HNORMAL1ZEO ENERGY VS. LATITUDES) 
PLOT RESONANCE DATA 

XMAXI-0. 

XMINI-1000. 

DO 65 J-1 .12 
DO 66 1-1,400 

IFITCI I,J ).GT. 1000. ) GOTO 67 
TCI I ,J >-ITCt I ,0 )-TRI J ) )“l.a00. 

THK J)-TCH I-l ),J) 

I F I TH II J ) . GT . XMAX I > XMAX I -TH II J ) 

CONTINUE 
DO 63 J-I ,12 
TLOI J )-TCI 1 ,J ) 

IFIA8SITLOIO ) ) .LT.ABSI XMINI > ) XM I N I -TLOI J ) 
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181 


969 


63 

WRITE(6,69> 

970 


69 

FORMAT!' RESCNAflCE * M 3 , ' TMI N- ' , F 1 0. 4 , ‘ TMAX- ' , F 10. 4 > 

971 



XMAXI»INT( XMAXI/10. )"10. 

972 



>!MINI-INT( XMINI/10. >«10. 

973 



IF{ABS( XMI.NI ) .'3T. 200.0) XMINI--200.0 

974 



WRITE (6, 130) XMUlI.XMAXr 

975 


130 

FORMAT!/' RESONANCE TIME WINDOW'/' TMI N- ' ,F 10. 4/ ' TMA!!»' 

976 


1 

.F10.4// ) 

977 

C 



970 

C 


SET XMIN AND XMAX 

979 

C 



980 



CALL AGSETI! 'X/NI 

981 



CALL AGSETF! 'V/MIN. ' ,0.0) 

982 



CALL AGSETF! 'V/MAX. ' ,360.0) 

983 



CALL AGSETI! 'LEFT/MAJOR/TYPE. ' , 1 ) 

984 



CALL AGSETF! 'LEFT/MAOOR/8ASE. ’ ,50.0) 

985 



CALL AGSETI! 'LEFT/MINOR/SPACING. ' ,5) 

986 



CALL AGSETF! 'X/MI .' .XMINI ) 

987 



CALL AGSETF! 'X/MA. ' ,XMAXI ) 

988 

C 



989 

c 



990 

c 


DO PHASE PLOTS 

991 

c 



992 



CALL AGSETF! ' LABEL/NAME . ' , ' L ' ) 

993 



CALL AGSETI! ' L I NE/NUMBER . ' , 100) 

994 



CALL AGSETP! 'LIME/TEXT. ', 'PHASE ! DEGREES )S ', 1 ) 

995 

c 



996 



CALL AGSETF! 'LABEL/NAME. ' ,1HB ) 

997 



CALL AGSETI! 'LINE/NUMBER. ' ,-100) 

993 



CALL AGSETP! 10HLINE/TEXT. , 12HTIME !MSEC)S,1) 

999 

c 



1000 

c 



1001 

c 


SET BOTTOM AXIS PARAMETERS 

1002 

c 



1003 



CALL AGSETI ! ' BOTTOM/MAOOR/TVPE . ' , 1 ) 

1004 



CALL AGSETF! ' BOTTOM/MAOOR/3ASE . ' ,50.0) 

1003 



CALL AGSETI! ' OOTTOM/MINOR/SPACING . ' ,4) 

1006 

c 



1007 

c 



1008 

c 



1009 

c 



1010 



DO 103 0-1 ,12 

101 1 



DO 102 I-1,C50 

1012 



XXI! I )-l .E36 

1013 


102 

ENER! I.l )-l.E36 

1014 



ICNT-1 

1015 



ENtR! ICMT, 1 )»CARGU! ICNT,0 ) 

1016 



XXl! ICNT)-TC! ICNT,0 ) 

1017 



ICNT-2 

1018 



DO 104 1-2,400 

1019 



DIFF-A8S!CARGU!!I-1 ),0)-CAR6U! 1,0)) 

1020 



IF!D1FF.LT. 180.0) GOTO 10E 

1021 



ENER! ICMT, 1 )-360.0+CARGU! 1,0) 

1022 



IF! CARGO! 1 , 0 ) . GT . CARGO! ! I-l ),0 ) ) ENER! ICNT, 1 )-CARGO 

1023 


1 

! I ,0 )-360.0 

1024 



XXI! ICNT)-TC! 1,0) 

1025 



ICNT-ICNT+1 

1026 



ENER! ICNT, 1 )-l .E36 

1027 



XXI! ICMT)-TC! I ,0 ) 

1028 



ICNT-ICNT-t-1 

1029 



ENER! ICNT, 1 )«CA.RGU! ! I-l ),0 )-360.0 

1030 



IF! CARGO! 1 ,0 > .GT.CARGO! ! I-l ),0 ) ) ENER! ICNT, 1 )-CARGU! ! I-l ),0 ) 

1031 


1 

+360.0 

1032 



XXI! ICNT)»TC! ! I-l > ,0 > 

1033 



ICNT-ICNT+1 


i 

1 


I 

i 

i 


f 



ORIGINAL PACji 353 
OF POOR QUALITY 


lfiT34 


105 

EMER{ ICNT, 1 >-CA!5.GU( I ,J ) 

li?3S 



>1X11 ICriT)-TC( I . J ) 

1036 


104 

ICNT-ICrJT»l 

1037 



CALL £2MXV< XXI ,ENER.8S0. 1 .850, 15HPHASE VS. TIMES) 

1.C33 


103 

COfjri.N'UE 

1039 



CALL AGSETFl ‘ V/MINIML'M. ' . 1 .E36 > 

1040 



CALL AGSETFl ‘ V/MAXIMUM. ‘ . 1 .E36 > 

1041 

C 



1042 

C 



1043 

C 



1044 



CALL AGSETFI ' LEFT/MAJOR/TVPE . ' ,1 .E3S) 

104S 



CALL AGSETFl ' LEFT/MAJOR/BASE . ' ,1 .E36) 

1046 



CALL AGSETFI ' L EFT/MI NOR 'SPAC ING . • .1 .E36> 

1047 

C 



1049 

C 



1049 

C 


PLOT VP AND VPA VS. TIME 

1050 

C 



1051 



CALL AGSETFI 'LABEL/NAME ’L ' ) 

1052 



CALL AGSETK ' L INE/NUMBER . ‘ . 100) 

1053 



CALL AGSETPI 'LINE/TEXT. •, ‘VELOCITY ( KM/SEC ) S ‘ , 1 ) 

1054 

C 



1055 



00 72 J-1,12 

1056 



DO 73 1-1,400 

1057 



XX2( I )-TC< I. J ) 

1058 



XX3( 1 , 1 )-VPHA( I ,J ) 

1059 



XX3< 1 ,2)-VPARA( J ,J ) 

1050 


73 

CONTINUE 

1061 



CALL EZMXV(XX2,XX3,400,2,400,18HVELOCITY VS. TIMES) 

1062 


72 

CONTINUE 

1063 

C 



1064 

C 


PLOT PHASE BUNCHING 

1065 

C 


SET X.Y AND LABELS 

1066 

C 



1067 



CALL AGSETFI • V/MI . ' ,40.0) 

1063 



CALL AGSETFI ‘Y/MA. • ,320.0) 

1069 



CALL AGSETK ‘ BOTTOM/MAJOR/TYPE . ' , 1 ) 

1070 



CALL AGSETFI ‘ BOTTOM/MAJOR/BASE . ' .0.05) 

1071 



CALL AGSETK ‘ BOTTOM/MINOR/SPACING . ' ,4) 

•y?2 

C 



773 

C 



1074 



CALL ACGETFI 'LABEL/NAME. 'L ■ ) 

1075 



CALL AGSETK 'LINE/NUMBER. 100) 

1076 



CALL AGSETPI 'LINE/TEXT. 'PHASE ( DEGREES )S‘ , 1 ) 

1077 

C 



1078 



CALL AGSETFC ‘LABEL/NAME . ' ,1HB ) 

1079 



CALL AGSETK ‘LINE/NUMBER. ■ ,-100) 

1080 



CALL AGSETPI 10HLINE/TEXT. .IIHTIME (SEC)S,1' 

1081 

C 



1082 



DO 450 J-1 ,12 

1083 



DO 401 I-l .830 

1084 



IF( I .GT.505) GOTO 402 

1085 



IF( INOEX(>) ).EQ.0) ENER( I ,J )-l .E36 

1086 



IF( INDEX! J).EQ.l ) ENERC I . J )-PBCARGU( I,J ) 

1087 



GOTO 403 

1088 


402 

ENER( 1 .0 )-l .E3G 

1089 


4«i;3 

CONTINUE 

1090 


401 

CONTINUE 

1091 


450 

CONTINUE 

1092 



DO 410 1-1,850 

1093 



IF( I.LE.505) XXK I )-TPB( I ) 

1094 



IF(I. GT.505) XXKD-1.E36 

1085 


410 

CONTINUE 

1056 



TMINI-TMIN 

1097 



TMAXI I-TMIN+0. 1 


183 
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1098 

1099 

1 100 
1101 
1 102 

1103 

1 104 

1105 

1106 

1107 

1108 

1109 

1110 
nil 
1112 

1113 

1114 

1115 

1116 

1117 

1118 

1119 

1120 
1121 
1122 

1123 

1124 

1125 

1126 

1127 

1128 

1129 

1130 

1131 

1132 

1133 

1134 

1135 

1136 

1137 

1138 

1139 

1140 

1141 

1142 

1143 

1144 

1145 

1146 

1147 

1148 

1149 

1150 

1151 

1152 

1153 

1154 

1155 

1156 

1157 

1158 

1159 

1160 
1161 


00 200 1-1,20 

IFdMAXI I .GT.TMAXI ) GOTO 201 
CALL AGSETF( ' ' .TWINI ) 

CALL AGSSTF( 'X/MAK. ' .TMAXII > 

CALL EZM>!V<XX1,£M«ER.650.12,O50.1SHPHAS2 VS. TIMES) 
TMINI-TMIMI+O. 1 
’00 TMAXII-TMAXI 1+0. 1 
201 CONTINUE 
C 

C RESET X 

C 

CALL AGSETFl ’X/MAX. ' ,TMAX) 

CALL AGSETFl 'X/MIN. ■ ,TMIN) 

CALL AGSETFl 'BOTTOM/MAJOR/TYPE. ' .1.E36) 

CALL AGSETFl ’ BOTTOM/MAJOR/BASE . ' , 1 .E36 > 

CALL ACSETFl 'BOTTOM/MINOR/SPACING. ‘ , 1 .36 > 

C 

C PLOT VPA4VPHASE VS. TIME 

C 

CALL AGSETFl • Y/Ml . * , 1 .E36) 

CALL AGSETFl 'Y/MA. • ,l.E36> 

CALL AGSETFl ‘LA8EL/NAME . • , 'L ' ) 

CALL AGSETK 'LINE/NUMBER. ‘ . 100) 

CALL AGSETPl ’LINE/TEXT. •, 'VELOCITY ( KM/SEC >S 1 ) 

CALL AGSETFl 'LABEL/NA.ME . ' . IHB > 

CALL AGSETK 'LINE/NUMBER. ' ,-100> 

CALL AGSETPl 10HLINE/TEXT. .IIHTIME (SEC)S.l) 

00 110 1-1,850 

IF< I.LE.505) XXK I )-TPB( I ) 

IF< I .GT.505) XXK 1 )-l .E36 

110 CONTINUE 

DO 111 J-1,2 
DO 112 1-1,850 

112 ENER( I . J }-l .E36 

111 CONTINUE 
C 

DO 113 J-1,12 
00 114 1-1,505 
ENER< 1 , 1 )-PBVPA( 1,3) 

114 ENER( 1 ,2 )-PBVPH< I ,3 ) 

CALL EZMXY<XX1,ENER,3S0,2,8S0,18HVELOCITV VS. TIMES) 

113 CONTINUE 
C 

C 

C PLOT VELOCITY VS. LATITUDE 

C 

DO 400 In 1.850 

IF(I.LE.505) XXKD-DISTANKI) 

IF< I. GT.505) XXK I )-l.E36 
400 CONTINUE 

CALL AGSETFl 'X/MAX. ■ ,1.E36) 

CALL AGSETF' 'X/MIN. • .1.E36) 

CALL AGSETK 'X/NI . ' ,0) 

CALL AGSETFl ' LABEL/NAME. ' ,1HB) 

CALL AGSETK ' L INE/NUMBER . ' ,-100) 

CALL AGSETPl 10HLINE/TEXT. , 19HLATITUDE ( DEGREES )$, 1 ) 

C 

DO 300 3-1,12 
DO 301 1-1,505 
ENER( I , 1 )-PBVPA( 1,3) 

301 ENERd .2 >-PBVPH( 1 .3 > 

CALL EZMXV<XXl,ENER,e50,2.850.22HVE’LOCITY VS. LATITUDES) 
300 CONTINUE 

CALL AGSETK 'X/MI . ' .-1 ) 



1162 

1163 

1164 

1165 

1166 
1167 

1138 

1139 
1 17fT 

1171 

1172 

1173 

1174 

1175 

1176 

1177 

1178 

1179 

1180 
1181 
1182 

1183 

1184 

1185 

1186 

1187 

1188 
1189 
1 1C0 

1191 

1192 

1193 

1194 

1195 

1196 

1197 

1198 

1199 

12r0 

1201 

1202 

1203 

1204 

1205 

1206 

1207 

1208 

1209 

1210 
1211 
1212 

1213 

1214 

1215 

1216 

1217 

1218 

1219 

1220 
1221 
1222 

1223 

1224 

1225 

1226 
1227 
1223 

1229 

1230 


ORJG/N/,’, 

OF POOP 




'QUALITY 


]84 


c 

c 

c 


121 


I 


1 


124 

123 


122 

101 


C 

c 

c 


1 

2 


PLOT EACH PHASE CHANGE SEPARATELY 

CALL AGSETFl ' LABEL/flAME . ’ , ' L ’ ) 

CALL AGSETK ' L I TIE /NUMSER . ' . 100) 

CALL AGSETP( 'LiriE/TEXT. 'PHASE ( DEGREES >S 1 ) 

CALL AGS£T7( 'LABEL/NAME IHB ) 

CALL AGSETK ' L I ME /NUMBER . ' .-100) 

CALL AGSETPl 10HL1N£/TE:<T. , 1 IHTIME (SEC)S.l) 

CALL AGSETF( "X/MI . • ,TMIN > 

CALL AGSETFl 'X/MA. • .TMAX) 

CALL AGSETFl • V/MIN. ‘ ,0.0) 

CALL AGSETFl 'V/MAX. ' .360.0/ 

CALL AGSETK ' LEFT/MAJOR/TVPE . * . 1 ) 

CALL AGSETFl ’ LEFT/MAOOR/BASE . ' .60.0) 

CALL AGSETK ' LEFT/MI flOR/SPAC I NG . ' .5 ) 

00 122 J>1,12 
00 121 I I« 1,850 
XXK II )*1 .E36 
ENER( II .1 )-l .E36 
ICNT-1 

ENER< ICMT. 1 >P3CARGU( ICNT.J ) 

XXK ICNT)«TP6< ICNT) 

ICNT-2 

DO 123 1-2.505 

0IFF-A3S(PLCARGU( ( I-l ) , 0 )-P3CARGU( I ,J ) ) 

IFCOIFF.LT. 180.0) GOTO 124 
EMER( ICNT, 1 )-360.*PBCARGU( I . J ) 

IF(PBCAR6U< I .0 ) .GT.PBCARGUl ( I-l ).J ) ) ENER( ICNT, I )■ 
PBCARGUl I ,J )-360.0 
XXK ICNT)-TPB( I ) 

ICNT-ICNT+l 
ENERC ICNT, 1 )-l .E36 
XXK ICNT)-TPB( I ) 

ICNT-ICNT+l 

ENER{ ICNT. 1 )-PBCARGU( ( I-l ) , J )-360.0 

IF(PBCARGU( I , J ) .GT.PBCARGUl < I-1), J)) ENER( ICNT, 1 )- 

PBCARGUl ( I-I ) ,J )+360.0 

XXK ICNT )-TPB( I-l ) 

ICNT-ICNT+l 

ENER< ICNT.l >-PBCARGU( I ,0 ) 

XXK ICNT)-TPB( I ) 

ICNT-ICNT+1 

CALL AGSETFl 'V/MI .' ,0.0) 

CALL AGSETFl ■ Y/MA. ' ,360.0) 

CALL EZMXYIXXI.ENER, 850,1, 850, 15HPHASE VS. TIMES) 

CONTINUE 

WRITE16.101 ) 

FORMATl///' ALL DONE II') 

RETURN 

END 


SUBROUTINE EQCONV 

COflMOM OVPA,EQALD,ALGRD,VPA.FVPA1400),SDIST,ALEQ.A,SVPA,FOIST1 13 
0.400) .EQAL .FPOISTl 130) ,P1 . EM , EL . RPH I ,VPE ,E .EV.KMAX .VMIN .VPMAX . 
ALIIIN.ALMAX.ALDCI 1 2 ) . R . RO . VPAEO . EPA . EVOCl 12) . IG.EPAGl 3000) 
SF-SQRTl 1 .f3.»SINlRPHI )**2) 

WPA-EM/a.^VPA^VPA 
WPE-EM/2. -VPE^VPE 
EV-lW?AHyPE )/E 
VyPEEQ»WPE/3F/l RO*EL/R )**3 
WPAEO^^WPAWPE-WPEEQ 
VPAEO-SORTl Z.n/PAEQ/EM) 

EQAL-ATANl SQRTl VyPEEQ/W?AEQ ) ) 

LC„LO-EOAL/A 

RETURN 

END 


ORIGINAL PAGK f9 
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1231 

1232 

1233 

1234 

1235 

1236 

1237 
1233 

1239 

1240 

1241 

1242 

1243 

1244 

1245 

1246 

1247 

1248 

1249 

1250 

1251 

1252 

1253 

1254 

1255 

1256 

1257 

1258 

1259 

1260 
1261 
1262 
1253 

1264 

1265 

1266 

1267 

1268 

1269 

1270 

1271 

1272 

1273 

1274 

1275 

1276 

1277 
1279 

1279 

1280 
1231 
1282 
1233 

1284 

1285 

1236 

1237 
1288 

1289 

1290 

1291 
1252 

1293 

1294 

1295 
U96 


C 

c 


c 


c 

c 


SUCROUTINE DFUNC 

COMMON DVPA.EQALO,ALGRD,VPA,FVPA( 400 ) . SD I ST , AL EQ . A , SVPA , F 0 I ST( IQ 

1 0.400) ,EQAI. ,FPDIST( 1 20 ) . P I , EM . E L . RPH I . VPE . E . E V . , VM I N . VPMAX . 

2 ALMIN.ALHA;<,ALDC< 12 > .R.RO,VPAEQ.EPA,EVOC< 12).IG.EPAG( 3000) 
COflMOM/BLOCKl/ KF D I ST( I 80 . 400 ) . I F D I 5T( 1 30 . 20 ) 

IDENTIFY SLOT FOR FVPA AHO EQALD 
J«INT( EQALO/0.5 ) + l 
ALGRO-J*0.5-0.2S 
K»INT( (VPAEQ-FVPAl 1 ) )/OVPA/ 10) + 1 
KFOISTl J ,K >«KFDIST( J .K )>l 
K1»INT( ( VPAEQ-VMIN )/DVPA ) + l 
IF<(K1.LT.1>.OR.(K1.GT.20)) GOTO 4 
IF0IST(J.K1)-IFDIST(J,KI)+1 

4 CONTINUE 

IF <ALE0.GE.5.5) SO I ST-( COSl ALEQ'A l/SVPA >**4 
IF (ALEQ.LT.S.S) SOIST"0. 

FDIST< 0 ,K >-FOIST< 0 ,K )*S0IST/12. *{ FVPA( K )/SVPA ) **2*S I H ( ALGRD*A ) 

1 /SIN(ALEQ-A)*<COS(ALEQ'*A)/COS(ALGRD’‘A> )»"3 
RETURN 
END 

SU3R0UTINE SUMARY 

COMMON OVPA, EQALD, ALGRO,VPA,FVPA< 400 > , SO I ST , ALEQ ,A , SVPA , FOI ST( 18 

1 0,400) ,EQAL,FPDIST( 180), PI , EM . E L , RP H I , VPE , E . EV , KMAX . VM I N , VPHAX . 

2 ALMIN,ALMAX,ALOC( 1 2 ) . R . RO . VP AEQ . E PA . E VDC< 12 > . IG.EPAG( 3000) 
COMMON/QLOCK2/ SFDIST( 18.0), HAS, I lAF . NVG , ALFALO . ALFAH I 

1 ,ALFA(35),JL0.3HI 

EMin-EM/2.*VMIN*VMIN 
EMAX=<EM/2.‘>VPMAX*VPMAX 
EFMIN-EM/2.*FVPA< 1 )*'FVPA( I ) 

EFMAX-EM/2 . *FVPA( KMAX )*F''PM KMAX ) 

EVMIN-EMIN/E 

EVMAX-EMAX/E 

EVFMIM-EFMIN/E 

EVFMAX-EFMAX/0 

if; IG ,NE. 1 ) WRITE(6,50)EPA 
IF<IG ,EQ. 1) UR1TE(6,51) EPAG( 1 ) 

51 F0RMAT(1H1,‘ EQ PAR E FIELD FOR GENORIN MODE- ' , 1 10. 4 , • V M-1’//) 
50 FORMAT (IHl,’ PARALLEL WAVE ELECTRIC F I ELD- ’ , E 10. 4 , • VOLT M-1'//) 
WRITE (6,6) 

5 FORMAT ( INTEGRATION RANGE'//) 

WRITE (6,5) VMIM,EMIN,EVMIN 

5 FORMAT (• MIN INITIAL VEL- ' ,E 10. 4 , ' M SEC- 1 ' , 3X , E 10. 4 , * JOULES', 

1 3X,Ei0.4,‘ EV'/) 

WRITE (6,4) VPMAX,EMAX,EVMAX 

4 FORMAT (• MAX INITIAL VEL- ' , E 10. 4 , ' M SEC- 1 ' , 3X , E 10. 1 , ' JOULES', 

1 3X,E10.4,' EV'/) 

WRITE (6,3) FVPA( 1 ),EFMIN,EVFMIN 

3 FORMAT (' MIN FINAL VEL - ' , E 10 . 4 , ' M SEC- 1 ' , 3X , E 10. 4 , ' JOULES', 

1 3X.E10.4,' EV'/) 

WRITE (6,2) FVPA(KItAX),EFMAX,EVFMAX 
2 FORMAT (' MAX FINAL VEL - ' . E 10. 4 , ' M SEC - 1 ' , 3X , E 10. 4 , ' JOULES', 

1 3X.E10.4,' EV'/> 

WRITE (6,1) ALMIN.ALMAX 

1 FORMAT (' INITIAL PITCH ANGLE RANGE - ' , 2F 6 . 2 , 3X , ' DEGREES'/) 

DO 60 J-1,180 
SFDIST(0 )-0. 

60 FPOIST(J)-0. 

DO 11 J-’ , 180 
CO 10 K-1,KMAX 

10 FPDIST( J )-2. •PI’'-FOIST( J ,IO'»FVPA( K ) '”'2’*DVPA* 10+FP D I ST( J ) 

11 CONTINUE 
K4-NVG+ 1 


OR/GjwAi r,: ? j 

OF POOR QUAUr/ 


1297 
1293 
1299 
1 320 

1301 

1302 
1 303 

1304 

1305 
1 306 

1307 

1308 

1309 

1310 

1311 

1312 

1313 

1314 
1313 

1316 

1317 
1313 

1319 

1320 

1321 
132? 

1323 

1324 

1325 
1325 

1327 

1328 

1329 

1330 

1331 

1332 

1333 

1334 

1335 

1336 

1337 

1338 

1339 

1340 


DO ’.00 J-IIAS.IIAF 
P ITr,,i3-O»0.5 + 5 .25 
DO 101 K« 1 . K4 

1F(PITCH3.GT.5.5) DIST = (C03( P I TCH 3 ’A ) / ( VM I N *0 VPA" ( K- 1 ) ) )*»4 
1F( P 1TCH3 . LE . 5 .5 > O1ST = 0. 

SrD!S''(J + ll)»2. *P I"0I ST"( VMIN + DVPA"( K - 1 > ) ■ "2'»D VPA + 3F 0 I ST J ♦ 1 1 ) 

101 CONTINUE 
100 COtITIMUE 

C PHIAL PITCH ANGLE OISTKIBUTION FUNCTION 

WRITE (6,20) 

20 FORMAT! ////• FINAL PITCH ANGLE 0 1 STR I B JT I ON • / / * PITCH ANGLE'. 5X. 

1 'NORM OIST FUNCT' ,8X . ' INIT NORM OIST FUNCT’//) 

JL0-H:T< AI.FAL0*2 ) 

JHI = INT(ALFAHI*2 ) + l 

IF! ( I IAS+1 1 ) .LT. JLO > JLO-I IAS+1 1 

IF! JH.I .LT. ! I IAF*1 1 > > JHI-I IAF+1 1 

DO 21 J»JLO,JHI 

ALGRO=>J«0.5-0.25 

21 WRITE!6,22) ALGRD . F PD 1 ST! J > . S F D I ST! J ) 

22 FORflAT!F7.2,8X.E12.4.3X.E12.4 ) 

C PRECIPITATED PARTICLE AND ENERGY FLUX 

JLOSS»iNT!3.25/0.5)*l 
PFLUX-0. 

EFLUX-0. 

DO 31 0-1 . JLOSS 
DO 30 K-1 ,KMAX 
EQAL-! J*0.5-0.25 )*A 

ACCUM-FDIST! J ,K >*FVPA! K >**2*SIN! EQAL >/C03( EQAL ) •»3»DVPA 
1 •10"J.5'-A 

PFLUX-PFLUX+ACCUM 

30 EFLU'X-£FLUX+ACCUM»0,5«EM*! FVPA! K > /COS! EQAL > )-'*2 

31 CONTINUE 

C CONVERT FLUXES TO IONOSPHERIC VALUES AT 100 KM 

PHII=ATAM!SQRT!S370.*EL/6470.-l. ) ) 

FAC -SORT! 1.4-3.*SIN!PHII )*«2)'‘EL*»3 

PFLUX=PFLUX*FAC 

EFLUX-EFLUX*FAC 

EVVLUX-EFLUX/E 

WRITE !6,40) PFLUX.EFLUX.EVFLUX 

40 FORMAT !//' PRECIPITATION FLUX- ' , E 10. 4 , ' M-2 SEC-l'//' ENERGY FLUX 
1 -'.E10.4.' JOULE M-2 SEC-1 OR ’.E10.4,' EV SEC-1') 

C FLUXES ARE NORMALIZED TO F-V---4 

RETURN 
END 


-3 


OR!G«r:AL TA ' ‘ 

187 OF POOR Q JA-iV/ 


1341 

C 



1342 



SUBROUTIfJE DIFCO 

1 343 

C 



1344 



COMMC!) DVPA.EOALO.ALGRD ,V?A.FVP/\( 4<?0 ) , 3 D I 3T , AL EQ , A . SVPA . F D I ST( IS 

1345 


1 

£/'.4,0),ECAL.FPDIST<ia0).PI,EM.:L.?'’HI.VPE.E.EV. KMAll .VMIN . VP MAX . 

13 46 


2 

ALMiri.ALMAX.ALDCl 1 2 ) . R . RO . VPAEQ . 4? ^ . E VOC ( 1 2 ) . I G , E F AG ( 30J0 ) 

1347 

C 


ALOC IS IN RADIANS, ALEQ IN DEG 

1348 



S*0 

1349 



52-0 

ISSCr 



CS-0 

1351 



CS2-0 

1352 



SS-0 

1353 



SCS-0 

1354 



SE-0. 

1355 



00 10 1-1.12 

1356 



S-S+(ALDC( I >-ALEQ"A)/12. 

1357 



S2-52h-< ALDCI I ) -AL EQ*A > »’*2/ 1 2 . 

1353 



C3-CS'MC0S<ALDC( I ) ) -COS ( ALEQ-A ) > / 1 2 . 

1359 


10 

CS2»CS2 + < COS(ALDC( I ) ) -COS ( AL EQ^A ) ) *»2/ 1 2 . 

136£f 



SD-S/A 

1361 



S2-SQRT< S2 )/A 

1362 



WRITE(6,20) S,S0,S2.CS,CS2 

1363 


20 

FORMAT!//' DEL AL-’,E10.4,' RAD OR ',F3.3,' DEG'. 3X, 'DEL AL R 

1364 


1 

MS-'. £10. 4.' DEG', 3X. 'DEL COS AL- ' ,E 10. 4 . 3X , ' DEL COS AL SO-', 

1365 


2 

E10.4 > 

1366 



DO 11 1=1,12 

1337 



SS = SLf< ALDC( I )-S-ALEQ*A)**2/12. 

1368 


1 1 

SCS-SCS-M COS(ALDC< I > )-COS< S+ALEO*A ) )**2/12. 

1369 



SS-SQRT< SS )/A 

13 70 



WRITE (6,21) SS.SCS 

1371 


21 

FORMAT (' REFE.RENCE CHANGED TO AVE SCATTERED PITCH ANGLE', 5X, 

1372 


1 

'DEL AL RMS-' ,E10.4,5X, 'DEL AL COS 5Q=',E12.4> 

1373 



DO 30 1-1,12 

1374 


30 

SE-SE+EVDC( I )/12. 

1375 



WRITE( 5,31 ) SE 

1378 


31 

FORMAT!' AVE FINAL ENERGY- ' , E 1 2 . 6 , ' EV) 

1377 



RETURN 

1378 



E -ID 
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